Abstract
Background: Phase 2a trials in tuberculosis typically use early bactericidal activity (EBA), the decline in sputum CFU over 14 days, as the primary end-point for testing the efficacy of drugs as monotherapy. However, the cost of phase 2a trials can range from USD 7 million to USD 19.6 million on average, while >30% of drugs fail to progress to phase 3. Better utilising pre-clinical data to predict and prioritise the most likely drugs to succeed will thus help to accelerate drug development and reduce costs. We aim to predict clinical EBA using pre-clinical in vivo pharmacokinetic (PK)-pharmacodynamic (PD) data and a model-based translational pharmacology approach.
Methods and findings: First, mouse PK, PD and clinical PK models were compiled. Second, mouse PK-PD models were built to derive an exposure–response relationship. Third, translational prediction of clinical EBA studies was performed using mouse PK-PD relationships and informed by clinical PK models and species-specific protein binding. Presence or absence of clinical efficacy was accurately predicted from the mouse model. Predicted daily decreases of CFU in the first 2 days of treatment and between day 2 and day 14 were consistent with clinical observations.
Conclusion: This platform provides an innovative solution to inform or even replace phase 2a EBA trials, to bridge the gap between mouse efficacy studies and phase 2b and phase 3 trials, and to substantially accelerate drug development.
Tweetable abstract
A translational pharmacology platform to predict early bactericidal activity https://bit.ly/3NoO89J
Introduction
Mycobacterium tuberculosis remains one of the deadliest infectious agents globally. Tuberculosis (TB) drug discovery and development activity has increased emphasis on shorter, more universal regimens to treat all TB cases independent of resistance status [1]. However, with an increasing number of new drugs and limited resources for clinical trials, further innovation of drug development is imperative to identify effective drugs and regimens more efficiently and with higher confidence [1]. A phase 2a early bactericidal activity (EBA) study is typically the first clinical evaluation of novel anti-TB drug efficacy, with the primary purpose of detecting the presence and magnitude of EBA and informing possible dose–response relationships [2]. However, the cost of phase 2a trials can range from USD 7 million to USD 19.6 million on average, while >30% of drugs fail to progress to phase 3 [3]. This highlights the challenges inherent in translating results in pre-clinical models such as in vivo mouse models or in vitro hollow-fibre systems, into successful clinical end-points and outcomes. United States Food and Drug Administration guidance for industry on drug development for pulmonary TB states that appropriate animal models can serve as an important bridge between the identification of in vitro antimycobacterial effects of an investigational drug and the initiation of clinical trials [4]. However, traditional translation of findings from pre-clinical in vivo models, by pharmacokinetic (PK) modelling and allometric scaling to identify the dosing regimen in humans that best matches the efficacious drug exposure in animals is insufficient, as it only covers exposure, but not response. Mechanistic mouse-to-human PK-pharmacodynamic (PD) models that describe the bacterial kill and PK-PD relationships are better at predicting clinical results, including the results of late-stage trials [5]. Therefore, our objective is to establish a relevant and robust model-based translational platform that can reliably link pre-clinical to clinical drug development and predict early efficacy trials for anti-TB drugs across different compound classes (figure 1). We compiled a comprehensive pre-clinical and clinical database of PK, PD and baseline bacterial growth data for 10 drugs. The drugs used to develop and validate our proposed platform consisted of a bacteriostatic antibiotic, namely ethambutol (EMB); five bactericidal antibiotics, namely isoniazid (INH), delamanid (DLM), pretomanid (PMD), linezolid (LZD), and moxifloxacin (MXF); and four sterilising antibiotics namely rifampin (RIF), rifapentine (RPT), pyrazinamide (PZA) and bedaquiline (BDQ). The translational platform in the present study intends to increase the accuracy of pre-clinical to clinical translation by enabling quantitative prediction of clinical studies from pre-clinical outputs and serves as a foundation for model-informed TB drug discovery and development.
The translational pharmacology approach to predicting early bactericidal activity (EBA) in patients. Components necessary for translation include mouse pharmacokinetics (PK)-pharmacodynamics (PD) and clinical PK (actual or scaled). The estimated relationship between drug concentration and bacterial kill is assumed to be portable after correction for protein binding and integrated with clinical PK. Using baseline bacterial burden from previous EBA trials as initial conditions, the EBA is simulated with the translational model. INH: isoniazid; RIF: rifampin; PZA: pyrazinamide; RPT: rifapentine; MXF: moxifloxacin; BDQ: bedaquiline; DLM: delamanid; PMD: pretomanid; LZD: linezolid; EMB: ethambutol; TB: tuberculosis; Q: flow between central and peripheral compartments; Vc: central compartment volume; Cp: peripheral concentration; Cc: central concentration; Vp: peripheral compartment volume; EC50: delayed concentration resulting in half of the maximal drug effect; Emax: the maximal level of drug effect; Kg: bacterial growth rate; fimm(B,t): bacteria- and time-dependent immune inhibition; Kd: bacterial death rate; PopPK: population pharmacokinetics.
Methods
Drug dataset for model building and validation
To build our model and evaluate its predictive accuracy for clinical EBA, 10 first- and second-line anti-TB drugs (BDQ, DLM, EMB, INH, LZD, MXF, PMD, PZA, RIF, RPT) were selected, for which mouse PK, mouse PD, human population PK models and human clinical EBA data were available.
Data required to assess pre-clinical drug efficacy
A large database of PK and PD data in mice for 10 TB drugs with clinical EBA data was collected (figure 2, supplementary table S1). Most experiments were performed at Johns Hopkins University (JHU; Baltimore, MD, USA), with the exception of DLM, for which PK data came both from JHU as well as from literature for one dose level [6]; EMB, for which PK data came from literature [7]; and LZD, which had data from TB Alliance (New York, NY, USA). PK experiments in BALB/c mice were dose-ranging (between two and 10 dose levels), single or multiple oral dosing for up to 8 weeks, with 29–238 observations of plasma concentration per drug. PD experiments in BALB/c mice infected through aerosol delivery were dose-ranging (between two and 15 dose levels), with treatment durations of 21–70 days, and 55–252 observations of lung CFU counts per drug. Lung CFU counts were measured by plating lung homogenates at designated time points. In the case where DLM mouse PK data showed the unexpected trend of a double peak with a single oral dose (figure 2a), we confirmed the trend with the data provider, JHU.
A rich dataset of mouse a) pharmacokinetic (PK) and b) pharmacodynamic (PD) and human c) PK and d) PD data for 10 first- and second-line tuberculosis drugs was compiled for model building. Only minimum, median and maximum doses are represented as median lines when multiple doses were present. Data points for all doses are plotted. Information on all doses is present in supplementary table S1. Mouse PK data are presented for the following doses: bedaquiline (BDQ) 12.5, 25 mg·kg−1; delamanid (DLM) 2.5, 3 mg·kg−1; ethambutol (EMB) 10, 30, 1000 mg·kg−1; isoniazid (INH) 1.56, 6.25, 25 mg·kg−1; linezolid (LZD) 5, 100, 500 mg·kg−1; moxifloxacin (MXF) 100, 200, 400 mg·kg−1; pretomanid (PMD) 6, 28.8, 486 mg·kg−1; pyrazinamide (PZA) 7, 100, 900 mg·kg−1; rifampin (RIF) 10, 15, 40 mg·kg−1; rifapentine (RPT) 5, 10, 20 mg·kg−1. All doses were given once daily unless otherwise stated. Mouse PD data are presented for the following doses: BDQ 12.5, 25, 50 mg·kg−1; DLM 3, 10, 100 mg·kg−1; EMB 100, 400, 1600 mg·kg−1; INH 0.1, 6.25, 100 mg·kg−1; LZD 100, 300, 1000 mg·kg−1; MXF 25, 50, 100 mg·kg−1; PMD 50, 100 mg·kg−1; PZA 3, 50, 900 mg·kg−1; RIF 2.5, 40, 640 mg·kg−1; RPT 5, 10, 20 mg·kg−1. All doses were given once daily, 5 days per week, unless otherwise stated. Human PK simulations from validated population PK models presented for the following doses: BDQ 25, 200, 400 mg; DLM 100, 200, 400 mg; EMB 15, 25, 50 mg·kg−1; INH 9, 75, 600 mg; LZD 600 mg once daily, 600 mg twice daily; MXF 400 mg; PMD 50, 200, 1200 mg; PZA 1500, 2000mg; RIF 600, 1350, 1950mg; RPT 300, 600, 1200 mg. All doses were given once daily, unless otherwise stated. Human phase 2a early bactericidal activity study data are presented for the following doses: BDQ 25, 200, 400 mg; DLM 100, 200, 400 mg; EMB 15, 25, 50 mg·kg−1; INH 9, 75, 600 mg; LZD 600 mg once daily, 600 mg twice daily; MXF 400 mg; PMD 50, 200, 1200 mg; PZA 200 mg; RIF 600, 1350, 1950mg; RPT 300, 600, 900, 1200 mg. All doses were given once daily, unless otherwise stated.
Mouse PK-PD model development
An integrated mouse PK-PD model, involving a PK model to describe drug exposure, a bacterial dynamics model to account for the mouse immune system and a PD model describing the combined effect of bacterial dynamics and drug effect was developed for each drug. PK data were described using one- or two-compartment models with first order absorption with or without delay, and saturable elimination when necessary. Bacterial growth dynamics without treatment were described using our previously published baseline model (supplementary equation S1) [8]. The baseline model captures the decreased rate of growth over time and attributes the decline to time- and bacteria-dependent immune control over the infection. The drug effect, measured as the log10 CFU drop independent of the immune effect over time, was incorporated using a (sigmoidal) Emax (maximal level of drug effect) relationship (supplementary equation S2). A delay effect (Kdelay) was tested where appropriate to mouse PK-PD models to establish an indirect relationship between plasma drug concentrations and drug effect at the site of action (supplementary equations S3 and S4). Detailed model development and model diagnostics can be found in the supplementary material.
Prediction of the clinical EBA
The PK-PD relationship quantified in mice was used to predict the clinical EBA. Drug concentrations in humans were simulated based on clinical population PK models (supplementary table S1) to drive the concentration–effect relationship in the clinical predictions. Where clinical population PK models were unavailable, allometric scaling from mouse PK model parameters was used [9]. Protein binding ratios between humans and mice (fu,human/mouse) were used to convert unbound plasma drug concentrations from human to mouse to translate the mouse PK-PD relationships (supplementary table S1) [10–16].
Clinical predictions for 10 drugs were simulated, with 14 unique studies at several dose levels used for validation by graphically overlaying simulated EBA from pre-clinical models with observed EBA from clinical trials. Predictions were done by simulating CFU decline in 1000 virtual patients treated with the same dose as reported in the clinical EBA study. The baseline (day 0) sputum values used were derived from the mean value for each arm reported in each study, and the variability in baseline bacterial burden between individuals used was the baseline variance among all clinical studies. The net growth and death of bacteria without treatment was assumed to be zero (supplementary equation S5). Predictions were reported as the mean±sd of the predicted time course of CFU decline. For drugs where observed data were available, the data were overlaid for visual inspection. Finally, quantitative predictions of commonly reported parameters (change from baseline to day 2 and from day 2 to day 14) were compared to the observed at various dose levels along a line of unity.
Software and statistical method
Pre-clinical and clinical PK-PD modelling was performed in NONMEM (7.4.3) through PsN (4.8.1). For LZD pre-clinical PK modelling, Monolix (5.0.0) was used. Models were developed following numerical and graphical diagnostics, assessing drop in objective function value through the likelihood ratio test and parameter precision, as well as goodness-of-fit plots and visual predictive checks, respectively, in addition to pharmacological relevance. Data transformation and graphical output were performed in R (4.1.3) through the RStudio (2022.02.3) interface using the xpose4 and tidyverse packages.
Results
Large pre-clinical and clinical PK and PD database of 10 TB drugs
We collated a rich longitudinal dataset of mouse PK (plasma concentrations, 1220 data points) and PD data (lung CFU counts, 1550 data points), as well as human population PK models and human PD data (sputum CFU counts) (supplementary table S1). PD experiments were done mostly in mouse infection models infected via aerosol with an inoculum size ≥3.5 log10 CFU·mL−1 and incubation periods of 13–17 days, prior to the start of treatment. Exceptions were EMB and LZD, which had incubation periods of 7 and 5 days, respectively, but had a similar inoculation size of >3.5 log10 CFU·mL−1, and RPT, which had an incubation period of 41 days, but an inoculation size of <3.5 log10 CFU·mL−1.
Human PK data were simulated using published models from literature (supplementary table S1 and figure 2c). Human PD data with a total of 287 human sputum CFU datapoints originating from phase 2a trials across 14 different studies ranging from 2 to 14 days were used to validate our phase 2a EBA predictions.
Pre-clinical PK and PK-PD models adequately described mouse data
The final PK and PK-PD model parameter estimates are shown in tables 1 and 2. A two-compartment model with saturated clearance described via the Michaelis–Menten equation best described the mouse plasma concentration data for INH, LZD, PMD, PZA and RIF. BDQ, EMB and MXF were best described using two-compartment models with linear elimination; RPT by a one-compartment model with saturated elimination; and DLM by a one-compartment model with linear elimination. Visual predictive checks of the final model for both mouse PK and PK-PD data showed good fits (supplementary figures S1 and S2). The exposure–response relationships for each drug in mouse infection models are summarised in table 1 and supplementary figure S3 and aligned with clinical knowledge of the efficacy of each drug.
Parameter estimates of final pharmacokinetic (PK) models for 10 tuberculosis drugs in mouse studies
Parameter estimates of final pharmacokinetic (PK)-pharmacodynamic (PD) models for 10 tuberculosis drugs in mouse studies
Clinical EBA was well predicted by the translational platform
The translational platform predicted clinical EBA in TB patients receiving monotherapy with the 10 drugs, as shown in figure 3. Our predictions overlapped well with the observed data across multiple doses and time points for most of the drugs. BDQ and LZD had slight over-predictions at the later time, and RPT showed activity up to 5 days after a single dose, whereas our model predicted limited declines in CFU.
Translational (mouse to human) pharmacokinetic (PK)-pharmacodynamic (PD) model predicts clinical early bactericidal activity (EBA) trial results well. Median (95% CI) of 1000 simulations from the translational model overlap with observed EBA data from clinical trials. BDQ: bedaquiline; DLM: delamanid; EMB: ethambutol; INH: isoniazid; LZD: linezolid; MXF: moxifloxacin; RIF: rifampin; RPT: rifapentine; PZA: pyrazinamide; PMD: pretomanid.
Agreement between predicted and observed quantitative change in CFU is shown in figure 4 as a correlation plot for EBA at time intervals of 0–2 days and 2–14 days. Most predictions for 0–2 days fell within 0.25 log10 CFU·mL−1·day−1 of the observed EBA as indicated by the line of unity and corresponding dotted lines. Predictions for 2–14 days were even closer to observed. Predictions were overall consistent with the observed data in the clinical EBA studies for all 10 drugs, except for RPT, where activity was under-predicted.
Model-based prediction of daily change in log10 CFU·mL−1 correlates well with clinically observed daily change in log10 CFU·mL−1 for 10 tuberculosis drugs at multiple dose levels of monotherapy between a) days 0–2 and b) days 2–14. For some drugs, day 14 data were not available. Data are presented as line of unity±0.25. BDQ: bedaquiline; DLM: delamanid; EMB: ethambutol; INH: isoniazid; LZD: linezolid; MXF: moxifloxacin; PMD: pretomanid; PZA: pyrazinamide; RIF: rifampin; RPT: rifapentine. #: regimen contained a loading dose.
Discussion
We established a mouse-to-human translational platform by integrating a bacterial dynamics model, mouse PK-PD relationships, clinical PK models and species-specific drug plasma protein binding, and validated the platform with clinical EBA data (figure 1). The changes in sputum CFU counts over the first 2 days and from day 2 to day 14 in TB patients receiving monotherapy with each of 10 TB drugs in 14 clinical EBA studies spanning more than four decades were predicted successfully (figures 3 and 4). Compared to the participants enrolled in more recent EBA studies (2007 to 2015) [17–23] at the same site, the participants enrolled between 1992 and 2005 [24–28] had more severe disease and therefore higher baseline CFU counts in their sputum samples (mean baseline 6.9 log10 CFU·mL−1). However, the predictive accuracy of our model was robust despite this large variation in baseline bacterial burden. For example, RIF had a good overlap of predicted and observed EBA (figure 3) despite the study being conducted in 2015 with the lowest median baseline of 4.58 log10 CFU·mL−1 [18].
A key component to our model accuracy was the addition of the bacterial dynamics model. Mouse and human immune activation against TB infection differ significantly; therefore, the underlying baseline of bacterial dynamics will differ. Subtracting the mouse immune effect on bacterial decline more accurately estimates the drug contribution to CFU decline. Without such consideration, the clinical CFU decline was over-predicted (supplementary figure S4). Despite inherent differences between species in terms of drug pharmacokinetics, sampling (whole lung homogenate versus sputum), and infecting bacterial strain, the relationship between drug effect on bacteria and the concentration to achieve the effect appear, based on this analysis, to be portable between mice and patients. In addition, although the mouse strain used in the studies (BALB/c) models intracellular bacteria, but not extracellular bacteria in caseous lesions [29], the PK-PD relationships observed in this model, when derived in comparison to the baseline bacterial dynamics, appear to accurately reflect those observed in EBA studies. Other approaches or more information may be needed to fully account for drug exposures at the site of infection in cavities or other caseous lesions or any PK-PD relationships unique to those microenvironments.
Murine TB models are routinely and often exclusively used as in vivo efficacy models in pre-clinical TB drug development [1]. As the inoculum size and incubation period for bacterial infection in the lung prior to treatment can affect drug response [8], we standardised our inclusion criteria to experiments using the most common design with the incubation duration of 13–17 days and inoculum size to ≥3.5 log10 CFU·mL−1. Incubation durations outside this range were considered only when data were otherwise not available, which was the case for EMB, LZD and RPT.
Clinical EBA studies are the only acceptable way to evaluate a drug as monotherapy in TB patients despite their limitations in predicting long-term efficacy. In addition to detecting the presence of an EBA response, the trial can inform on the dose–response curve (e.g. INH and RIF), which could be used in dose selection for future trials [18, 30, 31]. We have shown here that our translational platform can adequately predict these EBAs for different doses. With limited resources, this costly clinical study can be designed more efficiently or avoided altogether by using our approach to predict a reliable result regarding clinical dose–response effects, and to provide useful information about dose and/or drug candidate selection for further clinical development. This scenario is well exemplified by the nitroimidazole, PMD. PMD has a dose response at doses up to 192 mg·kg−1 in mice, which, following the conventional allometric scaling method, approximates 1500 mg in humans. However, such translation is problematic, as the clinical observations from two human EBA trials demonstrated no dose response above 200 mg in human EBA. Using our translational platform, we found that the drug effect of PMD reaches plateau after 200 mg, which is consistent with clinical observations (supplementary figure S3). Therefore, our translational platform could serve as a powerful tool for, but not limited to, better dose selection for clinical trials design. By better informing dose selection, the translational modelling platform may reduce the time and effort spent in early clinical development, and therefore, accelerate progress to trials that are more informative of long-term outcomes.
Phase 2a trials also gather information on initial safety and tolerability of compounds of interest, which our translational framework is not designed to predict. However, when anti-TB drug development progresses directly to regimen efficacy studies, these safety data can have already been captured during phase 1 healthy volunteer monotherapy studies, especially the multiple ascending dose study which is a dose-ranging study for up to 4 weeks. The trial population of healthy volunteers is different from the patient population, but safety signals are more pronounced in this healthy population that is more sensitive to adverse events and less burdened by other symptoms. For example, treatment of healthy volunteers with rifapentine resulted in safety signals, but has been proven safe in patients with TB [32]. Furthermore, we can also get guidance on the optimal efficacious dose versus the safety range using the dose–response curve by overlaying both exposures (supplementary figure S3). All our studied drugs have been approved previously and are used clinically. Having all clinical dose exposures lower than safety limits were thus unsurprising and reassuring. Such dose–response curves as visualisation would be helpful for determining the dose of new TB drugs too, as they provide a measure of both efficacy and safety.
Of the clinical EBA studies included in our analysis, the RPT EBA trial was the only one in which EBA was assessed for multiple days after a single dose. Our human population PK model indicated that RPT was mostly cleared from the body 2 days after a single dose, but the trial results indicated that RPT was still exerting an effect on bacterial load between 2 and 5 days post-dose. It is possible that RPT has a post-antibiotic effect that was not sufficiently captured by the model. The model over-predicted the EBA of BDQ. However, in the model, the active metabolite, BDQ-M2, was not considered. In mice, M2 is estimated to contribute ∼50% of the drug effect. One possible reason for the over-prediction are the parent-to-metabolite ratios between species that differ, where mice have higher M2-to-BDQ ratios than humans [33]. Future studies can account for these differences.
Our translational framework has been developed to predict clinical EBA in the typical adult patient population participating in EBA trials. Heterogeneity in the patient population is an important consideration in drug development and individualised medicine [34]. This is both from the perspective of representing patients in easy- or hard-to-treat phenotypes, as we have observed in our clinical projects as a result from different risk factors (baseline bacterial burden, disease phenotype such as cavitation, gender, comorbidities or co-medications), as well as from a diversity perspective to study pharmacology in under-represented individuals [35, 36]. Certain risk factors such as baseline bacterial burden are easily implementable in our clinical simulations, as well as, for example, the influence of comorbidities or co-medications on the pharmacokinetics. However, pre-clinical (mouse) models are traditionally more homogeneous to reduce noise in the data and be more sensitive to detect pharmacological signals. At the same time, the EBA trials with relatively small sample sizes (<15 per arm) will also not reflect clinical heterogeneity, and risk factors other than those described have limited relevance for the prediction of the EBA (e.g. cavitation). Similar considerations are applicable to the prediction of EBA in children. Throughout their development, infants and children show changes in their pharmacology that are well established and can be incorporated in our quantitative model-based framework. Pharmacodynamically, the bacterial dynamics and the disease phenotype differences can be modelled based on different animal disease models. Children with TB aged <1 year have limited immunity, which can be approximated through the immunocompromised athymic (nude) mouse model, while dose for those aged >1 year without lesion phenotypes can be approximated through the BALB/c mouse model. One limitation of these pre-clinical models is their reliance of bacterial load measurement, which is complicated in the paediatric population in the context of sputum collection. Alternative models are being developed which are part of future collaborative work in our group. Pharmacokinetically, the development of the metabolic pathways responsible for the elimination of anti-TB drugs can be incorporated through the use of maturation functions. As a result, paediatric dosing can be projected that will reach similar exposure as in adults given a chosen dose, based on the understanding of the maturation of the relevant elimination pathways and the adult pharmacokinetics.
Building on our translational framework, we aim to predict the efficacy of combination regimens of TB drugs in long-term TB clinical outcomes for phase 2b and 3 from pre-clinical mouse data. Being able to better understand the time to stable culture conversion and relapse 6 months post-treatment will better help us prioritise sterilising regimens. We hope to achieve this by including the characterisation of PK-PD relationships in combination regimens by accounting for PK-PD drug–drug interactions, as well as characterising lesion-specific PK-PD relationships. Technically, the bacterial dynamic parameters of the translational tool will be re-evaluated and possibly updated through Bayesian methods based on untreated control data of ongoing experiments with novel anti-TB drugs, benefiting from a larger data collection while keeping the structure of the translational tool. Clinical TB disease (e.g. caseation necrosis and cavitation) will be represented in the translational platform to include infection and efficacy data in animal models of TB with more human-like necrotic lesions, such as C3HeB/FeJ mice and New Zealand white rabbits [37]. This will allow us to have a comprehensive platform that informs us not only of monotherapy EBA, but also combination drug regimen efficacy in providing a stable cure and prevention of long-term unfavourable outcomes.
In summary, we established a foundation for translating the results from mouse efficacy models to clinical EBA studies through establishing quantitative relationships involving mouse PK and PD models, as well as drug dose response in vivo. In the future, our platform will be expanded to include combination regimens and longer durations of treatment by accounting for PK-PD drug–drug interactions and necrotic lesion penetration. This platform is an innovation to accelerate TB drug development and a good example of model-informed drug discovery and development.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material ERJ-00165-2023.Supplement
Shareable PDF
Supplementary Material
This one-page PDF can be shared freely online.
Shareable PDF ERJ-00165-2023.Shareable
Acknowledgements
We acknowledge TB Alliance for support and generously sharing the in-house data of anti-TB drugs linezolid and pretomanid.
Footnotes
Author contributions: The manuscript was written by J.P. Ernest, J.J.N. Goh, N. Strydom, Q. Wang, R.C. van Wijk and N. Zhang and commented on by all authors. J.P. Ernest, J.J.N. Goh, N. Strydom, Q. Wang, R.C. van Wijk and N. Zhang contributed to data collection, model development, data and model management, and code review. Contributions to data collection and model development were as follows. J.P. Ernest: BDQ and RPT; J.J.N. Goh: EMB; N. Strydom: LZD; Q. Wang: DLM, PMD and MXF; N. Zhang: PZA, RIF and INH. J.J.N. Goh and R.C. van Wijk revised the manuscript, including generation of figures and tables. J.P. Ernest and R.C. van Wijk carried out code review for all drugs. A. Deitchman worked on data collection, human PK model development, simulation, and preliminary model development. E. Nuermberger provided pre-clinical data used in our current study, provided substantial scientific context, and edited the manuscript. R.M. Savic supervised the whole research.
This article has an editorial commentary: https://doi.org/10.1183/13993003.01108-2023
Conflict of interest: R.M. Savic and E. Nuermberger report support for the present work from TB Alliance and National Institutes of Health. In addition, E. Nuermberger also reports advisory board participation with Janssen, outside the submitted work. The remaining authors have no potential conflicts of interest to disclose.
Support statement: This work was supported by NIH grant R01AI-111992. Funding information for this article has been deposited with the Crossref Funder Registry.
- Received January 25, 2023.
- Accepted June 2, 2023.
- Copyright ©The authors 2023.
This version is distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0. For commercial reproduction rights and permissions contact permissions{at}ersnet.org