Search

Purpose Tepotinib is a highly selective MET inhibitor approved for treatment of non-small cell lung cancer (NSCLC) harboring METex14 skipping alterations. Analyses presented herein evaluated the relationship between tepotinib exposure, and efficacy and safety outcomes. Methods Exposure–efficacy analyses included data from an ongoing phase 2 study (VISION) investigating 500 mg/day tepotinib in NSCLC harboring METex14 skipping alterations. Efficacy endpoints included objective response, duration of response, and progression-free survival. Exposure–safety analyses included data from VISION, plus four completed studies in advanced solid tumors/hepatocellular carcinoma (30–1400 mg). Safety endpoints included edema, serum albumin, creatinine, amylase, lipase, alanine aminotransferase, aspartate aminotransferase, and QT interval corrected using Fridericia’s method (QTcF). Results Tepotinib exhibited flat exposure–efficacy relationships for all endpoints within the exposure range observed with 500 mg/day. Tepotinib also exhibited flat exposure–safety relationships for all endpoints within the exposure range observed with 30–1400 mg doses. Edema is the most frequently reported adverse event and the most frequent cause of tepotinib dose reductions and interruptions; however, the effect plateaued at low exposures. Concentration-QTc analyses using data from 30 to 1400 mg tepotinib resulted in the upper bounds of the 90% confidence interval being less than 10 ms for the mean exposures at the therapeutic (500 mg) and supratherapeutic (1000 mg) doses. Conclusions These analyses provide important quantitative pharmacologic support for benefit/risk assessment of the 500 mg/day dosage of tepotinib as being appropriate for the treatment of NSCLC harboring METex14 skipping alterations. Registration Numbers NCT01014936, NCT01832506, NCT01988493, NCT02115373, NCT02864992.

AO_SCPLOWBSTRACTC_SCPLOWPharmacometric modeling establishes causal quantitative relationship between administered dose, tissue exposures, desired and undesired effects and patients risk factors. These models are employed to de-risk drug development and guide precision medicine decisions. Recent technological advances rendered collecting real-time and detailed data easy. However, the pharmacometric tools have not been designed to handle heterogeneous, big data and complex models. The estimation methods are outdated to solve modern healthcare challenges. We set out to design a platform that facilitates domain specific modeling and its integration with modern analytics to foster innovation and readiness to data deluge in healthcare. New specialized estimation methodologies have been developed that allow dramatic performance advances in areas that have not seen major improvements in decades. New ODE solver algorithms, such as coefficient-optimized higher order integrators and new automatic stiffness detecting algorithms which are robust to frequent discontinuities, give rise to up to 4x performance improvements across a wide range of stiff and non-stiff systems seen in pharmacometric applications. These methods combine with JIT compiler techniques and further specialize the solution process on the individual systems, allowing statically-sized optimizations and discrete sensitivity analysis via forward-mode automatic differentiation, to further enhance the accuracy and performance of the solving and parameter estimation process. We demonstrate that when all of these techniques are combined with a validated clinical trial dosing mechanism and non-compartmental analysis (NCA) suite, real applications like NLME parameter estimation see run times halved while retaining the same accuracy. Meanwhile in areas with less prior optimization of software, like optimal experimental design, we see orders of magnitude performance enhancements. Together we show a fast and modern domain specific modeling framework which lays a platform for innovation via upcoming integrations with modern analytics.

Editor’s Perspective What We Already Know about This Topic What This Article Tells Us That Is New Background Population-based, pharmacokinetic modeling can be used to describe variability in fluid distribution and dilution between individuals and across populations. The authors hypothesized that dilution produced by crystalloid infusion after hemorrhage would be larger in anesthetized than in awake subjects and that population kinetic modeling would identify differences in covariates. Methods Twelve healthy volunteers, seven females and five males, mean age 28 ± 4.3 yr, underwent a randomized crossover study. Each subject participated in two separate sessions, separated by four weeks, in which they were assigned to an awake or an anesthetized arm. After a baseline period, hemorrhage (7 ml/kg during 20 min) was induced, immediately followed by a 25 ml/kg infusion during 20 min of 0.9% saline. Hemoglobin concentrations, sampled every 5 min for 60 min then every 10 min for an additional 120 min, were used for population kinetic modeling. Covariates, including body weight, sex, and study arm (awake or anesthetized), were tested in the model building. The change in dilution was studied by analyzing area under the curve and maximum plasma dilution. Results Anesthetized subjects had larger plasma dilution than awake subjects. The analysis showed that females increased area under the curve and maximum plasma dilution by 17% (with 95% CI, 1.08 to 1.38 and 1.07 to 1.39) compared with men, and study arm (anesthetized increased area under the curve by 99% [0.88 to 2.45] and maximum plasma dilution by 35% [0.71 to 1.63]) impacted the plasma dilution whereas a 10-kg increase of body weight resulted in a small change (less than1% [0.93 to 1.20]) in area under the curve and maximum plasma dilution. Mean arterial pressure was lower in subjects while anesthetized (P < 0.001). Conclusions In awake and anesthetized subjects subjected to controlled hemorrhage, plasma dilution increased with anesthesia, female sex, and lower body weight. Neither study arm nor body weight impact on area under the curve or maximum plasma dilution were statistically significant and therefore no effect can be established.

