Population pharmacokinetic analysis of blood and joint synovial fluid concentrations of robenacoxib from healthy dogs and dogs with osteoarthritis.

PURPOSE: The purpose of this population analysis was to characterize the pharmacokinetic properties of robenacoxib in blood and stifle joint synovial fluid of dogs. METHODS: Data were obtained from two studies: 1) 8 healthy Beagle dogs in which an acute inflammation was induced by injection of urate crystals into one joint; 2) 95 dogs from various breeds diagnosed with osteoarthritis (OA). Robenacoxib concentrations in blood and synovial fluid were measured using a validated HPLC-UV and LC-MS method. Non-linear mixed effects modeling was performed using NONMEM6. RESULTS: A two-compartment pharmacokinetic model with linear elimination was developed to describe blood concentrations of robenacoxib. Blood clearance in healthy animals was found to be 75% higher than in OA dogs. Synovial fluid concentrations were modeled using an effect-compartment-type model predicting longer residence times in OA dogs compared to healthy Beagles (e.g. concentrations above the IC(50) for COX-2, respectively, 16 h vs. 10 h at 1.5 mg/kg). CONCLUSIONS: Robenacoxib was found to reside longer at the effect site (inflamed joint) compared to blood in both healthy and OA dogs. These results may explain the good efficacy observed with once-daily dosing in clinical trials and the high safety index of robenacoxib in dogs.

An integrated model for the glucose-insulin system.

The integrated glucose-insulin model was originally developed on a variety of intravenous glucose provocation experiments in healthy volunteers and type 2 diabetic patients. The model, which simultaneously describes time-courses of glucose and insulin based on mechanism-based components for production, elimination and homeostatic feedback, has been further extended to oral glucose provocations, meal tests and insulin administration. The model has been used to describe experiments ranging from 4 to 24 hr. Applications of the integrated glucose-insulin model include the clinical assessment of the mechanism of action of anti-diabetic drugs and the magnitude of their effects. Finally, the model was used for optimizing the design of provocation experiments.

An integrated glucose-insulin model to describe oral glucose tolerance test data in healthy volunteers.

The extension of the previously developed integrated models for glucose and insulin (IGI) to include the oral glucose tolerance test (OGTT) in healthy volunteers could be valuable to better understand the differences between healthy individuals and those with type 2 diabetes mellitus (T2DM). Data from an OGTT in 23 healthy volunteers were used. Analysis was based on the previously developed intravenous model with extensions for glucose absorption and incretin effect on insulin secretion. The need for additional structural components was evaluated. The model was evaluated by simulation and a bootstrap. Multiple glucose and insulin concentration peaks were observed in most individuals as well as hypoglycemic episodes in the second half of the experiment. The OGTT data were successfully described by the extended basic model. An additional control mechanism of insulin on glucose production improved the description of the data. The model showed good predictive properties, and parameters were estimated with good precision. In conclusion, a previously presented integrated model has been extended to describe glucose and insulin concentrations in healthy volunteers following an OGTT. The characterization of the differences between the healthy and diabetic stages in the IGI model could potentially be used to extrapolate drug effect from healthy volunteers to T2DM.

Optimization of the intravenous glucose tolerance test in T2DM patients using optimal experimental design.

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.

The impact of misspecification of residual error or correlation structure on the type I error rate for covariate inclusion.

It has been shown that when using the FOCE method in NONMEM, the likelihood ratio test (LRT) can be sensitive to the use of an inappropriate estimation method in that ignoring an existing eta-epsilon interaction leads to actual significance levels for type I errors being higher than the nominal levels. The objective of this study was to assess through simulations the LRT sensitivity to various types of residual error model misspecifications in both continuous and categorical data. The study contained two parts, simulations based on continuous and categorical data. Data sets containing 250 individuals with up to 24 observations per individual were simulated multiple times (1000) with different types of residual error models for the continuous data and different strength of correlation between observations for the categorical data. The data sets were analyzed using either the correct or a simpler (incorrect) model with or without addition of a covariate. The type I error rate of inclusion of the non-informative covariate on the 5% level was calculated as the number of runs where the drop in the objective function value (OFV) was larger than 3.84 when the covariate relationship was included in the model using the correct or the incorrect model. The difference in OFV between the model with the correct and the incorrect structure was also calculated as a measure of the residual error model misspecification. For continuous data the FOCE method was used in most cases (with interaction when appropriate). The Laplacian estimation method was used for one of the continuous models and for categorical data. The results showed that the residual error model misspecifications when the erroneous model was used were pronounced, as indicated by the OFV being substantially higher than for the corresponding correct models. The significance levels of the LRT with the incorrect model were appropriate in all cases but ignoring (serial) correlations between observations (continuous and categorical data) as well as when the eta-epsilon interaction was ignored (which has previously been shown, continuous data). When ignoring correlation, the type I error rates were shown to be sensitive to the correlation strength, the number of observations per individual and the magnitude of the inter-individual variability on clearance. We conclude that the LRT appears robust towards all tested cases, but ignoring (serial) correlations between observations and eta-epsilon interaction.

