Development of population pharmacokinetics model of isoniazid in Indonesian patients with tuberculosis

Machine Learning Approach in Dosage Individualization of Isoniazid for Tuberculosis

INTRODUCTION

Isoniazid is a first-line antituberculosis agent with high variability, which would profit from individualized dosing. Concentrations of isoniazid at 2 h (C2h), as an indicator of safety and efficacy, are important for optimizing therapy.

OBJECTIVE

The objective of this study was to establish machine learning (ML) models to predict the C2h, that can be used for establishing an individualized dosing regimen in clinical practice.

METHODS

Published population pharmacokinetic (PopPK) models for adults were searched based on PubMed and ultimately four reliable models were selected for simulating individual C2h datasets under different conditions (demographics, genotype, ethnicity, etc.). Machine learning models were trained on simulated C2h obtained from the four PopPK models. Five different algorithms were used for ML model building to predict C2h. Real-world data were used for predictive performance evaluations. Virtual trials were used to compare ML-optimized doses with PopPK model-optimized doses.

RESULTS

Categorical boosting (CatBoost) exhibited the highest prediction ability. Target C2h can be predicted using the ML model combined with the dosing regimen and three covariates (N-acetyltransferase 2 [NAT2] genotypes, weight and race [Asians and Africans]). Real-world data validation results showed that the ML model can achieve an overall prediction accuracy of 93.4%. Using the final ML model, the mean absolute prediction error value decreased by 45.7% relative to the average of PopPK models. Using the ML-optimized dosing regimen, the probability of target attainment increased by 43.7% relative to the PopPK model-optimized dosing regimens.

CONCLUSION

Machine learning models were developed with great predictive performance, which can be used to determine the individualized initial dose of isoniazid in adult patients.

Model-Informed Precision Dosing of Isoniazid: Parametric Population Pharmacokinetics Model Repository

Introduction

Isoniazid (INH) is a crucial first-line anti tuberculosis (TB) drug used in adults and children. However, various factors can alter its pharmacokinetics (PK). This article aims to establish a population pharmacokinetic (popPK) models repository of INH to facilitate clinical use.

Methods

A literature search was conducted until August 23, 2022, using PubMed, Embase, and Web of Science databases. We excluded published popPK studies that did not provide full model parameters or used a non-parametric method. Monte Carlo simulation works was based on RxODE. The popPK models repository was established using R. Non-compartment analysis was based on IQnca.

Results

Fourteen studies included in the repository, with eleven studies conducted in adults, three studies in children, one in pregnant women. Two-compartment with allometric scaling models were commonly used as structural models. NAT2 acetylator phenotype significantly affecting the apparent clearance (CL). Moreover, postmenstrual age (PMA) influenced the CL in pediatric patients. Monte Carlo simulation results showed that the geometric mean ratio (95% Confidence Interval, CI) of PK parameters in most studies were within the acceptable range (50.00-200.00%), pregnant patients showed a lower exposure. After a standard treatment strategy, there was a notable exposure reduction in the patients with the NAT2 RA or nonSA (IA/RA) phenotype, resulting in a 59.5% decrease in AUC0-24 and 83.2% decrease in Cmax (Infants), and a 49.3% reduction in AUC0-24 and 73.5% reduction in Cmax (Adults).

Discussion

Body weight and NAT2 acetylator phenotype are the most significant factors affecting the exposure of INH. PMA is a crucial factor in the pediatric population. Clinicians should consider these factors when implementing model-informed precision dosing of INH. The popPK model repository for INH will aid in optimizing treatment and enhancing patient outcomes.

Population pharmacokinetic model of rifampicin for personalized tuberculosis pharmacotherapy: Effects of SLCO1B1 polymorphisms on drug exposure

BACKGROUND

Rifampicin (RIF) exhibits high pharmacokinetic (PK) variability among individuals; a low plasma concentration might result in unfavorable treatment outcomes and drug resistance. This study evaluated the contributions of non- and genetic factors to the interindividual variability of RIF exposure, then suggested initial doses for patients with different weight bands.

