Physiologically-Based Pharmacokinetic-Pharmacodynamic Modeling of 1α,25-Dihydroxyvitamin D3 in Mice

1α,25-Dihydroxyvitamin D3 [1,25(OH)2D3] concentrations are regulated by renal CYP27B1 for synthesis and CYP24A1 for degradation. Published plasma and tissue 1,25(OH)2D3 concentrations and mRNA fold change expression of Cyp24a1 and Cyp27b1 following repetitive i.p. injections to C57BL/6 mice (2.5 μg × kg(-1) every 2 days for 4 doses) were fitted with a minimal and full physiologically-based pharmacokinetic-pharmacodynamic models (PBPK-PD). The minimal physiologically-based pharmacokinetic-pharmacodynamic linked model (mPBPK-PD) related Cyp24a1 mRNA fold changes to linear changes in tissue/tissue baseline 1,25(OH)2D3 concentration ratios, whereas the full physiologically-based pharmacokinetic-pharmacodynamic model (PBPK-PD) related measured tissue Cyp24a1 and Cyp27b1 fold changes to tissue 1,25(OH)2D3 concentrations with indirect response, sigmoidal maximal stimulatory effect/maximal inhibitory effect functions. Moreover, the intestinal segregated flow model (SFM) that describes a low and partial intestinal (blood/plasma) flow to enterocytes was nested within both models for comparison with the traditional model for intestine (TM) where the entire flow perfuses the intestine. Both the mPBPK(SFM)-PD and full PBPK(SFM)-PD models described the i.p. plasma and tissue 1,25(OH)2D3 concentrations and fold changes in mRNA expression significantly better than the TM counterparts with F test comparisons. The full PBPK(SFM)-PD fits showed estimates with good precision (lower percentage of coefficient of variation), and the model was more robust in predicting data from escalating i.v. doses (2, 60, and 120 pmol) and the rebound in 1,25(OH)2D3 tissue concentrations after dosing termination. The full PBPK(SFM)-PD model performed the best among the tested models for describing the complex pharmacokinetic-pharmacodynamic interplay among Cyp27b1, Cyp24a1, and 1,25(OH)2D3.