Volume 4 Supplement 1
The ISIBM International Joint Conferences on Bioinformatics, Systems Biology and Intelligent Computing (IJCBS)
Noncompartment model to compartment model pharmacokinetics transformation metaanalysis – a multivariate nonlinear mixed model
 Zhiping Wang†^{1},
 Seongho Kim†^{1},
 Sara K Quinney^{1},
 Jihao Zhou^{2} and
 Lang Li^{1}Email author
DOI: 10.1186/175205094S1S8
© Li et al; licensee BioMed Central Ltd. 2010
Published: 28 May 2010
Abstract
Background
To fulfill the model based drug development, the very first step is usually a model establishment from published literatures. Pharmacokinetics model is the central piece of model based drug development. This paper proposed an important approach to transform published noncompartment model pharmacokinetics (PK) parameters into compartment model PK parameters. This metaanalysis was performed with a multivariate nonlinear mixed model. A conditional firstorder linearization approach was developed for statistical estimation and inference.
Results
Using MDZ as an example, we showed that this approach successfully transformed 6 noncompartment model PK parameters from 10 publications into 5 compartment model PK parameters. In simulation studies, we showed that this multivariate nonlinear mixed model had little relative bias (<1%) in estimating compartment model PK parameters if all noncompartment PK parameters were reported in every study. If there missing noncompartment PK parameters existed in some published literatures, the relative bias of compartment model PK parameter was still small (<3%). The 95% coverage probabilities of these PK parameter estimates were above 85%.
Conclusions
This noncompartment model PK parameter transformation into compartment model metaanalysis approach possesses valid statistical inference. It can be routinely used for model based drug development.
Background
In recent decades, a new drug requires an average of 15 years and approaching a billion dollars in research and development [1]. Unfortunately, only one in 10 drugs that enter clinical testing receives eventual FDA approval [2, 3]. Scientists have become increasingly mechanistic in their approach to drug development [4]. The recent ability to integrate genetic mutations and altered protein expression to pharmacokinetics (PK) and pharmacodynamic (PD) models allow a deeper understanding of the mechanisms of disease and therapies that are genuinely targeted [5–8]. In 2004, the FDA released a report entitled: “Innovation or Stagnation, Challenge and Opportunity on the Critical Path to New Medical Products” [9]. Among its six general topic areas, three of them emphasized the importance of computational modeling and bioinformatics in biomarker development and streamlining clinical trials [10, 11]. In multiple followup papers, clinical researchers, experimental biologists, computational biologists, and biostatisticians from both academia and industry all supported the FDA leadership in this critical path, and pointed out the challenges and opportunities of the PK/PD model based approach in drug development [12][13–15].
Pharmacokinetics model is the central piece of model based drug development. Almost all of the published PK data were summarized without fitting a compartment model. They are usually called noncompartment model PK parameters. For example, area under the concentration curve (AUC) is calculated from drug plasma concentration data based on trapezoidrule [16]; clearance is calculated from dose and AUC; Cmax and Tmax are calculated from concentrations and their associated time points; terminal halflife is usually calculated from the last two to four sampling timepoints directly; and etc. All these parameters cannot be used directly in a compartment model, and their transformation to compartment model PK parameters is essential.
Methods
NonCompartment Model to OneCompartment Model Transformation
Similarly, if CL_{ iv } is reported, instead of AUC, then CL_{ IV } = V × k_{ e }. These onecompartmentmodel and noncompartment model parameters and transformation were defined and discussed in great detail by [16].
NonCompartment Model to TwoCompartment Model Transformation
A Multivariate Nonlinear Mixed Effect Model (Model Specification)
Model (9) also shows that the observed noncompartment model parameters, , are independent. This is a multivariate nonlinear regression model.
where is a J×1 ( ) observed noncompartment model PK parameter vector; is a J×p indicator matrix, and X_{ k } is a J_{ k }×p matrix indicating the corresponding transformation function; g(.) is a p×1 transformation function vector; is a study level compartmentmodel PK parameter vector; is a diagonal J×J covariance matrix for W, a nd ; is a Kp×p design matrix relating studyspecific parameter β to population parameter µ, and I_{ k } is an identity matrix; and is a Kp×Kp covariance matrix for studyspecific parameter β.
This multivariate nonlinear mixed model (11) is different from the conventional univariate nonlinear mixed model [17] structurally in the additional design matrix X in front of the nonlinear function ( i.e. transformation function g(.)). Model (11) is a metaanalysis approach, in which sample mean noncompartment model PK parameters are formulated. Among the existing nonlinear mixed model metaanalysis literatures, some dealt with the subjectlevel data from multiple studies [18, 19]; the others dealt with sample mean drug concentration data [20, 21]; and none of them discussed the metaanalysis on summarized PK parameters through the noncompartment model.
A Multivariate Nonlinear Mixed Effect Model (Estimation and Inference)
As a conditional first order linearization approach provides the least biased estimate in estimating the PK parameter with comparable efficiency [22, 23]), it is chosen as the estimation approach for this multivariate nonlinear mixed model. This conditional first order linearization approach was firstly introduced by Lindstrom and Bates [24]. We revise their derivation based on our special metaanalysis multivariate nonlinear mixed model (11). This twostep estimation scheme is described as following.
Hence, θ can be estimated through an iterative Fisher algorithm. An alternative derivation of this twostep first order linearization is through a second order Laplace’s approximation [25–27].
Results
Midazolam NonCompartment Model Parameters to Compartment Model Parameters Transformation Data Analysis
Summary of Published NonCompartment Model Midazolam Pharmacokinetics Parameters
NonCompartment PK Parameters  Reported  Missed 

