Abstract
Optimal experimental design can be used for optimizing new pharmacokinetic (PK)–pharmacodynamic (PD) studies to increase the parameter precision. Several methods for optimizing non-linear mixed effect models has been proposed previously but the impact of optimizing other continuous design parameters, e.g. the dose, has not been investigated to a large extent. Moreover, the optimization method (sequential or simultaneous) for optimizing several continuous design parameters can have an impact on the optimal design. In the sequential approach the time and dose where optimized in sequence and in the simultaneous approach the dose and time points where optimized at the same time. To investigate the sequential approach and the simultaneous approach; three different PK–PD models where considered. In most of the cases the optimization method did change the optimal design and furthermore the precision was improved with the simultaneous approach.
Similar content being viewed by others
References
Mentré F, Mallet A, Baccar D (1997) Optimal design in random-effects regression models. Biometrika 84(2):429–442. doi:10.1093/biomet/84.2.429
Retout S, Duffull S, Mentre F (2001) Development and implementation of the population Fisher information matrix for the evaluation of population pharmacokinetic designs. Comput Methods Prog Biomed 65(2):141–151. doi:10.1016/S0169-2607(00)00117-6
Waterhouse TH et al (2005) Optimal design for model discrimination and parameter estimation for itraconazole population pharmacokinetics in cystic fibrosis patients. J Pharmacokinet Pharmacodyn 32(3–4):521–545. doi:10.1007/s10928-005-0026-2
Duffull SB, Mentre F, Aarons L (2001) Optimal design of a population pharmacodynamic experiment for ivabradine. Pharm Res 18(1):83–89. doi:10.1023/A:1011035028755
Foracchia M et al (2004) POPED, a software for optimal experiment design in population kinetics. Comput Methods Prog Biomed 74(1):29–46. doi:10.1016/S0169-2607(03)00073-7
Mould DR, Denman NG, Duffull S (2007) Using disease progression models as a tool to detect drug effect. Clin Pharmacol Ther 82(1):81–86. doi:10.1038/sj.clpt.6100228
Zhang L, Beal SL, Sheiner LB (2003) Simultaneous vs. sequential analysis for population PK/PD data I: best-case performance. J Pharmacokinet Pharmacodyn 30(6):387–404. doi:10.1023/B:JOPA.0000012998.04442.1f
Gueorguieva I et al (2007) A program for individual and population optimal design for univariate and multivariate response pharmacokinetic-pharmacodynamic models. Comput Methods Prog Biomed 86(1):51–61. doi:10.1016/j.cmpb.2007.01.004
Retout S, Mentré F (2003) Further developments of the Fisher information matrix in nonlinear mixed effects models with evaluation in population pharmacokinetics. J Biopharm Stat 13(2):209–227. doi:10.1081/BIP-120019267
Wang Y (2007) Derivation of various NONMEM estimation methods. J Pharmacokinet Pharmacodyn 34(5):575–593. doi:10.1007/s10928-007-9060-6
Beal SL, Sheiner LB (1998) NONMEM User’s Guide. Technical report. University of California, San Francisco
Brazolli C, Retout S, Mentré F (2007) Population design in nonlinear mixed effects multiple responses models: extension of PFIM and evaluation by simulation with NONMEM and MONOLIX, p 16, Abstr 1176 [www.page-meeting.org/?abstract=1176]
MONOLIX: version 2.3 (2007) http://www.monolix.org
Atkinson AC, Donev AN (1992) Optimum experimental designs. Oxford University Press, Oxford
Hooker A, Vicini P (2005) Simultaneous population optimal design for pharmacokinetic-pharmacodynamic experiments. AAPS J. 7(4):E759–E785. doi:10.1208/aapsj070476
Hashimoto Y, Sheiner LB (1991) Designs for population pharmacodynamics: value of pharmacokinetic data and population analysis. J Pharmacokinet Biopharm 19(3):333–353
Beal SL (1983) Computation of the explicit solution to the Michaelis-Menten equation. J Pharmacokinet Biopharm 11(6):641–657. doi:10.1007/BF01059062
Pronzato L, Walter E (1985) Robust experiment design via stochastic approximation. Math Biosci 75:103–120. doi:10.1016/0025-5564(85)90068-9
Tod M et al (1998) Robust optimal design for the estimation of hyperparameters in population pharmacokinetics. J Pharmacokinet Biopharm 26(6):689–716. doi:10.1023/A:1020703007613
Dodds MG, Hooker AC, Vicini P (2005) Robust population pharmacokinetic experiment design. J Pharmacokinet Pharmacodyn 32(1):33–64. doi:10.1007/s10928-005-2102-z
Dokoumetzidis A, Aarons L (2007) Bayesian optimal designs for pharmacokinetic models: sensitivity to uncertainty. J Biopharm Stat 17(5):851–867. doi:10.1080/10543400701514007
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix 1: simultaneous versus iterative sequential optimization
A repeated sequential iterative optimization on the setup presented in Fig. 1 was performed where optimization was done on the PD sample and the Dose. This was repeated numerous times to assure that the highest determinant was found. Regardless of the initial values, the simultaneous approach always found the global optimum. This was repeated with 5 different combinations of initial values for the dose and sample time (randomly chosen).
For the iterative sequential approach different initial values were also tested. After the first iteration of Time first (T|D) or Dose first (D|T) optimization initial values were updated for the next iteration to the optimal value from the previous iteration. The search was considered to have converged when the design did not change in an iteration.
The results show that regardless of the initial values for the time and an initial value of the Dose > 1.625 the local optima in the T|D approach would be found (the circles area in Fig. 1).
Similarly, for the D|T approach, regardless of the initial value of the dose and with an initial sample time of >0.91 or <0.17 the local optimum would be found. Further, as an example, with an initial value of 0.2 the global optimum will only be found after ~15 iterations.
These results show that the iterative sequential approach does not give the same results as the simultaneous approach. Further, in this simple case assuming a randomly distributed initial value over the design space, it can be shown that in ~75% of the cases with the T|D approach the probability of finding the global optima is 0. Similarly the D|T approach will in ~26% of the cases also have a probability of 0 finding the true optima while the simultaneous approach will have a probability > 0 finding the global optimum.
Appendix 2: simultaneous ED-optimal design versus sequential ED-optimal design
An ED-optimal design was tested with the simultaneous approach against the iterative sequential approach.
Model 1 was used (see Fig. 1) with an ED-uncertainty assigned to the model. The magnitude of the ED-uncertainty for the effect versus time plot is presented in Fig. 7. This effect is presented as a visual predictive check with 1000 simulations given the model (and the uncertainty). The ED-FIM surface was calculated with 90 Latin Hyper Cube samples from the distribution around the parameters with Monte Carlo Integration.
Figure 8 shows the FIM surface for this ED optimal design. Although not tested because of runtime issues, it is believed that with an initial dose >1 and using the time first approach (T|D) the probability of finding the global optimum will be 0 as in Appendix 1. In this example the advantage of the simultaneous approach is even more pronounced than for D-optimal design because of the spikier surface of the ED-optimal FIM. However, a general rule as to when an ED-optimal design will suffer more than a D-optimal design when optimizing sequentially is not presented here.
Appendix 3
See Fig. 9.
Rights and permissions
About this article
Cite this article
Nyberg, J., Karlsson, M.O. & Hooker, A.C. Simultaneous optimal experimental design on dose and sample times. J Pharmacokinet Pharmacodyn 36, 125–145 (2009). https://doi.org/10.1007/s10928-009-9114-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10928-009-9114-z