Application of Individualized PBPK Modeling of Rate Data to Evaluate the Effect of Hemodialysis on Nonrenal Clearance...

The aim of this study was to apply individualized, physiologically based pharmacokinetic modeling of ¹⁴ CO₂ production rates (iPBPK-R) measured by the erythromycin breath test to characterize the effect of hemodialysis on the function of nonrenal clearance pathways in patients with end-stage renal disease. Twelve patients previously received ¹⁴ C-erythromycin intravenously pre- and post-hemodialysis. Serial breath samples were collected after each dose over 2 hours. Eight PBPK parameters were co-estimated across periods, whereas activity of cytochrome P450 (CYP) 3A4 clearance was independently estimated for each period. Inhibition coefficients for organic anion transporting polypeptide (OATP), P-glycoprotein, and multidrug resistance-associated protein 2 activities were also estimated. Nonrenal clearance parameter estimates were explored regarding sex differences and correlations with uremic toxins and were used in hierarchical cluster analysis (HCA). Relationships between the function of nonrenal clearance pathways and uremic toxin concentrations were explored. Mean CYP 3A4 clearance increased by 2.2% post-hemodialysis. Uptake transporter activity was highly intervariable across hemodialysis. Females had 22% and 30% higher median CYP3A4 activity than males pre- and post-hemodialysis, respectively. Exploratory HCA indicated that using both CYP3A4 and OATP activity parameters at pre- and post-hemodialysis best identifies heterogeneous patients. This is the first study to use the iPBPK-R approach to simultaneously estimate multiple in vivo nonrenal elimination pathways in individual patients with kidney disease and to assess the effect of hemodialysis.