The aim of the study is to develop a population pharmacokinetic (PPK) model, of Octreotide long acting repeatable (LAR) formulation in healthy volunteers, which describes the highly variable, multiple peak absorption pattern of the pharmacokinetics of the drug, in individual and population levels. An empirical absorption model, coupled with a one-compartment distribution model with linear elimination was found to describe the data well. Absorption was modelled as a weighted sum of a first order and three transit compartment absorption processes, with delays and appropriately constrained model parameters. Identifiability analysis verified that all twelve parameters of the structural model are identifiable. A machine learning method, i.e., cluster analysis, was performed as pre-processing of the PK profiles, to define subpopulations, before PPK modelling. It revealed that 13% of the patients deviated considerably from the typical absorption pattern and allowed better characterization of the observed heterogeneity and variability of the study, while the approach may have wider applicability in building PPK models. The final model was evaluated by goodness of fit plots, Visual Predictive Check plots and bootstrap. The present model is the first to describe the multiple-peak absorption pattern observed after octreotide LAR administration and may be useful to provide insights and validate hypotheses regarding release from PLGA-based formulations.
Keywords: long acting injectables; machine learning; octreotide; population pharmacokinetics.