Approaches to handling pharmacodynamic baseline responses.

A few approaches for handling baseline responses are available for use in pharmacokinetic (PK)-pharmacodynamic (PD) analysis. They include: (method 1-B1) estimation of the typical value and interindividual variability (IIV) of baseline in the population, (B2) inclusion of the observed baseline response as a covariate acknowledging the residual variability, (B3) a more general version of B2 as it also takes the IIV of the baseline in the population into account, and (B4) normalization of all observations by the baseline value. The aim of this study was to investigate the relative performance of B1-B4. PD responses over a single dosing interval were simulated from an indirect response model in which a drug acts through stimulation or inhibition of the response according to an Emax model. The performance of B1-B4 was investigated under 22 designs, each containing 100 datasets. NONMEM VI beta was used to estimate model parameters with the FO and the FOCE method. The mean error (ME, %) and root mean squared error (RMSE, %) of the population parameter estimates were computed and used as an indicator of bias and imprecision. Absolute ME (|ME|) and RMSE from all methods were ranked within the same design, the lower the rank value the better method performance. Average rank of each method from all designs was reported. The results showed that with B1 and FOCE, the average of |ME| and RMSE across all typical individual parameters and all conditions was 5.9 and 31.8%. The average rank of |ME| for B1, B2, B3, and B4 was 3.7, 3.8, 3.3, and 5.2 for the FOCE method, and 4.6, 4.3, 4.7, and 6.4 for the FO method. The smallest imprecision was noted with the use of B1 (rank of 3.1 for FO, and 2.9 for FOCE) and increased, in order, with B3 (3.9-FO and 3.6-FOCE), B2 (4.8-FO; 4.7-FOCE), and B4 (6.4-FO; 6.5-FOCE). We conclude that when considering both bias and imprecision B1 was slightly better than B3 which in turn was better than B2. Differences between these methods were small. B4 was clearly inferior. The FOCE method led to a smaller bias, but no marked reduction in imprecision of parameter estimates compared to the FO method.

An integrated glucose-insulin model to describe oral glucose tolerance test data in type 2 diabetics.

An integrated model for the glucose-insulin system describing oral glucose tolerance test data was developed, extending on a previously introduced model for intravenous glucose provocations. Model extensions comprised the description of glucose absorption by a chain of transit compartments with a mean transit time of 35 minutes, a bioavailability of 80%, and a representation of the incretin effect, expressed as a direct effect of the glucose absorption rate on insulin secretion. The ability of the model to predict the incretin effect was assessed by simulating the observed difference in insulin response following an oral glucose tolerance test compared with an isoglycemic glucose infusion mimicking an oral glucose tolerance test profile. The extension of the integrated glucose-insulin model to gain information from oral glucose tolerance test data considerably expands its range of applications because the oral glucose tolerance test is one of the most common glucose challenge experiments for assessing the efficacy of hypoglycemic agents in clinical drug development.

An integrated model for glucose and insulin regulation in healthy volunteers and type 2 diabetic patients following intravenous glucose provocations.

An integrated model for the regulation of glucose and insulin concentrations following intravenous glucose provocations in healthy volunteers and type 2 diabetic patients was developed. Data from 72 individuals were included. Total glucose, labeled glucose, and insulin concentrations were determined. Simultaneous analysis of all data by nonlinear mixed effect modeling was performed in NONMEM. Integrated models for glucose, labeled glucose, and insulin were developed. Control mechanisms for regulation of glucose production, insulin secretion, and glucose uptake were incorporated. Physiologically relevant differences between healthy volunteers and patients were identified in the regulation of glucose production, elimination rate of glucose, and secretion of insulin. The model was able to describe the insulin and glucose profiles well and also showed a good ability to simulate data. The features of the present model are likely to be of interest for analysis of data collected in antidiabetic drug development and for optimization of study design.