ABSTRACT
Conventional pharmacokinetic (PK) bioequivalence (BE) studies aim to compare the rate and extent of drug absorption from a test (T) and reference (R) product using non-compartmental analysis (NCA) and the two one-sided test (TOST). Recently published regulatory guidance recommends alternative model-based (MB) approaches for BE assessment when NCA is challenging, as for long-acting injectables and products which require sparse PK sampling. However, our previous research on MB-TOST approaches showed that model misspecification can lead to inflated type I error. The objective of this research was to compare the performance of model selection (MS) on R product arm data and model averaging (MA) from a pool of candidate structural PK models in MBBE studies with sparse sampling. Our simulation study was inspired by a real case BE study using a two-way crossover design. PK data were simulated using three structural models under the null hypothesis and one model under the alternative hypothesis. MB-TOST was applied either using each of the five candidate models or following MS and MA with or without the simulated model in the pool. Assuming T and R have the same PK model, our simulation shows that following MS and MA, MB-TOST controls type I error rates at or below 0.05 and attains similar or even higher power than when using the simulated model. Thus, we propose to use MS prior to MB-TOST for BE studies with sparse PK sampling and to consider MA when candidate models have similar Akaike information criterion.
Subject(s)
Computer Simulation , Cross-Over Studies , Models, Statistical , Therapeutic Equivalency , Humans , PharmacokineticsABSTRACT
Covariate analysis in population pharmacokinetics is key for adjusting doses for patients. The main objective of this work was to compare the adequacy of various modeling approaches on covariate clinical relevance decision-making. The full model, stepwise covariate model (SCM) and SCM+ PsN algorithms were compared in a clinical trial simulation of a 383-patient population pharmacokinetic study mixing rich and sparse designs. A one-compartment model with first-order absorption was used. A base model including a body weight effect on CL/F and V/F and a covariate model including 4 additional covariates-parameters relationships were simulated. As for forest plots, ratios between covariates at a specific value and that of a typical individual were calculated with their 90% confidence interval (CI90) using standard errors. Covariates on CL, V and KA were considered relevant if their CI90 fell completely outside the reference area [0.8-1.2]. All approaches provided unbiased covariate ratio estimates. For covariates with a simulated effect, the 3 approaches correctly identify their clinical relevance. However, significant covariates were missed in up to 15% of cases with SCM/SCM+. For covariate with no simulated effects, the full model mainly identified them as non-relevant or with insufficient information while SCM/SCM+ mainly did not select them. SCM/SCM+ assume that non-selected covariates are non-relevant when it could be due to insufficient information, whereas the full model does not make this assumption and is faster. This study must be extended to other methods and completed by a more complex high-dimensional simulation framework.