Model-based meta-analysis for development of a population-pharmacokinetic (PPK) model for Vitamin D3 and its 25OHD3 metabolite using both individual and arm-level data

Alanna S. Ocampo-Pelland, Marc R. Gastonguay, Jonathan F. French, Matthew M. Riggs

Clinical studies investigating relationships between D3 and 25OHD3 vary in dosing regimen, assays, demographics, and control of exogenous D3. This leads to uncertain and conflicting exposure-related associations with D3 and 25OHD3. To elucidate this parent-metabolite system, a PPK model was developed to predict mean D3 and 25OHD3 exposure from varied doses and administration routes. Sources of exposure variability related to metabolite baseline, weight, and assay type were explored. Specific search criteria were used in PUBMED to identify public source PK data pertaining to D3 and 25OHD3 in healthy or osteoporotic populations. Overall 57 studies representing 5395 individuals were selected, including 25 individual-level profiles and treatment-arm data. IV, oral, single and multiple dose data were used, with D3 and 25OHD3 dosing. A nonlinear mixed effects model was developed to simultaneously model PK dispositions of D3 and 25OHD3 (NONMEM v7.2), which were described by 2-compartment models with nonlinear and linear clearances, respectively. Proportional and additive assay variances were included on the 25OHD3 prediction. Unit-level random effects were weighted by treatment-arm size. D3 model estimates, relative to bioavailability were: maximum rate of metabolism ([Formula: see text], 1.62 nmol/h), Michaelis-Menten constant ([Formula: see text], 6.39 nmol/L), central volume of distribution ([Formula: see text], 15.5 L), intercompartmental clearance ([Formula: see text], 0.185 L/h), peripheral volume of distribution ([Formula: see text], 2333 L/h), and baseline concentration ([Formula: see text], 3.75 nmol/L). For 25OHD3 ([Formula: see text] = metabolite): [Formula: see text] = 0.0153 L/h, [Formula: see text] = 4.35 L, [Formula: see text] = 6.87 L, [Formula: see text] = 0.0507 L/h. Simulations of 25OHD3 concentration indicated an inverse relationship between 25OHD3 baseline and response, as well as a less than proportional 25OHD3 response. Estimation of assay bias parameters suggested that HPLC-MS and RIA produced similar measurement results, whereas CPBA and CHEMI are over-predictive of 25OHD3 concentration, relative to HPLC-MS.