Simulation-based systematic review of imatinib population pharmacokinetics and PK-PD relationships in chronic myeloid leukemia (CML) patients

Objectives: Several population pharmacokinetic (PPK) and pharmacokinetic-pharmacodynamic (PK-PD) analyses have been performed with the anticancer drug imatinib. Inspired by the approach of meta-analysis, we aimed to compare and combine results from published studies in a useful way - in particular for improving the clinical interpretation of imatinib concentration measurements in the scope of therapeutic drug monitoring (TDM). Methods: Original PPK analyses and PK-PD studies (PK surrogate: trough concentration Cmin; PD outcomes: optimal early response and specific adverse events) were searched systematically on MEDLINE. From each identified PPK model, a predicted concentration distribution under standard dosage was derived through 1000 simulations (NONMEM), after standardizing model parameters to common covariates. A "reference range" was calculated from pooled simulated concentrations in a semi-quantitative approach (without specific weighting) over the whole dosing interval. Meta-regression summarized relationships between Cmin and optimal/suboptimal early treatment response. Results: 9 PPK models and 6 relevant PK-PD reports in CML patients were identified. Model-based predicted median Cmin ranged from 555 to 1388 ng/ml (grand median: 870 ng/ml and inter-quartile range: 520-1390 ng/ml). The probability to achieve optimal early response was predicted to increase from 60 to 85% from 520 to 1390 ng/ml across PK-PD studies (odds ratio for doubling Cmin: 2.7). Reporting of specific adverse events was too heterogeneous to perform a regression analysis. The general frequency of anemia, rash and fluid retention increased however consistently with Cmin, but less than response probability. Conclusions: Predicted drug exposure may differ substantially between various PPK analyses. In this review, heterogeneity was mainly attributed to 2 "outlying" models. The established reference range seems to cover the range where both good efficacy and acceptable tolerance are expected for most patients. TDM guided dose adjustment appears therefore justified for imatinib in CML patients. Its usefulness remains now to be prospectively validated in a randomized trial.

Population pharmacokinetics in chronic myeloid leukaemia patients: observations under field-conditions

Introduction: Therapeutic drug monitoring (TDM) of imatinib has been increasingly proposed for chronic myeloid leukaemia (CML) patients, as several studies have found a correlation between trough concentrations (Cmin) >=1000ng/ml and improved response. The pharmacological monitoring project of EUTOS (European Treatment and Outcome Study) was launched to increase the availability of imatinib TDM, standardize labs, and validate proposed Cmin thresholds. Using the collected data, the objective of this analysis was to characterize imatinib Population pharmacokinetics (Pop-PK) in a large cohort of European patients, to quantify its variability and the influence of demographic factors and comedications, and to derive individual exposure variables suitable for further concentration-effect analyses.¦Methods: 4095 PK samples from 2478 adult patients were analyzed between 2006 and 2010 by LC-MS-MS and considered for Pop-PK analysis by NONMEM®. Model building used data from 973 patients with >=2 samples available (2590 samples). A sensitivity analysis was performed using all data. Available comedications (27%) were classified into inducers or inhibitors of P-glycoprotein, CYP3A4/5 and organic-cation-transporter-1 (hOCT-1).¦Results: A one-compartment model with linear elimination, zero-order absorption fitted the data best. Estimated Pop-PK parameters (interindividual variability, IIV %CV) for a 40-year old male patient were: clearance CL = 17.3 L/h (37.7%), volume V = 429L (51.1%), duration of absorption D1 = 3.2h. Outliers, reflecting potential compliance and time recording errors, were taken into account by estimating an IIV on the residual error (35.4%). Intra-individual residuals were 29.1% (proportional) plus ± 84.6 ng/mL (additive). Female patients had a 15.2% lower CL (14.6 L/h). A piece-wise linear effect of age estimated a CL of 18.7 L/h at 20 years, 17.3 L/h at 40 and 13.8 L/h at 60 years. These covariates explained 2% (CL) and 4.5% (V) of IIV variability. No effect of comedication was found. The sensitivity analysis expectedly estimated increased IIV, but similar fixed effect parameters.¦Conclusion: Imatinib PK was well described in a large cohort of CML patients under field conditions and results were concordant with previous studies. Patient characteristics explain only little IIV, confirming limited utility of prior dosage adjustment. As intra-variability is smaller than inter-patient variability, dose adjustment guided by TDM could however be beneficial in order to bring Cmin into a given therapeutic target.

