Development of a population pharmacokinetic model and optimal dosing regimen of leflunomide in Korean population

Abstract

Purpose: Leflunomide is an immunosuppressive drug indicated for the treatment of rheumatoid arthritis (RA). While the pharmacokinetics (PK) of its active metabolite A771726 reportedly show large interindividual variability, no efficient dose individualization strategy is currently available. The goal of this work was to develop a population PK model for A771726 and propose an optimal individualized dosing strategy.

Methods: A771726 plasma concentration data were collected from 50 healthy male volunteers participating in two leflunomide PK studies given a single oral dose of 40 mg. Concentrations were elevated in low body weight (WT) subjects and showed multiple peaks. Thus, A771726 PK modeling was conducted incorporating allometry scaling and enterohepatic circulation (EHC). For dose optimization, simulating a set of 1000 virtual subjects from the developed model and dividing the subjects into 5 groups with WT of 50, 60, 70, 80, 90 kg, respectively, the optimal dose was explored that achieves the drug concentration most similar to the target, which was defined as the concentration for the 70 kg subject treated with the current standard dosage regimen (the loading dose of 100 mg QD for 3 days, followed by the maintenance dose of 20 mg QD).

Results: The data were best described by a two compartment model with first order absorption incorporating EHC with the bile released into the intestine. None of the covariates tested was found to be significant other than WT used in allometry. Simulation showed that the optimal loading dose increased by 15 mg for every 10 kg increment in WT while the optimal maintenance dose was 15 and 25 mg for 50 and 90 kg groups, respectively, and the same (= 20 mg) for the others. Large concentration differences from the target observed in low and high WT groups disappeared when optimal doses were given.

Conclusions: This work demonstrates the importance of a population PK model-based dose optimization approach in maintaining drug therapeutic concentrations in leflunomide treatment.