METHODS

This multicenter prospective cohort study in Korea analyzed demographic and clinical data, the solute carrier organic anion transporter family member 1B1 (SLCO1B1) genotypes, and RIF concentrations. Population PK modeling and simulations were conducted using nonlinear mixed-effect modeling.

RESULTS

In total, 879 tuberculosis (TB) patients were divided into a training dataset (510 patients) and a test dataset (359 patients). A one-compartment model with allometric scaling for effect of body size best described the RIF PKs. The apparent clearance (CL/F) was 16.6% higher among patients in the SLCO1B1 rs4149056 wild-type group than among patients in variant group, significantly decreasing RIF exposure in the wild-type group. The developed model showed better predictive performance compared with previously reported models. We also suggested that patients with body weights of <40 kg, 40-55 kg, 55-70 kg, and >70 kg patients receive RIF doses of 450, 600, 750, and 1050 mg/day, respectively.

CONCLUSIONS

Total body weight and SLCO1B1 rs4149056 genotypes were the most significant covariates that affected RIF CL/F variability in Korean TB patients. We suggest initial doses of RIF based on World Health Organization weight-band classifications. The model may be implemented in treatment monitoring for TB patients.

Ionic Liquid-Based Polymer Matrices for Single and Dual Drug Delivery: Impact of Structural Topology on Characteristics and In Vitro Delivery Efficiency

Previously reported amphiphilic linear and graft copolymers, derived from the ionic liquid [2-(methacryloyloxy)ethyl]trimethylammonium chloride (TMAMA_Cl‾), along with their conjugates obtained through modification either before or after polymerization with p-aminosalicylate anions (TMAMA_PAS‾), were employed as matrices in drug delivery systems (DDSs). Based on the counterion type in TMAMA units, they were categorized into single drug systems, manifesting as ionic polymers with chloride counterions and loaded isoniazid (ISO), and dual drug systems, featuring ISO loaded in self-assembled PAS conjugates. The amphiphilic nature of these copolymers was substantiated through the determination of the critical micelle concentration (CMC), revealing an increase in values post-ion exchange (from 0.011-0.063 mg/mL to 0.027-0.181 mg/mL). The self-assembling properties were favorable for ISO encapsulation, with drug loading content (DLC) ranging between 15 and 85% in both single and dual systems. In vitro studies indicated ISO release percentages between 16 and 61% and PAS release percentages between 20 and 98%. Basic cytotoxicity assessments using the 2,5-diphenyl-2H-tetrazolium bromide (MTT) test affirmed the non-toxicity of the studied systems toward human non-tumorigenic lung epithelial cell line (BEAS-2B) cell lines, particularly in the case of dual systems bearing both ISO and PAS simultaneously. These results confirmed the effectiveness of polymeric carriers in drug delivery, demonstrating their potential for co-delivery in combination therapy.

Parametric population pharmacokinetics of isoniazid: a systematic review

INTRODUCTION

Isoniazid (INH) plays an important role in prevention and treatment of tuberculosis (TB). However, large pharmacokinetic (PK) variations are observed in patients receiving standard INH dosages. Considering the influence of PK variations on INH efficacy or adverse reactions, we reviewed the population PK studies of INH and explored significant covariates that influence INH PK.

METHODS

The PubMed and Embase databases were systematically searched from their inception to 30 January 2023. PPK studies on INH using a parametric nonlinear mixed-effect approach were included in this review. The characteristics and identified significant covariates of the included studies were summarized.

RESULTS

Twenty-one studies conducted in adults, and seven in pediatrics were included. A two-compartment model with first-order absorption and elimination was the frequently used structural model for INH. NAT2 genotype, body size, and age were identified as significant covariates affecting INH PK variation. The median clearance (CL) value in the fast metabolizers was 2.55-fold higher than that in the slow metabolizers. Infants and children had higher CL per weight values than adults with the same metabolic phenotype. In pediatric patients, CL value increased with postnatal age.