Limited sampling strategies for estimation of cyclosporine exposure in pediatric hematopoietic stem cell transplant recipients : methodological improvement and introduction of sampling time deviation analysis.

Le suivi thérapeutique est recommandé pour l’ajustement de la dose des agents immunosuppresseurs. La pertinence de l’utilisation de la surface sous la courbe (SSC) comme biomarqueur dans l’exercice du suivi thérapeutique de la cyclosporine (CsA) dans la transplantation des cellules souches hématopoïétiques est soutenue par un nombre croissant d’études. Cependant, pour des raisons intrinsèques à la méthode de calcul de la SSC, son utilisation en milieu clinique n’est pas pratique. Les stratégies d’échantillonnage limitées, basées sur des approches de régression (R-LSS) ou des approches Bayésiennes (B-LSS), représentent des alternatives pratiques pour une estimation satisfaisante de la SSC. Cependant, pour une application efficace de ces méthodologies, leur conception doit accommoder la réalité clinique, notamment en requérant un nombre minimal de concentrations échelonnées sur une courte durée d’échantillonnage. De plus, une attention particulière devrait être accordée à assurer leur développement et validation adéquates. Il est aussi important de mentionner que l’irrégularité dans le temps de la collecte des échantillons sanguins peut avoir un impact non-négligeable sur la performance prédictive des R-LSS. Or, à ce jour, cet impact n’a fait l’objet d’aucune étude. Cette thèse de doctorat se penche sur ces problématiques afin de permettre une estimation précise et pratique de la SSC. Ces études ont été effectuées dans le cadre de l’utilisation de la CsA chez des patients pédiatriques ayant subi une greffe de cellules souches hématopoïétiques. D’abord, des approches de régression multiple ainsi que d’analyse pharmacocinétique de population (Pop-PK) ont été utilisées de façon constructive afin de développer et de valider adéquatement des LSS. Ensuite, plusieurs modèles Pop-PK ont été évalués, tout en gardant à l’esprit leur utilisation prévue dans le contexte de l’estimation de la SSC. Aussi, la performance des B-LSS ciblant différentes versions de SSC a également été étudiée. Enfin, l’impact des écarts entre les temps d’échantillonnage sanguins réels et les temps nominaux planifiés, sur la performance de prédiction des R-LSS a été quantifié en utilisant une approche de simulation qui considère des scénarios diversifiés et réalistes représentant des erreurs potentielles dans la cédule des échantillons sanguins. Ainsi, cette étude a d’abord conduit au développement de R-LSS et B-LSS ayant une performance clinique satisfaisante, et qui sont pratiques puisqu’elles impliquent 4 points d’échantillonnage ou moins obtenus dans les 4 heures post-dose. Une fois l’analyse Pop-PK effectuée, un modèle structural à deux compartiments avec un temps de délai a été retenu. Cependant, le modèle final - notamment avec covariables - n’a pas amélioré la performance des B-LSS comparativement aux modèles structuraux (sans covariables). En outre, nous avons démontré que les B-LSS exhibent une meilleure performance pour la SSC dérivée des concentrations simulées qui excluent les erreurs résiduelles, que nous avons nommée « underlying AUC », comparée à la SSC observée qui est directement calculée à partir des concentrations mesurées. Enfin, nos résultats ont prouvé que l’irrégularité des temps de la collecte des échantillons sanguins a un impact important sur la performance prédictive des R-LSS; cet impact est en fonction du nombre des échantillons requis, mais encore davantage en fonction de la durée du processus d’échantillonnage impliqué. Nous avons aussi mis en évidence que les erreurs d’échantillonnage commises aux moments où la concentration change rapidement sont celles qui affectent le plus le pouvoir prédictif des R-LSS. Plus intéressant, nous avons mis en exergue que même si différentes R-LSS peuvent avoir des performances similaires lorsque basées sur des temps nominaux, leurs tolérances aux erreurs des temps d’échantillonnage peuvent largement différer. En fait, une considération adéquate de l'impact de ces erreurs peut conduire à une sélection et une utilisation plus fiables des R-LSS. Par une investigation approfondie de différents aspects sous-jacents aux stratégies d’échantillonnages limités, cette thèse a pu fournir des améliorations méthodologiques notables, et proposer de nouvelles voies pour assurer leur utilisation de façon fiable et informée, tout en favorisant leur adéquation à la pratique clinique.