SummaryTwo full model approaches was compared with respect to their ability to handle missing covariate information. The reference data analysis approach was the full model method in which the covariate effects are estimated conventionally using fixed effects, and missing covariate data is imputed with the median of the non-missing covariate information. This approach was compared to a novel full model method which treats the covariate data as observed data and estimates the covariates as random effects. A consequence of this way of handling the covariates is that no covariate imputation is required and that any missingness in the covariates is handled implicitly. The comparison between the two analysis methods was based on simulated data from a model of height for age z-scores as a function of age. Data was simulated with increasing degrees of randomly missing covariate information (0-90%) and analyzed using each of the two analysis approaches. Not surprisingly, the precision in the parameter estimates from both methods decreased with increasing degrees of missing covariate information. However, while the bias in the parameter estimates increased in a similar fashion for the reference method, the full random effects approach provided unbiased estimates for all degrees of covariate missingness.

Edoxaban exposure-response relationships from the phase III study evaluating edoxaban for prevention and treatment of venous thromboembolism (VTE) in patients with acute deep vein thrombosis (DVT) and/or pulmonary embolism (PE) were assessed by parametric time-to-event analysis. Statistical significant exposure-response relationships were recurrent VTE with hazard ratio (HR) based on average edoxaban concentration at steady state (C-av) (HRCav) 50.98 (i.e., change in the HR with every 1 ng/mL increase of C-av); the composite of recurrent DVT and nonfatal PE with HRC(av)50.99; and the composite of recurrent DVT, nonfatal PE, and all-cause mortality HRC(av)50.98, and all death using maximal edoxaban concentration (C-max) with HR (C-max) 50.99. No statistical significant exposure-response relationships were found for clinically relevant bleeding or major adverse cardiovascular event. Results support the recommendation of once-daily edoxaban 60 mg, and a reduced 30 mg dose in patients with moderate renal impairment, body weight

Aims This study characterized the population pharmacokinetics of edoxaban in patients with symptomatic deep-vein thrombosis and/or pulmonary embolism in the Hokusai-VTE phase 3 study. The impact of the protocol-specified 50% dose reductions applied to patients with body weight ≤ 60 kg, creatinine clearance (CLcr) of 30 to 50 ml min–1 or concomitant P-glycoprotein inhibitor on edoxaban exposure was assessed using simulations. Methods The sparse data from Hokusai-VTE, 9531 concentrations collected from 3707 patients, were pooled with data from 13 phase 1 studies. In the analysis, the covariate relationships used for dose reductions were estimated and differences between healthy subjects and patients as well as additional covariate effects of age, race and gender were explored based on statistical and clinical significance. Results A linear two-compartment model with first order absorption preceded by a lag time best described the data. Allometrically scaled body weight was included on disposition parameters. Apparent clearance was parameterized as non-renal and renal. The latter increased non-linearly with increasing CLcr. Compared with healthy volunteers, inter-compartmental clearance and the CLcr covariate effect were different in patients (+64.6% and +274%). Asian patients had a 22.6% increased apparent central volume of distribution. The effect of co-administration of P-glycoprotein inhibitors seen in phase 1 could not be confirmed in the phase 3 data. Model-based simulations revealed lower exposure in dose-reduced compared with non-dose-reduced patients. Conclusions The adopted dose-reduction strategy resulted in reduced exposure compared with non-dose-reduced, thereby overcompensating for covariate effects. The clinical impact of these differences on safety and efficacy remains to be evaluated.

Population Pharmacokinetic (PK)-Pharmacodynamic (PD) (PKPD) models are increasingly used in drug development and in academic research. Hence designing efficient studies is an important task. Follow ...