CONCLUSIONS

Compared with slow metabolizers, the daily dose of INH should be increased by 200-600 mg in fast metabolizers. To achieve effective treatment, pediatric patients need a higher dose per kilogram than adults. Further PPK studies of anti-tuberculosis drugs are needed to comprehensively understand the covariates that affect their PK characteristics and to achieve accurate dose adjustments.

Characteristics of Previous Tuberculosis Treatment History in Patients with Treatment Failure and the Impact on Acquired Drug-Resistant Tuberculosis

Tuberculosis (TB) treatment failure is a health burden, as the patient remains a source of infection and may lead to the development of multi-drug resistance (MDR). Information from cases of treatment failure that develop into MDR, which is related to a history of previous TB treatment, in accordance with the pharmacokinetic aspect, is one important thing to prevent TB treatment failure and to prevent drug resistance. This was an observational descriptive study in an acquired MDR-TB patient who had a prior history of treatment failure. A structured questionnaire was used to collect information. The questionnaire consisted of a focus on the use of TB drug formulas during the treatment period, as well as when and how to take them. This study included 171 acquired MDR-TB patients from treatment failure cases. An amount of 64 patients received the separated TB drug, and 107 patients received the fixed dose combination (FDC) TB drug. An amount of 21 (32.8%) patients receiving separated TB drug and six (5.6%) patients receiving FDC TB drug took their drug in divided doses. In addition, three (4.7%) patients receiving separated TB drug and eight (7.5%) patients receiving FDC TB drug took their drug with food. An amount of 132 out of 171 (77.2%) patients had a history of incorrect treatment that developed into MDR-TB. Education on how to take the correct medication, both the separate version and the FDC TB drug, according to the pharmacokinetic aspect, is important before starting TB treatment.

Development of population pharmacokinetics model and Bayesian estimation of rifampicin exposure in Indonesian patients with tuberculosis

BACKGROUND

Interindividual variability in the pharmacokinetics (PK) of anti-tuberculosis (TB) drugs is the leading cause of treatment failure. Herein, we evaluated the influence of demographic, clinical, and genetic factors that cause variability in RIF PK parameters in Indonesian TB patients.

METHODS

In total, 210 Indonesian patients with TB (300 plasma samples) were enrolled in this study. Clinical data, solute carrier organic anion transporter family member-1B1 (SLCO1B1) haplotypes *1a, *1b, and *15, and RIF concentrations were analyzed. The population PK model was developed using a non-linear mixed effect method.

RESULTS

A one-compartment model with allometric scaling adequately described the PK of RIF. Age and SLCO1B1 haplotype *15 were significantly associated with variability in apparent clearance (CL/F). For patients in their 40s, each 10-year increase in age was associated with a 10% decrease in CL/F (7.85 L/h). Patients with the SLCO1B1 haplotype *15 had a 24% lower CL/F compared to those with the wild-type. Visual predictive checks and non-parametric bootstrap analysis indicated good model performance.

CONCLUSION

Age and SLCO1B1 haplotype *15 were significant covariates of RIF CL/F. Geriatric patients with haplotype *15 had significantly greater exposure to RIF. The model could optimize TB pharmacotherapy through its application in therapeutic drug monitoring (clinical trial no. NCT05280886).

Advancing personalized medicine for tuberculosis through the application of immune profiling

While early and precise diagnosis is the key to eliminating tuberculosis (TB), conventional methods using culture conversion or sputum smear microscopy have failed to meet demand. This is especially true in high-epidemic developing countries and during pandemic-associated social restrictions. Suboptimal biomarkers have restricted the improvement of TB management and eradication strategies. Therefore, the research and development of new affordable and accessible methods are required. Following the emergence of many high-throughput quantification TB studies, immunomics has the advantages of directly targeting responsive immune molecules and significantly simplifying workloads. In particular, immune profiling has been demonstrated to be a versatile tool that potentially unlocks many options for application in TB management. Herein, we review the current approaches for TB control with regard to the potentials and limitations of immunomics. Multiple directions are also proposed to hopefully unleash immunomics' potential in TB research, not least in revealing representative immune biomarkers to correctly diagnose TB. The immune profiles of patients can be valuable covariates for model-informed precision dosing-based treatment monitoring, prediction of outcome, and the optimal dose prediction of anti-TB drugs.

