A Physiologically Based Pharmacokinetic Model of Isoniazid and Its Application in Individualizing Tuberculosis Chemotherapy

Due to its high early bactericidal activity, isoniazid (INH) plays an essential role in tuberculosis treatment. Genetic polymorphisms of N-acetyltransferase type 2 (NAT2) cause a trimodal distribution of INH pharmacokinetics in slow, intermediate, and fast acetylators. The success of INH-based chemotherapy is associated with acetylator and patient health status. Still, a standard dose recommended by the FDA is administered regardless of acetylator type or immune status, even though adverse effects occur in 5 to 33% of all patients. Slow acetylators have a higher risk of development of drug-induced toxicity, while fast acetylators and immune-deficient patients face lower treatment success rates. To mechanistically assess the trade-off between toxicity and efficacy, we developed a physiologically based pharmacokinetic (PBPK) model describing the NAT2-dependent pharmacokinetics of INH and its metabolites. We combined the PBPK model with a pharmacodynamic (PD) model of antimycobacterial drug effects in the lungs. The resulting PBPK/PD model allowed the simultaneous simulation of treatment efficacies at the site of infection and exposure to toxic metabolites in off-target organs. Subsequently, we evaluated various INH dosing regimens in NAT2-specific immunocompetent and immune-deficient virtual populations. Our results suggest the need for acetylator-specific dose adjustments for optimal treatment outcomes. A reduced dose for slow acetylators substantially lowers the exposure to toxic metabolites and thereby the risk of adverse events, while it maintains sufficient treatment efficacies. Vice versa, intermediate and fast acetylators benefit from increased INH doses and a switch to a twice-daily administration schedule. Our analysis outlines how PBPK/PD modeling may be used to design and individualize treatment regimens.