PAGE. Abstracts of the Annual Meeting of the Population Approach Group in Europe.
PAGE 27 (2018) Abstr 8637 [www.page-meeting.org/?abstract=8637]
Click to open
Poster: Methodology - Model Evaluation
Felix Held (1, 2), Carl Ekstrand (3), Edmund Hoppe (4), Marija Cvijovic (2), Johan Gabrielsson (3), and Mats Jirstrand (1)
(1) Fraunhofer-Chalmers Centre, Gothenburg, Sweden (2) Mathematical Sciences, Chalmers University of Technology and Gothenburg University, Sweden (3) Swedish University of Agricultural Sciences, Uppsala, Sweden (4) GrŁnenthal GmbH, Aachen, Germany
Introduction: Oscillating biomarker response-time courses challenge modelling of drug intervention. A periodically recurring pattern is typically seen for the stress hormone cortisol. This pattern can be captured by mechanism-based turnover models. However, analysing experimental data requires new mathematical techniques. Bayesian hierarchical modelling allows for full quantification of parameter uncertainty while also capturing the population aspects typical to nonlinear mixed effects modelling. Inter-occasion variability (IOV) is incorporated in addition to inter-individual variability (IIV).
Methods: Cortisol- and dexamethasone time courses were collected . Four different doses of dexamethasone were given (no drug and 0.1, 1, 10 µg/kg bolus + 0.07, 0.7, 7 µg/kg infusion over three hours). The pharmacokinetic/pharmacodynamic model was adapted from . Cortisol was described by a turnover model with oscillating turnover rate (average baseline kavg, amplitude α, phase-shift t0) and fractional turnover rate kout. Drug intervention was modelled with Hill-type suppression (maximum inhibition Imax, potency IC50, hill coefficient n). Dexamethasone exposure was described by a two-compartment model. The model was then extended to a population model by introduction of inter-individual and inter-occasion effects. The final model was inferred from data using a Bayesian framework with the Hamiltonian Monte Carlo algorithm in Stan . Ordinary differential equations were solved analytically for the case of constant drug exposure. The performance of the two-sample DST-protocol was studied by calculation of the specificity of the test. Specificity was predicted by Monte Carlo simulations and compared to two previously published experimental results.
Results: The proposed model described the data well. Estimated ranges for pharmacodynamic parameters were estimated as median (95% credible intervals): kavg = 12.7 (6.44, 23.5) µg L-1 h-1, α = 5.40 (1.38, 17.9) µg L-1 h-1, t0 = -3.71 (-7.54, 0.494) h, kout = 0.315 (0.221, 0.493) h-1, Imax = 0.923 (0.874, 0.965), IC50 = 0.0298 (0.00490, 0.155) µg L-1, n = 1.57 (1.03, 2.61 ). Low precision was found in the standard deviations of the random effect parameters. IIV and IOV present in the data were captured by the model. The average cortisol response level and its amplitude are suppressed with respect to magnitude and variability with increasing exposure to dexamethasone. The maximum and minimum levels of cortisol response were also suppressed by increasing exposure to dexamethasone. Mathematical expressions were derived describing cortisol oscillations with inhibition and were consistent with experimental data. Dependence of predicted specificity on drug administration time and time until measurement was observed. Different levels of variability (IIV and IOV) led to a fraction of healthy subjects with positive test results. The oscillatory behaviour of cortisol response led to an oscillatory pattern in predicted specificity.