[Development and identifiability analysis of parent-metabolite pharmacokinetic model for risperidone and its main...

To develop a parent-metabolite pharmacokinetic model for risperidone (RIP) and its major active metabolite (9-hydroxyrisperidone) and investigate their pharmacokinetics characteristics in healthy male volunteers, twenty-two healthy volunteers were orally given a single dose of 2 mg RIP. Plasma samples were collected in the period of 96 hours and concentrations of RIP and 9-hydroxyrisperidone were measured by a validated HPLC/MS method. CYP2D6 phenotypes were identified by the T1/2 of RIP and 9-hydroxyrisperidone according to the literature. Model structure identifiability analysis was performed by the similarity transformation approach to investigate whether the unknown parameters of the proposed model could be estimated from the designed experiment. Pharmacokinetics parameters were estimated using weighted least squares method, and the final pharmacokinetics model were tested and evaluated by Monte Carlo simulation. Eighteen volunteers were phenotyped as extensive metabolizers (EM) and four volunteers were identified as intermediate metabolizers (IM). The final model included central and peripheral compartment for both parent (RIP) and metabolite (9-hydroxyrisperidone) respectively. Model structure identifiability analysis indicated that the proposed model was local identifiable. However, if the ratio of RIP converted to 9-hydroxyrisperidone was assumed to be 32% in EM, and 22% in IM, the model could be globally identifiable. The predicted time-concentration curve and AUC(0-t), C(max), T(max) of RIP and 9-hydroxyrisperidone estimated by the established model were in agreement with the observations and noncompartment analysis. Rate constant of RIP conversion to 9-hydroxyrisperidone was (0.12 +/- 0.08) h(-1) and (0.014 +/- 0.007) h(-1) for EM and IM, respectively. Elimination rate constants of RIP were (0.25 +/- 0.18) and (0.05 +/- 0.23) h(-1) for EM and IM, respectively. Model validation result showed that all parameters derived from the concentration data fitted well with the theoretical value, with mean prediction error of most PK parameter within +/- 15%. The established model well defined the disposition of RIP and 9-hydroxyrisperidone simultaneously and showed large inter-individual pharmacokinetics variation in different CYP2D6 phenotype. The model also provide a useful approach to characterize pharmacokinetics of other parent-metabolite drugs.