Thermodynamic Analysis of the Solubility of Isoniazid in (PEG 200 + Water) Cosolvent Mixtures from 278.15 K to 318.15 K

The solubility of drugs in cosolvent systems of pharmaceutical interest is of great importance for understanding and optimizing a large number of processes. Here, we report the solubility of isoniazid in nine (PEG 200 + water) cosolvent mixtures at nine temperatures (278.15, 283.15, 288.15, 293.15, 298.15, 303.15, 308.15, and 318.15 K) determined by UV-vis spectrophotometry. From the solubility data, the thermodynamic solution, mixing, and transfer functions were calculated in addition to performing the enthalpy-entropy compensation analysis. The solubility of isoniazid depends on the concentration of PEG 200 (positive cosolvent effect) and temperature (endothermic process) reaching its maximum solubility in pure PEG 200 at 318.15 K and the lowest solubility in pure water at 278.15 K. The solution process is favored by the solution entropy and according to the enthalpy-entropy compensation analysis it is driven by entropy in mixtures rich in water and by enthalpy in mixtures rich in PEG 200.

Influence of N-acetyltransferase 2 (NAT2) genotype/single nucleotide polymorphisms on clearance of isoniazid in tuberculosis patients: a systematic review of population pharmacokinetic models

Purpose

Significant pharmacokinetic variabilities have been reported for isoniazid across various populations. We aimed to summarize population pharmacokinetic studies of isoniazid in tuberculosis (TB) patients with a specific focus on the influence of N-acetyltransferase 2 (NAT2) genotype/single-nucleotide polymorphism (SNP) on clearance of isoniazid.

Methods

A systematic search was conducted in PubMed and Embase for articles published in the English language from inception till February 2022 to identify population pharmacokinetic (PopPK) studies of isoniazid. Studies were included if patient population had TB and received isoniazid therapy, non-linear mixed effects modelling, and parametric approach was used for building isoniazid PopPK model and NAT2 genotype/SNP was tested as a covariate for model development.

Results

A total of 12 articles were identified from PubMed, Embase, and hand searching of articles. Isoniazid disposition was described using a two-compartment model with first-order absorption and linear elimination in most of the studies. Significant covariates influencing the pharmacokinetics of isoniazid were NAT2 genotype, body weight, lean body weight, body mass index, fat-free mass, efavirenz, formulation, CD4 cell count, and gender. Majority of studies conducted in adult TB population have reported a twofold or threefold increase in isoniazid clearance for NAT2 rapid acetylators compared to slow acetylators.

Conclusion

The variability in disposition of isoniazid can be majorly attributed to NAT2 genotype. This results in a trimodal clearance pattern with a multi-fold increase in clearance of NAT2 rapid acetylators compared to slow acetylators. Further studies exploring the generalizability/adaptability of developed PopPK models in different clinical settings are required.

Semi-Automated Therapeutic Drug Monitoring as a Pillar toward Personalized Medicine for Tuberculosis Management

Standard tuberculosis (TB) management has failed to control the growing number of drug-resistant TB cases worldwide. Therefore, innovative approaches are required to eradicate TB. Model-informed precision dosing and therapeutic drug monitoring (TDM) have become promising tools for adjusting anti-TB drug doses corresponding with individual pharmacokinetic profiles. These are crucial to improving the treatment outcome of the patients, particularly for those with complex comorbidity and a high risk of treatment failure. Despite the actual benefits of TDM at the bedside, conventional TDM encounters several hurdles related to laborious, time-consuming, and costly processes. Herein, we review the current practice of TDM and discuss the main obstacles that impede it from successful clinical implementation. Moreover, we propose a semi-automated TDM approach to further enhance precision medicine for TB management.