An algorithm for estimating the terminal half-life in pharmacokinetic studies

Pharmacokinetics and tissue distribution of vigabatrin enantiomers in rats

Purpose

To investigate the pharmacokinetics and tissue distribution of VGB racemate and its single enantiomers, and explore the potential of clinic development for single enantiomer S-VGB.

Methods

In the pharmacokinetics study, male Sprague-Dawley rats were gavaged with VGB racemate or its single enantiomers dosing 50, 100 or 200 mg/kg, and the blood samples were collected during 12 h at regular intervals. In the experiment of tissue distribution, VGB and its single enantiomers were administered intravenously dosing 200 mg/kg, and the tissues including heart, liver, spleen, lung and kidney, eyes, hippocampus, and prefrontal cortex were separated at different times. The concentrations of R-VGB and S-VGB in the plasma and tissues were measured using HPLC.

Results

Both S-VGB and R-VGB could be detected in the plasma of rats administered with VGB racemate, reaching Cmax at approximately 0.5 h with t1/2 2-3 h. There was no significant pharmacokinetic difference between the two enantiomers when VGB racemate was given 200 mg/kg and 100 mg/kg. However, when given at the dose of 50 mg/kg, S-VGB presented a shorter t1/2 and a higher Cl/F than R-VGB, indicating a faster metabolism of S-VGB. Furthermore, when single enantiomer was administered respectively, S-VGB presented a slower metabolism than R-VGB, as indicated by a longer t1/2 and MRT but a lower Cmax. Moreover, compared with the VGB racemate, the single enantiomers S-VGB and R-VGB had shorter t1/2 and MRT, higher Cmax and AUC/D, and lower Vz/F and Cl/F, indicating the stronger oral absorption and faster metabolism of single enantiomer. In addition, regardless of VGB racemate administration or single enantiomer administration, S-VGB and R-VGB had similar characteristics in tissue distribution, and the content of S-VGB in hippocampus, prefrontal cortex and liver was much higher than that of R-VGB.

Conclusions

Although there is no transformation between S-VGB and R-VGB in vivo, those two enantiomers display certain disparities in the pharmacokinetics and tissue distribution, and interact with each other. These findings might be a possible interpretation for the pharmacological and toxic effects of VGB and a potential direction for the development and optimization of the single enantiomer S-VGB.

Pharmacokinetic parameters estimation using adaptive Bayesian P-splines models

In preclinical and clinical experiments, pharmacokinetic (PK) studies are designed to analyse the evolution of drug concentration in plasma over time i.e. the PK profile. Some PK parameters are estimated in order to summarize the complete drug's kinetic profile: area under the curve (AUC), maximal concentration (C(max)), time at which the maximal concentration occurs (t(max)) and half-life time (t(1/2)).Several methods have been proposed to estimate these PK parameters. A first method relies on interpolating between observed concentrations. The interpolation method is often chosen linear. This method is simple and fast. Another method relies on compartmental modelling. In this case, nonlinear methods are used to estimate parameters of a chosen compartmental model. This method provides generally good results. However, if the data are sparse and noisy, two difficulties can arise with this method. The first one is related to the choice of the suitable compartmental model given the small number of data available in preclinical experiment for instance. Second, nonlinear methods can fail to converge. Much work has been done recently to circumvent these problems (J. Pharmacokinet. Pharmacodyn. 2007; 34:229-249, Stat. Comput., to appear, Biometrical J., to appear, ESAIM P&S 2004; 8:115-131).In this paper, we propose a Bayesian nonparametric model based on P-splines. This method provides good PK parameters estimation, whatever be the number of available observations and the level of noise in the data. Simulations show that the proposed method provides better PK parameters estimations than the interpolation method, both in terms of bias and precision. The Bayesian nonparametric method provides also better AUC and t(1/2) estimations than a correctly specified compartmental model, whereas this last method performs better in t(max) and C(max) estimations.We extend the basic model to a hierarchical one that treats the case where we have concentrations from different subjects. We are then able to get individual PK parameter estimations. Finally, with Bayesian methods, we can get easily some uncertainty measures by obtaining credibility sets for each PK parameter.

Proposal for a standardised identification of the mono-exponential terminal phase for orally administered drugs

The area under the plasma concentration-time curve from time zero to infinity (AUC(0-inf)) is generally considered to be the most appropriate measure of total drug exposure for bioavailability/bioequivalence studies of orally administered drugs. However, the lack of a standardised method for identifying the mono-exponential terminal phase of the concentration-time curve causes variability for the estimated AUC(0-inf). The present investigation introduces a simple method, called the two times t(max) method (TTT method) to reliably identify the mono-exponential terminal phase in the case of oral administration. The new method was tested by Monte Carlo simulation in Excel and compared with the adjusted r squared algorithm (ARS algorithm) frequently used in pharmacokinetic software programs. Statistical diagnostics of three different scenarios, each with 10,000 hypothetical patients showed that the new method provided unbiased average AUC(0-inf) estimates for orally administered drugs with a monophasic concentration-time curve post maximum concentration. In addition, the TTT method generally provided more precise estimates for AUC(0-inf) compared with the ARS algorithm. It was concluded that the TTT method is a most reasonable tool to be used as a standardised method in pharmacokinetic analysis especially bioequivalence studies to reliably identify the mono-exponential terminal phase for orally administered drugs showing a monophasic concentration-time profile.

Pharmacokinetics, safety and bioavailability of HPMPC (cidofovir) in human immunodeficiency virus-infected subjects

We conducted a single-center, double-blind, placebo-controlled phase I study in HIV-positive subjects to ascertain the safety, tolerance, bioavailability, pharmacokinetics, and maximum tolerated dose of HPMPC (cidofovir). Five subjects were randomized to receive drug and two to receive placebo at each of three dosage tier (1, 3, and 10 mg/kg) with a 2-week washout period doses. Subjects at 1 and 3 mg/kg received single doses of HPMPC by subcutaneous (s.c.) intravenous (i.v.), and oral (p.o.) routes, while subjects at 10 mg/kg received only i.v. and p.o. doses. For subjects already taking zidovudine, zidovudine AUC values are determined before and then with HPMPC administration for each route. The AUC values of HPMPC were dose-proportional. Subcutaneous bioavailability was essentially equivalent to that of the intravenous route, but the development of transient local fibrosis ad the volumes needed for subcutaneous dosing precluded higher subcutaneous dosing than 3 mg/kg. Oral bioavailability was poor, estimated to be less than 5%. Drug elimination was predominantly renal. Nephrotoxicity in one subject was the only serious adverse event observed. This subject had a significant lag period prior to oral absorption and also had the highest AUC values for both HPMPC and zidovudine. We found no consistent effect on zidovudine AUC concomitant HPMPC.

An interactive computer program for determining areas bounded by drug concentration curves using lagrange interpolation

The LAGRAN method of Rocci and Jusko for determining the area of plasma concentration curves has been implemented in a user-friendly form. Drug concentration versus time data may be entered using the keyboard or imported in the form of simple text files from spreadsheets or other software. The MSDOS program allows prompt graphic observation of the data. The effect of selecting different fitting modes for each segment of the curve may be viewed interactively using this graphic display. Pharmacokinetic parameters that are provided by the program include mean residence time, variance of the residence time, plasma clearance, and steady-state volume of distribution.