With the increasing popularity of optimal design in drug development it is important to understand how the approximations and implementations of the Fisher information matrix (FIM) affect the resulting optimal designs. The aim of this work was to investigate the impact on design performance when using two common approximations to the population model and the full or block-diagonal FIM implementations for optimization of sampling points. Sampling schedules for two example experiments based on population models were optimized using the FO and FOCE approximations and the full and block-diagonal FIM implementations. The number of support points was compared between the designs for each example experiment. The performance of these designs based on simulation/estimations was investigated by computing bias of the parameters as well as through the use of an empirical D-criterion confidence interval. Simulations were performed when the design was computed with the true parameter values as well as with misspecified parameter values. The FOCE approximation and the Full FIM implementation yielded designs with more support points and less clustering of sample points than designs optimized with the FO approximation and the block-diagonal implementation. The D-criterion confidence intervals showed no performance differences between the full and block diagonal FIM optimal designs when assuming true parameter values. However, the FO approximated block-reduced FIM designs had higher bias than the other designs. When assuming parameter misspecification in the design evaluation, the FO Full FIM optimal design was superior to the FO block-diagonal FIM design in both of the examples.

Positron emission tomography (PET) is an imaging technique that is used to investigate ligand-receptor binding in the living brain and to determine the time course of plasma concentration/receptor occupancy (RO). The purpose of this work was to demonstrate the added value of an adaptive-optimal design for PET scan timings and dose selection over traditional study designs involving fixed or educated selections of timings and doses. A k(on)-k(off) model relating plasma concentration to PET data was applied to generate the simulated data. Optimization was performed on scanning timings and doses using the D-optimality criterion. Optimal designs as applied to scanning timings provided unbiased estimates and improved the accuracy of results relative to those of fixed and educated designs. Optimization of both timings and dose provided improvements in accuracy and precision when the initial dose selection was noninformative regarding the time course of RO. These results indicate that adaptive-optimal designs can provide an efficient experimental design for RO studies using PET, by minimizing the number of subjects required and maximizing information related to the plasma concentration-RO relationship.

Inter occasion variability (IOV) is of importance to consider in the development of a design where individual pharmacokinetic or pharmacodynamic parameters are of interest. IOV may adversely affect the precision of maximum a posteriori (MAP) estimated individual parameters, yet the influence of inclusion of IOV in optimal design for estimation of individual parameters has not been investigated. In this work two methods of including IOV in the maximum a posteriori Fisher information matrix (FIMMAP) are evaluated: (i) MAP(occ)-the IOV is included as a fixed effect deviation per occasion and individual, and (ii) POPocc-the IOV is included as an occasion random effect. Sparse sampling schedules were designed for two test models and compared to a scenario where IOV is ignored, either by omitting known IOV (Omit) or by mimicking a situation where unknown IOV has inflated the IIV (Inflate). Accounting for IOV in the FIMMAP markedly affected the designs compared to ignoring IOV and, as evaluated by stochastic simulation and estimation, resulted in superior precision in the individual parameters. In addition MAP(occ) and POPocc accurately predicted precision and shrinkage. For the investigated designs, the MAP(occ) method was on average slightly superior to POPocc and was less computationally intensive.

International audience; : Population Pharmacokinetic (PK)-Pharmacodynamic (PD) (PKPD) models are increasingly used in drug development and in academic research. Hence designing efficient studies is an important task. Following the first theoretical work on optimal design for nonlinear mixed effect models, this research theme has grown rapidly. There are now several different software tools that implement an evaluation of the Fisher information matrix for population PKPD. We compared and evaluated five software tools: PFIM, PkStaMP, PopDes, PopED, and POPT. The comparisons were performed using two models: i) a simple one compartment warfarin PK model; ii) a more complex PKPD model for Pegylated-interferon (peg-interferon) with both concentration and response of viral load of hepatitis C virus (HCV) data. The results of the software were compared in terms of the standard error values of the parameters (SE) predicted from the software and the empirical SE values obtained via replicated clinical trial simulation and estimation. For the warfarin PK model and the peg-interferon PKPD model all software gave similar results. Of interest it was seen, for all software, that the simpler approximation to the Fisher information matrix, using the block diagonal matrix, provided predicted SE values that were closer to the empirical SE values than when the more complicated approximation was used (the full matrix). For most PKPD models, using any of the available software tools will provide meaningful results, avoiding cumbersome simulation and allowing design optimisation.

