Modeling and Simulation for Estimating the Influence of Renal Dysfunction on the Hypouricemic Effect of Febuxostat in Hyperuricemic Patients Due to Overproduction or Underexcretion of Uric Acid

Whether renal dysfunction influences the hypouricemic effect of febuxostat, a xanthine oxidase (XO) inhibitor, in patients with hyperuricemia due to overproduction or underexcretion of uric acid (UA) remains unclear. We aimed to address this question with a modeling and simulation approach. The pharmacokinetics (PK) of febuxostat were analyzed using data from the literature. A kinetic model of UA was retrieved from a previous human study. Renal UA clearance was estimated as a function of creatinine clearance (CLcr) but non-renal UA clearance was assumed constant. A reversible inhibition model for bovine XO was adopted. Integrating these kinetic formulas, we developed a PK-pharmacodynamic (PK-PD) model for estimating the time course of the hypouricemic effect of febuxostat as a function of baseline UA level, febuxostat dose, treatment duration, body weight, and CLcr. Using the Monte Carlo simulation method, we examined the performance of the model by comparing predicted UA levels with those reported in the literature. We also modified the models for application to hyperuricemia due to UA overproduction or underexcretion. Thirty-nine data sets comprising 735 volunteers or patients were retrieved from the literature. A good correlation was observed between the hypouricemic effects of febuxostat estimated by our PK-PD model and those reported in the articles (observed) (r=0.89, p