Population pharmacokinetics and clinical response for artemether-lumefantrine in pregnant and nonpregnant women with uncomplicated Plasmodium falciparum malaria in Tanzania.

Artemether-lumefantrine (AL) is the first-line treatment for uncomplicated malaria in the second and third trimesters of pregnancy. Its efficacy during pregnancy has recently been challenged due to altered pharmacokinetic (PK) properties in this vulnerable group. The aim of this study was to determine the PK profile of AL in pregnant and nonpregnant women and assess their therapeutic outcome. Thirty-three pregnant women and 22 nonpregnant women with malaria were treated with AL (80/480 mg) twice daily for 3 days. All patients provided five venous plasma samples for drug quantification at random times over 7 days. Inter- and intraindividual variability was assessed, and the effects of covariates were quantified using a nonlinear mixed-effects modeling approach (NONMEM). A one-compartment model with first-order absorption and elimination with linear metabolism from drug to metabolite fitted the data best for both arthemether (AM) and lumefantrine (LF) and their metabolites. Pregnancy status and diarrhea showed a significant influence on LF PK. The relative bioavailability of lumefantrine and its metabolism rate into desmethyl-lumefantrine were, respectively, 34% lower and 78% higher in pregnant women than in nonpregnant patients. The overall PCR-uncorrected treatment failure rates were 18% in pregnant women and 5% in nonpregnant women (odds ratio [OR] = 4.04; P value of 0.22). A high median day 7 lumefantrine concentration was significantly associated with adequate clinical and parasitological response (P = 0.03). The observed reduction in the relative bioavailability of lumefantrine in pregnant women may explain the higher treatment failure in this group, mostly due to lower posttreatment prophylaxis. Hence, a modified treatment regimen of malaria in pregnancy should be considered.

Population pharmacokinetics of teicoplanin in children.

Teicoplanin is frequently administered to treat Gram-positive infections in pediatric patients. However, not enough is known about the pharmacokinetics (PK) of teicoplanin in children to justify the optimal dosing regimen. The aim of this study was to determine the population PK of teicoplanin in children and evaluate the current dosage regimens. A PK hospital-based study was conducted. Current dosage recommendations were used for children up to 16 years of age. Thirty-nine children were recruited. Serum samples were collected at the first dose interval (1, 3, 6, and 24 h) and at steady state. A standard 2-compartment PK model was developed, followed by structural models that incorporated weight. Weight was allowed to affect clearance (CL) using linear and allometric scaling terms. The linear model best accounted for the observed data and was subsequently chosen for Monte Carlo simulations. The PK parameter medians/means (standard deviation [SD]) were as follows: CL, [0.019/0.023 (0.01)] × weight liters/h/kg of body weight; volume, 2.282/4.138 liters (4.14 liters); first-order rate constant from the central to peripheral compartment (Kcp), 0.474/3.876 h(-1) (8.16 h(-1)); and first-order rate constant from peripheral to central compartment (Kpc), 0.292/3.994 h(-1) (8.93 h(-1)). The percentage of patients with a minimum concentration of drug in serum (Cmin) of <10 mg/liter was 53.85%. The median/mean (SD) total population area under the concentration-time curve (AUC) was 619/527.05 mg · h/liter (166.03 mg · h/liter). Based on Monte Carlo simulations, only 30.04% (median AUC, 507.04 mg · h/liter), 44.88% (494.1 mg · h/liter), and 60.54% (452.03 mg · h/liter) of patients weighing 50, 25, and 10 kg, respectively, attained trough concentrations of >10 mg/liter by day 4 of treatment. The teicoplanin population PK is highly variable in children, with a wider AUC distribution spread than for adults. Therapeutic drug monitoring should be a routine requirement to minimize suboptimal concentrations. (This trial has been registered in the European Clinical Trials Database Registry [EudraCT] under registration number 2012-005738-12.).