C _{ max }  9  1 
AUC  10  0 
T _{ max }  7  3 
T _{ 1/2,fast }  2  8 
T _{ 1/2,slow }  8  2 
V _{ d }  5  5 
CL _{ iv }  4  6 
A multivariate nonlinear mixed effect model is fitted to these published noncompartment PK parameters to estimate their compartment model PK parameters. The NONMEM code is reported in Appendix I. In this metaanalysis, between study variances are assumed for (V_{ 1, } k_{ a, } k _{ e }). (k_{ 12, } k_{ 21 }) were assumed to be the fixed effects across different studies without random effects, because only two papers published the MDZ distribution information, i.e. T_{ 1/2,fast }. All of the noncompartment model parameters were logtransformed. They were assumed to have the same within study variance in logscale (i.e. same coefficient of variance in the raw scale). All of the compartment model PK parameters were also logtransformed, and their between study standard deviations can be interpreted as coefficient of variance in raw scale.
Midazolam Compartment Model Pharmacokinetics Parameter Estimates
Compartment Model PK Parameters  NonCompartment Model to Compartment Model Transformation  

FixedEffect  
logscale  rawscale  Between Study CV*  
V _{ 1 }  3.5  33.11  10% 
k _{ a }  0.68  1.97  84% 
k _{ 12 }  1.1  0.33   
k _{ 21 }  1.32  0.27   
k _{ e }  0.403  0.67  23% 
WithinStudy CV**  27% 
Simulation Studies
Simulation Schemes
The primary concern of this noncompartment PK parameter transformation to compartment model PK parameter is the bias of PK parameter estimates. Two simulation studies were designed to investigate this problem. In the first simulation, every noncompartment PK parameter was observed for each study. In the second simulation, the same amount of missing data as our MDZ example was assumed to be present.
In each simulation, 1000 simulated data sets were generated. Each data set had 10 studies, and each study reported either all (C_{ max, } AUC, T_{ 1/2,slow }, T_{ 1/2,fast } , V_{ d } , CL_{ iv }) in simulation 1, or a partial amount of (C_{ max, } AUC, T_{ 1/2,slow }, T_{ 1/2,fast } , V_{ d } , CL_{ iv }) in simulation 2. These noncompartment model PK parameters were simulated based on the twocompartment model transformation relationship (5) and (6), their metaanalysis multivariate nonlinear mixed model (9) and (10), and MDZ PK parameter estimates and variances from Table 2.
Simulation Evaluation Criteria
Both fixed effect and variance components were evaluated in the simulation studies. The bias was calculated as the relative bias: abs(trueest)/est; and their 95% coverage probabilities were also reported based on model based 95% confidence interval. Coverage probabilities outside of (92.93, 97.07) were highlighted. The halfwidth of this interval is three times the binomial stand error, which is [(95%)(5%)/1000]^{1/2}=0.6892%. Standard error was also reported based on 1000 simulation results.
Simulation 1 (All Reported and No Missing Data)
Simulation Results with No Missing Data
Estimate  

TRUE(logscale)  Mean  SE  RelativeBias (%)  95% CP  
V _{ 1 }  3.5  3.505  0.110  0.14  0.89  
k _{ a }  0.68  0.680  0.115  0.02  0.87  
FixedEffect  k _{ e }  0.403  0.397  0.112  0.75  0.89  
k _{ 12 }  1.1  1.097  0.088  0.24  0.93  
k _{ 21 }  1.32  1.322  0.053  0.15  0.99  
V _{ 1 }  0.09  0.083  0.045  8.05  0.97  
Between Study Variance  k _{ a }  0.09  0.078  0.053  12.8  0.93  
k _{ e }  0.09  0.085  0.045  5.33  0.96  
Sigma ^{ 2 }  0.01  0.01  0.003  4.43  0.95 
Simulation 2 (With Missing Data)
Simulation Results with Missing Data
Estimate  

