Multi-Objective Optimization for Population Pharmacokinetic Model Selection: Evaluating Non-Dominated Sorting Genetic Algorithm III Performance

ORCID

0009-0009-9641-1325 (Doncel)

Document Type

Conference Paper

Publication Date

2025

DOI

10.70534/BACR9053

Publication Title

American Conference of Pharmacometrics (ACoP2025)

Conference Name

American Conference of Pharmacometrics (ACoP2025), October 18, 2025-October 21, 2025, Aurora, Colorado

Abstract

Objectives: Non-dominated sorting genetic algorithm III (NSGA-III) is an evolutionary algorithm intended to solve multi-objective optimization (MOO) problems, particularly those with 3 or more objectives, by applying a reference point based non-dominated sorting approach [1]. This study aims to evaluate the performance of NSGA-III in the context of PopPK model selection, assessing its ability to optimize multiple competing objectives.

Methods: We use MOO to select this set of non-dominated solutions and present that set of solutions to the user for final model selection. We compare that set to the results of a traditional PopPK model selection.

Emtricitabine (FTC) and emtricitabine triphosphate (FTC-TP) PK data from the CONRAD 137 study [2] were used in this analysis. The model search space included the number of compartments for plasma FTC (1, 2, or 3), presence or absence of an absorption lag time, formation kinetics of FTC-TP from plasma FTC (linear or Michaelis-Menten), elimination kinetics of FTC-TP (linear or Michaelis-Menten), inclusion or exclusion of between-subject variability on V2, Q2, V3, and Q3, and the choice of residual error model (additive, proportional, or combined) for both plasma FTC and PBMC FTC-TP. NSGA-III algorithm was used to conduct multi-objective optimization with 3 optimization criteria: OFV (as a measure of goodness-of-fit), the total number of estimated parameters (representing model parsimony), and the prediction bias in steady-state PBMC FTC-TP trough concentration (a clinically relevant exposure metric). Inequality constraints removed crashed NONMEM runs. We ran NSGA-III with 12 partitions for 15 generations, using a population size of 92 in each generation. The final Pareto fronts from NSGA-III were compared with the model developed using the traditional stepwise method.

Results: The traditional stepwise method final model has 3 compartments for plasma FTC, with saturable formation kinetics of PBMC FTC-TP using Michaelis-Menten equation and a linear elimination from PBMC. The OFV of the final traditional selection model is 3543.015, with 18 parameters estimated. There is a 9.03% bias in steady-state PBMC FTC-TP tough concentration prediction.

The NSGA-III algorithm identified a Pareto front consisting of 30 models from the search space, with the OFV ranging from 3520.126 to 10134.822 and the total number of estimated parameters varying between 12 and 24. Among the identified Pareto fronts, the bias in steady-state PBMC FTC-TP trough concentration ranged from 0.24% to 98.09%. This Pareto front illustrates a trade-off between model complexity and predictive performance.

Conclusions: The Pareto front identified by NSGA-III provides a broader view of the optimal solution space and offers insights into the trade-offs between the competing objectives. The NSGA-III algorithm with one of the optimization criteria as bias in PBMC FTC-TP concentration was able to identify a set of models in which some of the models had less bias than traditional methods. The selection of the final model(s) from among these non-dominated models is left to the pharmacometrician as a subjective decision based on the objectives of the analysis, biological plausibility, and examination of diagnostic graphics.

Rights

© 2025 The Authors

This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0) License. This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Original Publication Citation

Yu, Y., Sale, M. E., Mazur, A., Craig, J. W., Nieforth, K., Doncel, G. F., Hendrix, C., Scott, R., & Bies, R. R. (2025). Multi-objective optimization for population pharmacokinetic model selection: Evaluating non-dominated sorting genetic algorithm III performance [Meeting Abstract]. American Conference of Pharmacometrics (ACop2025), Aurora, Colorado. https://doi.org/10.70534/BACR9053

Link to Full Text

Share

COinS