Population pharmacokinetics of hydroxychloroquine in patients with rheumatoid arthritis

Hydroxychloroquine (HCQ) is an antimalarial drug that is also used as a second-line treatment of rheumatoid arthritis (RA). Clinically, the use of HCQ is characterized by a long delay in the onset of action, and withdrawal of treatment is often a result of inefficacy rather than from toxicity. The slow onset of action can be attributed to the pharmacokinetics (PK) of HCQ, and wide interpatient variability is evident. Tentative relationships between concentration and effect have been made, but to date, no population PK model has been developed for HCQ. This study aimed to develop a population PK model including an estimation of the oral bioavailability of HCQ. In addition, the effects of the coadministration of methotrexate on the PK of HCQ were examined. Hydroxychloroquine blood concentration data were combined from previous pharmacokinetic studies in patients with rheumatoid arthritis. A total of 123 patients were studied, giving the data cohort from four previously published studies. Two groups of patients were included: 74 received hydroxychloroquine (HCQ) alone, and 49 received HCQ and methotrexate (MTX). All data analyses were carried out using the NONMEM program. A one-compartment PK model was supported, rather than a three-compartment model as previously published, probably because of the clustering of concentrations taken at the end of a dosing interval. The population estimate of bioavailability of 0.75 (0.07), n = 9, was consistent with literature values. The parameter values from the final model were: (Cl) over bar = 9.9 +/- 0.4 L/h, (V) over bar 605 +/- 91 L, (k(d)) over bar = 0.77 +/- 0.22 hours(-1), (t(tag)) over bar = 0.44 +/- 0.02 hours. Clearance was not affected by the presence of MTX, and, hence, steady-state drug concentrations and maintenance dosage requirements were similar. A population PK model was successfully developed for HCQ.

Optimal design for model discrimination and parameter estimation for itraconazole population pharmacokinetics in cystic fibrosis patients

Optimal sampling times are found for a study in which one of the primary purposes is to develop a model of the pharmacokinetics of itraconazole in patients with cystic fibrosis for both capsule and solution doses. The optimal design is expected to produce reliable estimates of population parameters for two different structural PK models. Data collected at these sampling times are also expected to provide the researchers with sufficient information to reasonably discriminate between the two competing structural models.

Paediatric population pharmacokinetics of itraconazole and its active metabolite hydroxy-itraconazole in cystic fibrosis and bone marrow transplant patients