D-optimal designs for discrete-type responses have been derived using generalized linear mixed models, simulation based methods and analytical approximations for computing the fisher information matrix (FIM) of non-linear mixed effect models with homogeneous probabilities over time. In this work, D-optimal designs using an analytical approximation of the FIM for a dichotomous, non-homogeneous, Markov-chain phase advanced sleep non-linear mixed effect model was investigated. The non-linear mixed effect model consisted of transition probabilities of dichotomous sleep data estimated as logistic functions using piecewise linear functions. Theoretical linear and nonlinear dose effects were added to the transition probabilities to modify the probability of being in either sleep stage. D-optimal designs were computed by determining an analytical approximation the FIM for each Markov component (one where the previous state was awake and another where the previous state was asleep). Each Markov component FIM was weighted either equally or by the average probability of response being awake or asleep over the night and summed to derive the total FIM (FIM(total)). The reference designs were placebo, 0.1, 1-, 6-, 10- and 20-mg dosing for a 2- to 6-way crossover study in six dosing groups. Optimized design variables were dose and number of subjects in each dose group. The designs were validated using stochastic simulation/re-estimation (SSE). Contrary to expectations, the predicted parameter uncertainty obtained via FIM(total) was larger than the uncertainty in parameter estimates computed by SSE. Nevertheless, the D-optimal designs decreased the uncertainty of parameter estimates relative to the reference designs. Additionally, the improvement for the D-optimal designs were more pronounced using SSE than predicted via FIM(total). Through the use of an approximate analytic solution and weighting schemes, the FIM(total) for a non-homogeneous, dichotomous Markov-chain phase advanced sleep model was computed and provided more efficient trial designs and increased nonlinear mixed-effects modeling parameter precision.

Efficient power calculation methods have previously been suggested for Wald test-based inference in mixed-effects models but the only available alternative for Likelihood ratio test-based hypothesis testing has been to perform computer-intensive multiple simulations and re-estimations. The proposed Monte Carlo Mapped Power (MCMP) method is based on the use of the difference in individual objective function values (ΔiOFV) derived from a large dataset simulated from a full model and subsequently re-estimated with the full and reduced models. The ΔiOFV is sampled and summed (∑ΔiOFVs) for each study at each sample size of interest to study, and the percentage of ∑ΔiOFVs greater than the significance criterion is taken as the power. The power versus sample size relationship established via the MCMP method was compared to traditional assessment of model-based power for six different pharmacokinetic and pharmacodynamic models and designs. In each case, 1,000 simulated datasets were analysed with the full and reduced models. There was concordance in power between the traditional and MCMP methods such that for 90% power, the difference in required sample size was in most investigated cases less than 10%. The MCMP method was able to provide relevant power information for a representative pharmacometric model at less than 1% of the run-time of an SSE. The suggested MCMP method provides a fast and accurate prediction of the power and sample size relationship.

We propose a new method for optimal experimental design of population pharmacometric experiments based on global search methods using interval analysis; all variables and parameters are represented as intervals rather than real numbers. The evaluation of a specific design is based on multiple simulations and parameter estimations. The method requires no prior point estimates for the parameters, since the parameters can incorporate any level of uncertainty. In this respect, it is similar to robust optimal design. Representing sampling times and covariates like doses by intervals gives a direct way of optimizing with rigorous sampling and dose intervals that can be useful in clinical practice. Furthermore, the method works on underdetermined problems for which traditional methods typically fail.

Intravenous glucose tolerance test (IVGTT) provocations are informative, but complex and laborious, for studying the glucose-insulin system. The objective of this study was to evaluate, through optimal design methodology, the possibilities of more informative and/or less laborious study design of the insulin modified IVGTT in type 2 diabetic patients. A previously developed model for glucose and insulin regulation was implemented in the optimal design software PopED 2.0. The following aspects of the study design of the insulin modified IVGTT were evaluated; (1) glucose dose, (2) insulin infusion, (3) combination of (1) and (2), (4) sampling times, (5) exclusion of labeled glucose. Constraints were incorporated to avoid prolonged hyper- and/or hypoglycemia and a reduced design was used to decrease run times. Design efficiency was calculated as a measure of the improvement with an optimal design compared to the basic design. The results showed that the design of the insulin modified IVGTT could be substantially improved by the use of an optimized design compared to the standard design and that it was possible to use a reduced number of samples. Optimization of sample times gave the largest improvement followed by insulin dose. The results further showed that it was possible to reduce the total sample time with only a minor loss in efficiency. Simulations confirmed the predictions from PopED. The predicted uncertainty of parameter estimates (CV) was low in all tested cases, despite the reduction in the number of samples/subject. The best design had a predicted average CV of parameter estimates of 19.5%. We conclude that improvement can be made to the design of the insulin modified IVGTT and that the most important design factor was the placement of sample times followed by the use of an optimal insulin dose. This paper illustrates how complex provocation experiments can be improved by sequential modeling and optimal design.