Development of a genetic algorithm and NONMEM workbench for automating and improving population pharmacokinetic/pharmacodynamic model selection

Abstract

The current approach to selection of a population PK/PD model is inherently flawed as it fails to account for interactions between structural, covariate, and statistical parameters. Further, the current approach requires significant manual and redundant model modifications that heavily lend themselves to automation. Within the discipline of numerical optimization it falls into the “local search” category. Genetic algorithms are a class of algorithms inspired by the mathematics of evolution. GAs are general, powerful, robust algorithms and can be used to find global optimal solutions for difficult problems even in the presence of non-differentiable functions, as is the case in the discrete nature of including/excluding model components in search of the best performing mixed-effects PK/PD model. A genetic algorithm implemented in an R-based NONMEM workbench for identification of near optimal models is presented. In addition to the GA capabilities, the workbench supports modeling efforts by: (1) Organizing and displaying models in tabular format, allowing the user to sort, filter, edit, create, and delete models seamlessly, (2) displaying run results, parameter estimates and precisions, (3) integrating xpose4 and PsN to facilitate generation of model diagnostic plots and run PsN scripts, (4) running regression models between post-hoc parameter estimates and covariates. This approach will further facilitate the scientist to shift efforts to focus on model evaluation, hypotheses generation, and interpretation and applications of resulting models.