Objectives: The aim of the study was to characterise the population pharmacokinetics (popPK) properties of itraconazole (ITRA) and its active metabolite hydroxy-ITRA in a representative paediatric population of cystic fibrosis (CF) and bone marrow transplant (BMT) patients. The goals were to determine the relative bioavailability between the two oral formulations, and to explore improved dosage regimens in these patients. Methods: All paediatric patients with CF taking oral ITRA for the treatment of allergic bronchopulmonary aspergillosis and patients undergoing BMT who were taking ITRA for prophylaxis of any fungal infection were eligible for the study. A minimum of two blood samples were drawn after the capsules and also after switching to oral solution, or vice versa. ITRA and hydroxy-ITRA plasma concentrations were measured by HPLC[1]. A nonlinear mixed-effect modelling approach (NONMEM 5.1.1) was used to describe the PK of ITRA and hydroxy-ITRA simultaneously. Simulations were used to assess dosing strategies in these patients. Results: Forty-nine patients (29CF, 20 BMT) were recruited to the study who provided 227 blood samples for the population analysis. A 1-compartment model with 1st order absorption and elimination best described ITRA kinetics, with 1st order conversion to hydroxy-ITRA. For ITRA, the apparent clearance (ClItra/F) and volume of distribution (Vitra/F) was 35.5L/h and 672L, respectively; the absorption rate constant for the capsule formulation was 0.0901 h-1 and for the oral solution formulation it was 0.959 h-1. The capsule comparative bioavailability (vs. solution) was 0.55. For hydroxy-ITRA, the apparent volume of distribution and clearance were 10.6 L and 5.28 L/h, respectively. Of several screened covariates only allometrically scaled total body weight significantly improved the fit to the data. No difference between the two populations was found. Conclusion: The developed popPK model adequately described the pharmacokinetics of ITRA and hydroxy-ITRA in paediatric patients with CF and patients undergoing BMT. High inter-patient variability confirmed previous data in CF[2], leukaemia and BMT[3] patients. From the population model, simulations showed the standard dose (5 mg/kg/day) needs to be doubled for the solution formulation and even 4 times more given of the capsules to achieve an adequate target therapeutic trough plasma concentration of 0.5 mg/L[4] in these patients.

Population pharmacokinetics of single-dose intravenous paracetamol in children

Background To determine the pharmacokinetics (PK) of a new i.v. formulation of paracetamol (Perfalgan) in children ≤15 yr of age. Methods After obtaining written informed consent, children under 16 yr of age were recruited to this study. Blood samples were obtained at 0, 15, 30 min, 1, 2, 4, 6, and 8 h after administration of a weight-dependent dose of i.v. paracetamol. Paracetamol concentration was measured using a validated high-performance liquid chromatographic assay with ultraviolet detection method, with a lower limit of quantification (LLOQ) of 900 pg on column and an intra-day coefficient of variation of 14.3% at the LLOQ. Population PK analysis was performed by non-linear mixed-effect modelling using NONMEM. Results One hundred and fifty-nine blood samples from 33 children aged 1.8–15 yr, weight 13.7–56 kg, were analysed. Data were best described by a two-compartment model. Only body weight as a covariate significantly improved the goodness of fit of the model. The final population models for paracetamol clearance (CL), V1 (central volume of distribution), Q (inter-compartmental clearance), and V2 (peripheral volume of distribution) were: 16.51×(WT/70)0.75, 28.4×(WT/70), 11.32×(WT/70)0.75, and 13.26×(WT/70), respectively (CL, Q in litres per hour, WT in kilograms, and V1 and V2 in litres). Conclusions In children aged 1.8–15 yr, the PK parameters for i.v. paracetamol were not influenced directly by age but were by total body weight and, using allometric size scaling, significantly affected the clearances (CL, Q) and volumes of distribution (V1, V2).

The enantioselective population pharmacokinetics of intravenous ketorolac in children using a stereoselective assay suitable for microanalysis

