620 likes | 1.52k Views
Pharmacometrics Introduction . Yaning Wang, Ph.D. Team Leader, Pharmacometrics Office of Clinical Pharmacology Center for Drug Evaluation and Research Food and Drug Administration yaning.wang@fda.hhs.gov. Disclaimer: My remarks today do not necessarily reflect the official views of FDA.
E N D
Pharmacometrics Introduction Yaning Wang, Ph.D. Team Leader, Pharmacometrics Office of Clinical Pharmacology Center for Drug Evaluation and Research Food and Drug Administration yaning.wang@fda.hhs.gov Disclaimer: My remarks today do not necessarily reflect the official views of FDA NONMEM Estimation Methods
Outline • Issues and opportunities in drug development • Model-based drug development • Bayesian statistics and its relationship with NONMEM • NONMEM estimation methods NONMEM Estimation Methods
Issues in Drug Development • Low efficiency • NME IND = NDA <20% of time • Reported >50% failure rate in Phase 3 (Carl Peck, CDDS) • Decreased NME NDAs despite increased INDs • Cost per NME approved estimated at $1.7B NONMEM Estimation Methods
Low Success Rate Across Different Diseases Kola I, Landis J.Can the pharmaceutical industry reduce attrition rates? Nat.Rev.Drug.Disc. Aug 2004. NONMEM Estimation Methods
High Failure Rate Even in Late Development Kola I, Landis J.Can the pharmaceutical industry reduce attrition rates? Nat.Rev.Drug.Disc. Aug 2004. NONMEM Estimation Methods
10-Year Trends of Major Submissions to FDA NONMEM Estimation Methods
Investment Escalation Source: Windhover’s in Vivo: The Business & Medicine Report, Bain drug economics model, 2003 NONMEM Estimation Methods
Various Initiatives • National Institutes of Health (NIH) • Roadmap initiative • National Cancer Institute (NCI) • Specialized Programs of Research Excellence (SPOREs) • European Organization for the Treatment of Cancer (EORTC) • make translational research a part of all cancer clinical trials • National Translational Cancer Research Network • facilitate and enhance translational research in the United Kingdom • FDA • Critical Path Initiative NONMEM Estimation Methods
Application of quantitative disease models to drug development Pharmacogenomics in drug development New imaging technologies may contribute biomarkers in drug development 2004 FDA Critical Path Initiative Application of New Scientific Knowledge to Drug Development: http://www.fda.gov/oc/initiatives/criticalpath/ NONMEM Estimation Methods
2006 FDA Critical Path Update NONMEM Estimation Methods
Model-Based Drug Development “The concept of model-based drug development, in which pharmaco-statistical models of drug efficacy and safety are developed from preclinical and available clinical data, offers an important approach to improving drug development knowledge management and development decision-making”Adapted from Lewis B. Sheiner, “Learning vs Confirming in Clinical Drug Development”, Clin. Pharmacol. Ther., 1997, 61:275-291. NONMEM Estimation Methods
Terminology • Model-based drug development • Pharmacokinetics/Pharmacodynamics (PK/PD) • Modeling and simulation • Exposure-response • Quantitative clinical pharmacology • Quantitative disease and drug models • Pharmacometrics-Pharmacometricians • Areas involved: clinical pharmacology, statistics, pathophysiology,biology, bioengineering NONMEM Estimation Methods
DRUG MODEL PLACEBO/DISEASE MODEL Morbidity#1 Relative Risk Morbidity#2 Surrogate Toxicity Mortality Surrogate Exposure Disease Drug Trial Models Surrogate Exposure High dose Low dose TIME TIME PATIENT/CLINICAL TRIAL MODEL % Drop-out % Compliance QD TID TIME Surrogate Body Weight Patient Population NONMEM Estimation Methods
Applications in FDA Review • Exposure-response analysis of efficacy and safety data in NDA review for the choice of dosing regimen (s) • Dose adjustment in special populations (hepatic, renal, gender, age and drug interactions) based on intersubject variability and risk assessment • Routine use of population PK and PD data analysis to understand variability and to provide evidence for label claims • Issuing guidance to assist the industry in using these tools • Case studies • Leveraging Prior Quantitative Knowledge to Guide Drug Development Decisions and Regulatory Science Recommendations: Impact of FDA Pharmacometrics During 2004-2006, Journal of Clinical Pharmacology (2008) 48: 146-157 • Impact of Pharmacometric Reviews on New Drug Approval and Labeling Decisions - A Survey of 31 New Drug Applications Submitted Between 2005-2006, Clinical Pharmacology and Therapeutics (2007) 81: 213-221 • Impact of Pharmacometrics on Drug Approval and Labeling Decisions - A Survey of 42 New Drug Applications, The AAPS Journal, 7(3): E503-E512, 2005 NONMEM Estimation Methods
Tools for Modeling and Simulation • NONMEM (UCSF, Globomax) • SAS (SAS Institute Inc) • Splus (Insightful Corporation) or R (Free) • WinBUGS (MRC Biostatistics, Free) • ADAPT II (USC, Free) • WinNonLin/WinNonMix (Pharsight) • Trial Simulator (Pharsight) NONMEM Estimation Methods
0.6 0.5 0.4 0.3 0.2 0.1 0.0 Basic Statistics 1.0 • Random Variable (Y) • Discrete (gender, pain score) • Continuous (body weight, clearance) • Distribution (histogram) • Binomial • Normal • Lognormal • Probability Function (P(Y=y)) • Probability Mass Function (PMF) • Probability Density Function (PDF) 0.8 0.6 0.4 0.2 0.0 0 1 0.02 0.015 0.010 0.005 0.0 20 40 60 80 100 120 140 1 2 3 4 5 6 7 NONMEM Estimation Methods
Normal Distribution 0.04 0.03 P(y) 0.02 Y~N(100, 102) 0.01 0.0 • 100 120 • Y NONMEM Estimation Methods
Joint, Marginal and Conditional Probabilities P(Y2) P(Y1) Joint probability: P(Y1=0,Y2=0)=0.72 Marginal probability: P(Y1=0)= 0.90 Conditional probability: P(Y2=0|Y1=0)= 0.72/0.90=0.80 NONMEM Estimation Methods
Marginal and Conditional Probabilities P(Y2) P(Y1) Conditional probability: Marginal probability: For continuous variables: Marginal probability: weighted average of conditional probability NONMEM Estimation Methods
Bayes’ Theorem NONMEM Estimation Methods
Parameter • Unknown but fixed (Frequentist) • Unknown and random (Bayesian) • Unknown hyperparameters • Hyperprior (Full Bayesian) • Estimation (Empirical Bayesian) (Prior distribution of ) Hyperparameters NONMEM Estimation Methods
Prior, Likelihood and Posterior Likelihood Posterior Distribution Prior Distribution Marginal Distribution NONMEM Estimation Methods
Prior, Likelihood and Posterior NONMEM Estimation Methods
Posterior Likelihood Prior 2 n Shrinkage NONMEM Estimation Methods
A Simple Example • Unknown parameter (): long-term systolic blood pressure (SBP) of one particular 60-year-old female • 4 measurements with a mean and a standard deviation =5 • Survey of the same population (60-year-old female): mean SBP =120 and standard deviation =10 NONMEM Estimation Methods
Estimation of • Frequentist • Point estimate: • Interval estimate (95%CI) • Bayesian • Posterior distribution P(|Y) • Posterior mean • 95% credible interval NONMEM Estimation Methods
Individual parameter (posterior) Individual data (likelihood) Population Distribution (prior) Prior, Likelihood and Posterior Long-term systolic blood pressure (SBP) of 60-year-old woman NONMEM Estimation Methods
Population PK Variability Sheiner, 1992 NONMEM Estimation Methods
NONMEM and Bayesian i: Posthoc estimate (Empirical Bayesian Estimate) because , , 2 are all estimated based on MLE NONMEM Estimation Methods
A Simple Example yij: the jth observation for the ith individual Assume: only k is the unknown parameter to be estimated Goal: search for the estimate, , that can minimize the following objective function where Li(k) is the marginal likelihood of k for ith individual NONMEM Estimation Methods
What Is Marginal Likelihood ? P(i) P(Yi|k) P(Yi|k,i) P(Yi,i|k) This is just the area under the curve (AUC) of h(,k) versus at a fixed k • Each fixed k leads to different AUC • Assuming only 1 subject, is associated with maximum AUC. is called the maximum likelihood estimator (MLE) of k • Even though AUC is a function of k, it cannot be expressed as a closed-form equation of k (difficulty of nonlinear mixed effect modeling) k=1 h () NONMEM Estimation Methods
Grid Search to Find • Imagine we can try every possible k and evaluate AUC with numerical integration (similar to trapezoidal rule with very dense data) k=0.3 k=0.7 k=1.1 h () y AUC (k) k=0.3 k=0.7 k=1.1 Time (hr) k NONMEM Estimation Methods
Estimation Methods • NONMEM • Laplacian • First order conditional estimate (FOCE) • First order (FO) • SAS • Adaptive Gaussian Quadrature • Importance sampling • FO • Splus • Lindstrom and Bates Algorithm NONMEM Estimation Methods
What is LAPLACIAN doing? • Approximate h() with another function LAP() so that AUC is a close-form function of k k=0.3 k=0.7 k=1.1 LAP () or h () NONMEM Estimation Methods
What is FOCE doing? • After Laplacian approximation, the second derivative, , is further approximated by , a function of the first derivative k=0.3 k=0.7 k=1.1 FOCE () or h () NONMEM Estimation Methods
What is FO doing? • Approximate h() with another function FO() and also approximate the second derivative, , with a function of the first derivative k=0.3 k=0.7 k=1.1 FO () or h () NONMEM Estimation Methods
Which one is the best? True LAP FOCE FO Marginal Likelihood LAP is the best even though occasionally FOCE, even FO, is closer to the true marginal likelihood NONMEM Estimation Methods
Proportional Residual Error Model yij: the jth observation for the ith individual Assume: only k is the unknown parameter to be estimated NONMEM Estimation Methods
LAPI* FOCEI** LAP FOCE FO k=0.3 k=0.7 k=1.1 *: LAPI, Laplacian with interaction is available in NONMEM VI **: FOCEI, first-order conditional estimate with interaction NONMEM Estimation Methods
Which one is the best? True LAPI FOCEI LAP FOCE FO Marginal Likelihood LAPI is the best and LAP is worse than FOCEI for proportional error model k NONMEM Estimation Methods
A Two-Dimensional Example yij: the jth observation for the ith individual Assume: k and V are the unknown parameters to be estimated Goal: search for the estimates, and , that can minimize the following objective function where Li(k, V) is the marginal likelihood of k and V for ith individual NONMEM Estimation Methods
What Is Marginal Likelihood ? This is just the volume under the surface (VUS) of h(1, 2, , k, V) versus 1 and 2 at a fixed (k,V) k=1 and V=1 • Each fixed pair (k,V) leads to different VUS • If and are associated with maximum VUS. They are the maximum likelihood estimators (MLEs) of k and V • Even though VUS is a function of k and V, it cannot be expressed as a closed-form equation of k and V h (1, 2) 2 1 NONMEM Estimation Methods
LAPLACIAN, FOCE, FO Volumes vs True Volume LAPLACIAN FOCE FO h (1, 2) 1 1 1 2 2 2 h (1, 2) 2 1 2 2 1 1 NONMEM Estimation Methods
Summary for Estimation Methods • The approximation of true marginal likelihood is analogous to approximating an irregular shape with a symmetric shape whose AUC or VUS can be easily calculated • The difference among LAPLACIAN, FOCE and FO is mainly the thinness and height of the symmetric shapes • In general, the shape generated by LAPLACIAN method is the closest to the true shape, but the shape from FOCE is almost equally close • For models with - interaction, taking the interaction into account seems more important than avoiding first-order approximation (FOCEI vs LAP) Reference: Derivation of Various NONMEM Estimation Methods, Yaning Wang, Journal of Pharmacokinetics and Pharmacodynamics, 34: 575-93, 2007 NONMEM Estimation Methods
Full Bayesian for PopPK NONMEM Estimation Methods
Acyclic Structure H R a b MVN W G MVN Dose Model Time N Data NONMEM Estimation Methods