bioRxivPharmacology & ToxicologyPreprint
Accelerating drug development in infectious diseases using zebrafish disease models supported by pharmacokinetic pharmacodynamic modeling as new approach methodology (NAM)
Gabriel Forn‐CuníiD, Bart van LieshoutiD, Bjørn E. V. KochiD, Cristina VillellasiD, Stijn Van Asten, Dina BerlinaiD, Ellen Lanckacker, Bart StoopsiD, Rob J. VreekeniD, Dirk RoymansiD, Elke H. J. KrekelsiD, J. G. Coen van HasseltiD, Herman P. SpainkiD, Rob Christiaan van WijkiD
Abstract
The development of novel therapeutics for infectious diseases remains a global health priority. To accelerate the treatment development, innovative strategies through new approach methodology (NAM) are needed to bridge speed of in vitro with predictive power of in vivo studies, while reducing mammalian experiments. The zebrafish (Danio rerio), particularly the embryo/larva, has been established as a valuable non-mammalian in vivo model in biomedical research. We developed a standardized and streamlined workflow for the zebrafish as NAM, which consisted of 3 steps: drug selection and efficacy evaluation, internal exposure assessment, and PKPD modelling. Compounds with higher tolerated doses than minimum inhibitory concentration were selected. Drug efficacy was quantified through longitudinal individual fluorescence microscopy at baseline and 24 and 48h on treatment. Drug exposure was quantified in larval homogenates and exposure medium from 0-48h on treatment. The PKPD relationship was quantified by non-linear mixed effects modelling. For case study bedaquiline, PKPD was quantified using a one-compartment model with age-depending elimination, and an Emax concentration-response relationship on the delayed logistic bacterial growth function, with an EC50 of 26.5 g/mL and an Emax of 1.09-1.37. In the case of clarithromycin, in contrast, negligible internal exposure after waterborne treatment was observed, illustrating the risk of false negatives without internal exposure assessments. Bactericidal efficacy was confirmed by intravenous drug injections, showing a clear dose dependent antimycobacterial effect. The standardized zebrafish NAM workflow presented here facilitates the translation of drug efficacy to higher vertebrates, reducing rodent studies to confirmatory or replacing them completely, thus accelerating drug development.
- organism: zebrafish
- organism: bacteria (other)
- assay: imaging
- research focus: drug response / drug discovery
- approach: mixed model
- approach: experimental assay
- approach: statistical model
In-depth summary
Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 2 (30 Oct 2026) · CC BY 4.0
Zebrafish larva infection experiments plus pharmacokinetic-pharmacodynamic (PKPD) modelling gave bedaquiline's in vivo Emax and EC50 against M. marinum in about 2 months; waterborne clarithromycin barely entered larvae.
Zebrafish drug studies usually dose through the water and assume the internal concentration equals the medium's, which ignores absorption and elimination. This work adds measured internal exposure and PKPD modelling to zebrafish infection assays in a standardized workflow with go/no-go criteria, presented as a new approach methodology (NAM) to reduce rodent studies in anti-mycobacterial drug development.
Key results
ShownBedaquiline's exposure-response in M. marinum-infected larvae was described by delayed logistic bacterial growth with an Emax drug effect on the growth rate. The maximal effect was larger after 3 days post infection (dpi).
Growth rate constant 0.0853 h^-1. Emax, expressed as a multiple of the growth rate constant, was 1.09 up to 3 dpi and 1.37 later (increase p<0.01). EC50 26.5 µg/mL; all relative standard errors below 50%. A sigmoid Emax model fitted significantly better but its estimates were too imprecise. Fit was checked with visual predictive checks (500 simulations). · Bedaquiline: exposure-response modelling; Exposure-response modelling
ShownBedaquiline inside the larvae rose quickly, reached a dose-dependent plateau and fell in two phases after drug removal. A one-compartment model described these data, with elimination increasing at 5 days post fertilization (dpf).
45 medium measurements (0.5, 5, 50 µM) and 120 homogenates of 5 larvae each (0.5 to 50 µM), n=3 per timepoint. Absorption rate constant 0.000941 h^-1; elimination rate constant 0.089 h^-1, rising by 0.122 h^-1 at 5 dpf (p=0.025); all PK parameters had relative standard errors below 50%. · Bedaquiline: internal exposure; Bedaquiline: exposure-response modelling; Exposure-response modelling
ShownClarithromycin failed the internal-exposure criterion under waterborne dosing, yet killed bacteria when injected intravenously. It is therefore a false negative of the waterborne workflow.
Internal amounts were under 0.1 % of the external medium, although medium concentrations were stable and washing was optimized. Injecting 0.5 or 1.0 ng (vs PBS; Kruskal-Wallis with Dunn's tests) lowered M. marinum and M. avium MAC101 burden in a dose-dependent way. Internal exposure after injection was not measured, so no exposure-response was derived. · Clarithromycin did not pass workflow criteria; Figure legends
ShownBedaquiline passed the potency-versus-toxicity screen. Under waterborne treatment it reduced M. marinum burden in a dose-dependent way once differences in starting inocula were accounted for, and M. avium burden also showed a significant dose-dependent effect (Figure S2).
Minimum inhibitory concentration (MIC) 0.1 mg/L (0.18 µM), with no toxicity up to 50 µM (>25-fold margin). 60 larvae, at least n=5 per timepoint and dose (0, 0.5, 1, 5, 10, 50 µM), were imaged repeatedly. Bacterial growth slowed between the third and fourth day of infection regardless of dose. · Bedaquiline: drug efficacy
ClaimedThe authors state that the workflow gives in vivo efficacy and potency within 1 to 2 orders of magnitude of mammalian values, and that bedaquiline's exposure-response was close to that seen in rodents.
This is stated in the abstract and discussion. The text shows no rodent or mammalian parameter values and no side-by-side comparison. The workflow is said to take about 2 months, with the internal-exposure stage taking three experiments (about 1 month). · Abstract; Workflow development; Results; Discussion and conclusion
Methods, data and reproducibility
Zebrafish embryo/larva infection workflow: compare minimum inhibitory concentration with maximum tolerated dose, measure dose-ranging efficacy by repeated fluorescence imaging, quantify internal drug in whole-larva homogenates by mass spectrometry, then fit nonlinear mixed-effects PKPD models in NONMEM. Bedaquiline and clarithromycin were tested against M. marinum and M. avium.
- ABTL wild-type zebrafish embryos/larvae (up to 5 dpf) infected with mWasabi-expressing M. marinum; bedaquiline dose-ranging efficacy with repeated imaging · zebrafish (Danio rerio); Mycobacterium marinum · n = 60 larvae
- Bedaquiline internal exposure over time: incubation medium samples and larval homogenates at 0.5 to 50 µM, triplicates, plus washout phase · zebrafish (Danio rerio) · n = 45 medium measurements; 120 homogenates of 5 larvae each
- Zebrafish larvae infected with M. avium complex 101 (MAC101) and treated with bedaquiline (dose-dependent effect; Figure S2) · zebrafish; Mycobacterium avium
- Clarithromycin waterborne internal exposure, and intravenous injection (0.5 and 1.0 ng vs PBS) in M. marinum- and MAC101-infected larvae · zebrafish; M. marinum; M. avium
- Final dataset and model code (Zenodo deposit) · 10.5281/zenodo.23063496
- Methods & tools
- NONMEM v7.5.1, Pirana v23.10.1, PsN v5.3.1, R v4.3.3, RStudio v2023.12.1+402, vpc v1.2.2, QuantiFish, LC/MS (liquid chromatography-mass spectrometry), Leica M205 FA fluorescence stereomicroscope, Fish Embryo Acute Toxicity Test (OECD Test Guideline 236), EUCAST broth microdilution MIC, Kruskal-Wallis test with Dunn's post-hoc comparisons
- Reagents & environment
- NONMEM v7.5.1 (first-order conditional estimation with interaction), Pirana v23.10.1, PsN v5.3.1, R v4.3.3, RStudio v2023.12.1+402, vpc v1.2.2
- Reproducibility gaps
- Dataset and model code are on Zenodo (DOI 10.5281/zenodo.23063496, see methods.data). The protocols for MIC, infection, drug stability and homogenate quantification are in supplementary methods that are not in the text. Table 1 (parameter estimates) appears only as a caption.
Limitations (6)
- Exposure is measured as total drug per whole-larva homogenate. Converting it to a concentration assumes the drug spreads evenly, with volume of distribution equal to body volume. Blood sampling was not feasible for bedaquiline, and accumulation in tissues or at the infection site cannot be resolved. · stated by the authors
- Waterborne dosing can give false negatives when uptake is poor, as with clarithromycin. Injection avoids this but has much lower throughput, so the authors reserve it for selected promising compounds. · stated by the authors
- Larvae lack adaptive immunity and complex lung architecture. Translation to other species also depends on plasma protein binding, acute versus chronic infection, mycobacterial subpopulations and differences in body temperature. · stated by the authors
- Variability between larvae in PK could not be estimated because sampling destroys the larvae. A sigmoid Emax model was too imprecise, some dose groups had shifted prediction intervals, and the short PK window cannot capture full metabolic maturation. · stated by the authors
- Only one drug-pathogen pair (bedaquiline, M. marinum) reached PKPD modelling. The text gives no rodent or mammalian parameter values to support the stated agreement with mammalian exposure-response. · noted in this summary
- The text says clarithromycin's MIC was higher than its maximum tolerated dose (MTD), 'similar to' bedaquiline. This conflicts with bedaquiline's reported MIC/MTD margin and with the first workflow criterion, and no clarithromycin MIC or MTD values are given for these strains. · noted in this summary
Read this if you run zebrafish infection or toxicity drug screens and want to turn waterborne-dose efficacy into a measured exposure-response that can inform mammalian dose selection.
- Citations
- 1
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Pseudocapillaria tomentosa Infections in Laboratory Larval and Adult Zebrafish (Danio rerio): Development and Advances in an In Vivo Anthelmintic Drug Discovery Model
Connor Leong, Ruby Scanlon, … Michael L. Kent·Veterinary ParasitologyPeer-reviewed·Jan 20, 2026
- Long-term non-invasive drug treatments for adult zebrafish
Y. Lu, E. E. Patton·bioRxiv·Dec 21, 2021
- Screening of anti-mycobacterial compounds in a naturally infected zebrafish embryo model
J.P. Dalton, B. Uy, … S. Wiles·bioRxiv·Feb 4, 2016
- Development of zebrafish (Danio rerio) in vivo model for Borrelia burgdorferi infection
- A gelatin-based feed for individually tailored drug delivery to adult zebrafish
A. J. Ochocki, J. W. Kenney·bioRxiv·Oct 11, 2022
- Identification of potential sepsis therapeutic drugs using a zebrafish rapid screening approach
Mark W. Widder, Chance M Carbaugh, … Charlotte Anne Lanteri·LifePeer-reviewed·Dec 20, 2024