Objective: To describe the effect of age and body size on enantiomer selective pharmacokinetic (PK) of intravenous ketorolac in children using a microanalytical assay. Methods: Blood samples were obtained at 0, 15 and 30 min and at 1, 2, 4, 6, 8 and 12 h after a weight-dependent dose of ketorolac. Enantiomer concentration was measured using a liquid chromatography tandem mass spectrometry method. Non-linear mixed-effect modelling was used to assess PK parameters. Key findings: Data from 11 children (1.7–15.6 years, weight 10.7–67.4 kg) were best described by a two-compartment model for R(+), S(−) and racemic ketorolac. Only weight (WT) significantly improved the goodness of fit. The final population models were CL = 1.5 × (WT/46)0.75, V1 = 8.2 × (WT/46), Q = 3.4 × (WT/46)0.75, V2 = 7.9 × (WT/46), CL = 2.98 × (WT/46), V1 = 13.2 × (WT/46), Q = 2.8 × (WT/46)0.75, V2 = 51.5 × (WT/46), and CL = 1.1 × (WT/46)0.75, V1 = 4.9 × (WT/46), Q = 1.7 × (WT/46)0.75 and V2 = 6.3 × (WT/46)for R(+), S(−) and racemic ketorolac. Conclusions: Only body weight influenced the PK parameters for R(+) and S(−) ketorolac. Using allometric size scaling significantly affected the clearances (CL, Q) and volumes of distribution (V1, V2).

Population pharmacokinetics and pharmacodynamics: An underutilized resource

A meeting was convened in Canberra, Australia, at the request of the Australian Drug Evaluation Committee (ADEC), on December 3-4, 1997 to discuss the role of population pharmacokinetics and pharmacodynamics in drug evaluation and development. The ADEC was particularly concerned about registration of drugs in the pediatric age group. The population approach could be used more often than is currently the case in pharmacokinetic and pharmacodynamic studies to provide valuable information for the safe and effective use of drugs in neonates, infants, and children. The meeting ultimately broadened to include discussion about other subgroups. The main conclusions of the meeting were: 1. The population approach, pharmacokinetic and pharmacodynamic analysis, is a valuable tool both for drug registration purposes and for optimal dosing of drugs in specific groups of patients, 2. Population pharmacokinetic and pharmacodynamic studies are able to fill in the gaps' in registration of drugs, for example, to provide information on optimal pediatric dosing. Such studies provide a basis for enhancing product information to improve rational prescribing, 3. Expertise is required to perform the population studies and expertise, with a clinical perspective, is also required to evaluate such studies if they are to be submitted as part of a drug registration dossier Such expertise is available in the Australasian region and is increasing. Centers of excellence with the appropriate expertise to advise and assist should be encouraged to develop and grow in the region, 4. The use of the population approach by the pharmaceutical industry needs to be encouraged to provide valuable information not obtainable by other techniques. The acceptance of population pharmacokinetic and pharmacodynamic analyses by regulatory agencies also needs to be encouraged, and 5. Development of the population approach to pharmacokinetics and pharmacodynamics is needed from a public health perspective to ensure that all available information is collected and used to improve the way drugs are used. This important endeavor needs funding and support at the local and international levels.

Population pharmacokinetics of perhexiline from very sparse, routine monitoring data

Using NONMEM, the population pharmacokinetics of perhexiline were studied in 88 patients (34 F, 54 M) who were being treated for refractory angina. Their mean +/- SD (range) age was 75 +/- 9.9 years (46-92), and the length of perhexiline treatment was 56 +/- 77 weeks (0.3-416). The sampling time after a dose was 14.1 +/- 21.4 hours (0.5-200), and the perhexiline plasma concentrations were 0.39 +/- 0.32 mg/L (0.03-1.56). A one-compartment model with first-order absorption was fitted to the data using the first-order (FO) approximation. The best model contained 2 subpopulations (obtained via the $MIXTURE subroutine) of 77 subjects (subgroup A) and 11 subjects (subgroup B) that had typical values for clearance (CL/F) of 21.8 L/h and 2.06 L/h, respectively. The volumes of distribution (V/F) were 1470 L and 260 L, respectively, which suggested a reduction in presystemic metabolism in subgroup B. The interindividual variability (CV%) was modeled logarithmically and for CL/F ranged from 69.1% (subgroup A) to 86.3% (subgroup B). The interindividual variability in V/F was 111%. The residual variability unexplained by the population model was 28.2%. These results confirm and extend the existing pharmacokinetic data on perhexiline, especially the bimodal distribution of CL/F manifested via an inherited deficiency in hepatic and extrahepatic CYP2D6 activity.