TRUE(logscale)  Mean  SD  RelativeBias (%)  95% CP  
V_{1}  3.5  3.494  0.129  0.17  0.92  
k_{a}  0.68  0.672  0.159  1.13  0.87  
FixedEffect  k_{e}  0.403  0.389  0.141  2.84  0.90  
k_{12}  1.1  1.09  0.172  0.59  0.84  
k_{21}  1.32  1.323  0.070  0.19  0.99  
V1  0.09  0.081  0.052  9.82  0.98  
Between Study Variance  ka  0.09  0.082  0.066  9.43  0.95  
ke  0.09  0.087  0.055  3.60  0.97  
Sigma ^{ 2 }  0.01  0.01  0.003  13.9  0.96 
Conclusions
This paper proposed an important approach to transform published noncompartment model pharmacokinetics parameters into compartment model PK parameters. This metaanalysis was performed with a multivariate nonlinear mixed model. A conditional firstorder linearization approach was developed for statistical estimation and inference, and it was implemented in R. Using MDZ as an example, we have shown that this approach transformed 6 noncompartment model PK parameters from 10 publications into 5 compartment model PK parameters, and the conditional first order linearization approach converged to the maximum likelihood. In the followup simulation studies, we have shown that our metaanalysis multivariate nonlinear mixed model had little relative bias (<1%) in estimating compartment model PK parameters if all noncompartment PK parameters were reported in every study. If there existed missing noncompartment PK parameters, the relative bias of compartment model PK parameter was still small (<3%). The 95% coverage probabilities of these PK parameter estimates were usually above 85% or more. Therefore, this approach possesses adequately valid inference.
Although this paper only showed the transformation performance of noncompartment model PK parameters to twocompartment model with oral dose PK parameters, we think it is probably the most complicated case among published drug PK studies. One compartment models and twocompartment model with IV dose have simpler transformation function and less computational expense.
Sometimes, not all of the required noncompartment model PK parameters are available in the literature. Whether it is feasible to transform these data into compartment model is an interesting and important question. In this paper, MDZ was chosen as an example. Because MDZ has been a well studied probe drug, its published noncompartment model PK parameters were expected to be rich. Other rarely studied drugs may not have all these published information, and their compartment model developments from literature need further investigations.
Authors’ information
ZW is currently a Ph.D. Computer Science student in the Indiana University; SK is an assistant professor in the University of Louisville; SKQ is an assistant professor in the Indiana University; JZ is a PhD student in the University of Michigan; and LL is an association professor in the Indiana University.
Notes
List of Abbreviations
 AUC:

area under the concentration curve
 MDZ:

Midazolam
 PK:

Pharmacokinetics.
Declarations
Acknowledgements
Dr. Lang Li is supported by NIH grants, R01 GM74217. Dr. Seongho Kim is partially supported by DOE grants, DEEM0000197, and an Intramural Research Incentive Grant from the University of Louisville.
This article has been published as part of BMC Systems Biology Volume 4 Supplement 1, 2010: Proceedings of the ISIBM International Joint Conferences on Bioinformatics, Systems Biology and Intelligent Computing (IJCBS). The full contents of the supplement are available online at http://www.biomedcentral.com/17520509/4?issue=S1.
Authors’ Affiliations
References
 DiMasi JA, Hansen RW, Grabowski HG: The price of innovation: new estimates of drug development costs. Journal of Health Economics. 2003, 22 (2): 151185. 10.1016/S01676296(02)001261View ArticlePubMed
 Kola I, Landis J: Can the pharmaceutical industry reduce attrition rates?. Nature Reviews Drug Discovery. 2004, 3 (8): 711715. 10.1038/nrd1470View ArticlePubMed
 Woosley RL: Drug development and the FDA’s critical path initiative. Clinical Pharmacology and Therapeutics. 2007, 129133. 81View ArticlePubMed
 Veit M: New strategies for drug development. Berl. Munch. Tierarztl. Wochenschr. 276287. 117
 D'Andrea G e a: A polymorphism in the VKORC1 gene is associated with an interindividual variability in the doseanticoagulant effect of warfarin. Blood. 2005, 645649. 105View ArticlePubMed
 Kirchheiner JB: J. Clinical consequences of cytochrome P450 2C9 polymorphisms. Clin. Pharmacol. Ther. 2005, 116. 77View ArticlePubMed
 Badagnani I e a: Interaction of methotrexate with organicanion transporting polypeptide 1A2 and its genetic variants. J. Pharmacol. Exp. Ther. 521529. 318View Article
 Hung S I e a: Genetic susceptibility to carbamazepineinduced cutaneous adverse drug reactions. Pharmacogenet. Genom. 2006, 297306. 10.1097/01.fpc.0000199500.46842.4a. 16View Article
 Cleary JD, Taylor JW, Chapman SW: Itraconazole in antifungal therapy. Ann Pharmacother. 1992, 26 (4): 5029.PubMed
 Yamazaki M, Nishigaki R, Suzuki H, Sugiyama Y: [Kinetic analysis of hepatobiliary transport of drugs, importance of carriermediated transport]. Yakugaku Zasshi. 1995, 115 (12): 95377.PubMed
 Okuda H, Nishiyama T, Ogura K, Nagayama S, Ikeda K, Yamaguchi S, Nakamura Y, Kawaguchi Y, Watabe T: Lethal drug interactions of sorivudine, a new antiviral drug, with oral 5fluorouracil prodrugs. Drug Metab Dispos. 1997, 25 (5): 2703.
 Lalonde RL, Hutmacher MM, Ewy W, Nichols DJ, Milligan PA, Corrigan BW, Lockwood PA, Marshall SA, Benincosa LJ, Tensfeldt TG, Parivar K, Amantea M, Glue P, Koide H, Miller R: Modelbased drug development. Clin Pharmacol Ther. 2007, 82 (1): 2132. 10.1038/sj.clpt.6100235View ArticlePubMed
 Chang MKS, Bull J, Chiu YY, Wang W, Wakeford C, McCarthy K: Innovative approaches in drug development. J Biopharm Stat. 2007, 17 (5): 77589. 10.1080/10543400701513926View ArticlePubMed
 RT ON: FDA's critical path initiative: a perspective on contributions of biostatistics. Biom J. 2006, 48 (4): 55964. 10.1002/bimj.200510237View Article
 Chien JY, Heathman MA, de Alwis DP, Sinha V: Pharmacokinetics/Pharmacodynamics and the stages of drug development: role of modeling and simulation. The AAPS Journal. 2005, 7 (3): E544E559. 10.1208/aapsj070355PubMed CentralView ArticlePubMed
 Rowland M, Tozer TN: Clinical Pharmacokinetics Concept and Applications. 1995, Lippincott Williams & Wilkins, Third
 Davidian M, Giltinan DM: Nonlinear models for repeated measurment data. 1995, Chapman and Hall
 Wakefield JC, Rahman N: The combination of population pharmacokinetic studies. Biometrics. 2000, 56 (1): 26370. 10.1111/j.0006341X.2000.00263.xView ArticlePubMed
 Lopes HF, Muller P, Rosner GL: Bayesian metaanalysis for longitudinal data models using multivariate mixture priors. Biometrics. 2003, 59 (1): 6675. 10.1111/15410420.00008View ArticlePubMed
 Li L, Yu M, Chin R, Lucksiri A, Flockhart D, Hall S: DrugDrug Interaction Prediction: A Bayesian MetaAnalysis Approach. Statistis in Medicine. 2007, 26 (20): 37003721. 10.1002/sim.2837.View Article
 Yu M, et al.: A Bayesian metaanalysis on published sample mean and variance pharmacokinetic data with application to drugdrug interaction prediction. J Biopharm Stat. 2008, 18 (6): 106383. 10.1080/10543400802369004PubMed CentralView ArticlePubMed
 Wolfinger RD: Two Taylorseries approximations methods for nonlinear mixed models. Journal of Computational Statistics and Data Analysis. 1997, 25: 465490. 10.1016/S01679473(97)000121.View Article
 Wolfinger R, O'Connell M: Generalized linear mixed models: a pseudolikelihood approach. Journal of statistical Computation and Simulation. 1993, 48: 233243. 10.1080/00949659308811554.View Article
 Lindstrom ML, Bates DM: Nonlinear mixed effects models for repeated measures data. Biometrics. 1990, 46 (3): 67387. 10.2307/2532087View ArticlePubMed
 Wolfinger RD: Laplace's Approximation for Nonlinear Mixed Models. Biometrika. 1993, 80: 791795. 10.1093/biomet/80.4.791.View Article
 Westfall P, Hochberg Y, Wolfinger RD, Rorn D, Tobias RD: Multiple comparisons and Multiple Tests. 1999
 Vonesh EF: A Note on Laplace's Approximation in Nonlinear Mixed Effects Models. Biometrika. 1996, 83: 447452. 10.1093/biomet/83.2.447.View Article
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.