Uncertainty Computation at Finite Distance in Nonlinear Mixed Effects Models-a New Method Based on Metropolis-Hastings Algorithm.

The standard errors (SE) of the maximum likelihood estimates (MLE) of the population parameter vector in nonlinear mixed effect models (NLMEM) are usually estimated using the inverse of the Fisher information matrix (FIM). However, at a finite distance, i.e. far from the asymptotic, the FIM can underestimate the SE of NLMEM parameters. Alternatively, the standard deviation of the posterior distribution, obtained in Stan via the Hamiltonian Monte Carlo algorithm, has been shown to be a proxy for the SE, since, under some regularity conditions on the prior, the limiting distributions of the MLE and of the maximum a posterior estimator in a Bayesian framework are equivalent. In this work, we develop a similar method using the Metropolis-Hastings (MH) algorithm in parallel to the stochastic approximation expectation maximisation (SAEM) algorithm, implemented in the saemix R package. We assess this method on different simulation scenarios and data from a real case study, comparing it to other SE computation methods. The simulation study shows that our method improves the results obtained with frequentist methods at finite distance. However, it performed poorly in a scenario with the high variability and correlations observed in the real case study, stressing the need for calibration.

Pharmacometrics enhanced Bayesian borrowing approach to improve clinical trial efficiency: Case of empagliflozin in type 2 diabetes.

We report use of a pharmacometrics enhanced Bayesian borrowing (PEBB) approach to leverage historical clinical trial data on a drug product to build models, project the outcome of future clinical trials, and borrow information from these projections to improve the efficiency of future target trials. This design takes a two-stage approach. First, a design phase is performed before target trial data are available to determine the operating characteristics and an appropriate tuning parameter that will be used in the subsequent analysis phase of a chosen target trial. Second, once the target trial data are available, the analysis phase is performed with the determined tuning parameter. This step is where borrowing is applied from these projections to inform the results for the target trial. To illustrate how a PEBB could improve the efficiency of clinical trials, we apply our design to trials with empagliflozin for treating patients with type 2 diabetes. We performed a retrospective evaluation applying the method to a phase III target trial and a hypothetical smaller trial. The type I error could be kept below 10% while increasing the trial power and effective sample size. Our findings suggest that a PEBB has the potential to increase the power of clinical trials, while controlling for type I error, by leveraging the information from previous trials through population pharmacokinetic/pharmacodynamic modeling and simulation.

Robust and Adaptive Two-stage Designs in Nonlinear Mixed Effect Models.

To get informative studies for nonlinear mixed effect models (NLMEM), design optimization can be performed based on Fisher Information Matrix (FIM) using the D-criterion. Its computation requires knowledge about models and parameters, which are often prior guesses. Thus, adaptive designs composed of several stages may be used. Robust approach can also be used to account for various candidate models. In the estimation step of a given stage, model selection (MS) or model averaging (MA) can be performed. In this work we propose a new two-stage adaptive design strategy, based on the robust expected FIM and MA over several candidate models. The methodology is applied to a clinical trial simulation in ophthalmology to optimize doses and time measurements. A set of dose-response candidate models is defined, and one-stage designs are compared to two-stage 50/50 designs (i.e., each stage performed with half of the available subjects), using either local optimal design or robust design, and performing analysis with one model, MS or MA. Performing a two-stage design with MS at the interim analysis can correct the choice of a wrong model for designing the first stage. Overall, starting from a robust design (1- or 2-stage) is valuable and leads to reasonable bias and precision. The proposed robust adaptive design strategy is a new tool to design longitudinal studies that could be used in different therapeutic areas.

Pharmacokinetic Model of Tenofovir and Emtricitabine and Their Intracellular Metabolites in Patients in the ANRS 134-COPHAR 3 Trial Using Dose Records.

Tenofovir (TFV) and emtricitabine (FTC) are part of the recommended highly active antiretroviral therapy (ART). Both molecules show a large interindividual pharmacokinetic (PK) variability. Here, we modeled the concentrations of plasma TFV and FTC and their intracellular metabolites (TFV diphosphate [TFV-DP] and FTC triphosphate [FTC-TP]) collected after 4 and 24 weeks of treatment in 34 patients from the ANRS 134-COPHAR 3 trial. These patients received daily (QD) atazanavir (300 mg), ritonavir (100 mg), and a fixed-dose combination of coformulated TFV disoproxil fumarate (300 mg) and FTC (200 mg). Dosing history was collected using a medication event monitoring system. A three-compartment model with absorption delay (Tlag) was selected to describe the PK of, respectively, TFV/TFV-DP and FTC/FTC-TP. TFV and FTC apparent clearances, 114 L/h (relative standard error [RSE] = 8%) and 18.1 L/h (RSE = 5%), respectively, were found to decrease with age. However, no significant association was found with the polymorphisms ABCC2 rs717620, ABCC4 rs1751034, and ABCB1 rs1045642. The model allows prediction of TFV-DP and FTC-TP concentrations at steady state with alternative regimens.