Population pharmacokinetics of magnesium in preeclampsia

OBJECTIVE: The purpose of this study was to determine the population pharmacokinetics of magnesium from sparse observational data in patients with preeclampsia. STUDY DESIGN: Serum magnesium concentrations (1-11 per patient) were obtained retrospectively from the records of 116 patients with preeclampsia who had a loading dose of magnesium sulfate (16 or 20 mmol), followed by a maintenance dose (1 mmol/h) over an average of 28 hours. Population clearance, volume of distribution, and the baseline magnesium concentration were estimated using the NONMEM program. RESULTS: The following population typical values, together with the interpatient variability (expressed as coefficient of variation) were obtained with the use of a 1-compartment model: systemic clearance, 4.28 L/h (37.3%); volume of distribution, 32.3 L (32.1%); baseline concentration, 0.811 mmol/L (18.5%). The average half-life was 5.2 hours. Clonus was not obtunded in 4 patients whose serum magnesium concentrations were similar to the average concentration of 1.7 mmol/L. The variability remaining unexplained after the population model was fitted to the data was 6.5% to 10.8%. CONCLUSION: This study extended knowledge of the pharmacokinetic disposition of magnesium in preeclampsia. The results are potentially useful for the calculation of loading and maintenance doses, particularly when the relationship between serum concentration and effect in preeclampsia is clarified.

Population pharmacokinetics of lamotrigine

The present study estimated the population pharmacokinetics of lamotrigine in patients receiving oral lamotrigine therapy with drug concentration monitoring, and determined intersubject and intrasubject variability. A total of 129 patients were analyzed from two clinical sites. Of these, 124 patients provided spare data (198 concentration-time points); nine patients (four from a previous group plus five from the current group) provided rich data (431 points). The population analysis was conducted using P-PHARM (TM) (SIMED Scientific Software, Cedex, France), a nonlinear mixed-effect modeling program. A single exponential elimination model (first-order absorption) with heteroscedastic weighting was used. Apparent clearance (CL/F) and volume of distribution (V/F) were the pharmacokinetic parameters estimated. Covariate analysis was performed to determine which factors explained any of the variability associated with lamotrigine clearance. Population estimates of CL/F and V/F for lamotrigine generated in the final model were 2.14 +/- 0.81 L/h and 78.1 +/- 5.1 L/kg. Intersubject and intrasubject variability for clearance was 38% and 38%, respectively. The covariates of concomitant valproate and phenytoin therapy accounted for 42% of the intersubject variability of clearance. Age, gender, clinic site, and other concomitant antiepileptic drugs did not influence clearance. This study of the population pharmacokinetics of lamotrigine in patients using the drug clinically provides useful data and should lead to better dosage individualization for lamotrigine.

Population pharmacokinetics of tacrolimus in children who receive cut-down or full liver transplants

Background. The aim of this study was to investigate the population pharmacokinetics of tacrolimus in pediatric liver transplant recipients and to identify factors that may explain pharmacokinetic variability. Methods. Data were collected retrospectively from 35 children who received oral immunosuppressant therapy with tacrolimus. Maximum likelihood estimates were sought for the typical values of apparent clearance (CL/F) and apparent volume of distribution (V/F) with the program NONMEM. Factors screened for influence on the pharmacokinetic parameters were weight, age, gender, postoperative day, days since commencing tacrolimus therapy, transplant type (whole child liver or cut-down adult liver), liver function tests (bilirubin, alkaline phosphatase [ALP], aspartate aminotransferase [AST], gamma -glutamyl transferase [GGT], alanine aminotransferase [ALT]), creatinine clearance, hematocrit, corticosteroid dose, and concurrent therapy with metabolic inducers and inhibitors of tacrolimus. Results. No clear correlation existed between tacrolimus dosage and blood concentrations (r(2) =0.003). Transplant type, age, and liver function test values were the most important factors (P