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

Population pharmacokinetics of tacrolimus in adult kidney transplant recipients

Objectives: The aims of this study were to investigate the population pharmacokinetics of tacrolimus in adult kidney transplant recipients and to identify factors that explain variability. Methods: Population analysis was performed on retrospective data from 70 patients who received oral tacrolimus twice daily. Morning blood trough concentrations were measured by liquid chromatography-tandem mass spectrometry. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F), with the use of NONMEM (GloboMax LLC, Hanover, Md). Factors screened for influence on these parameters were weight, age, gender, postoperative day, days of tacrolimus therapy, liver function tests, creatinine clearance, hematocrit fraction, corticosteroid dose, and potential interacting drugs. Results. CL/F was greater in patients with abnormally low hematocrit fraction (data from 21 patients only), and it decreased with increasing days of therapy and AST concentrations (P

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.

Toward better outcomes with tacrolimus therapy: Population pharmacokinetics and individualized dosage prediction in adult liver transplantation

Patient outcomes in transplantation would improve if dosing of immunosuppressive agents was individualized. The aim of this study is to develop a population pharmacokinetic model of tacrolimus in adult liver transplant recipients and test this model in individualizing therapy. Population analysis was performed on data from 68 patients. Estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F) using the nonlinear mixed effects model program (NONMEM). Factors screened for influence on these parameters were weight, age, sex, transplant type, biliary reconstructive procedure, postoperative day, days of therapy, liver function test results, creatinine clearance, hematocrit, corticosteroid dose, and interacting drugs. The predictive performance of the developed model was evaluated through Bayesian forecasting in an independent cohort of 36 patients. No linear correlation existed between tacrolimus dosage and trough concentration (r(2) = 0.005). Mean individual Bayesian estimates for CL/F and V/F were 26.5 8.2 (SD) L/hr and 399 +/- 185 L, respectively. CL/F was greater in patients with normal liver function. V/F increased with patient weight. CL/F decreased with increasing hematocrit. Based on the derived model, a 70-kg patient with an aspartate aminotransferase (AST) level less than 70 U/L would require a tacrolimus dose of 4.7 mg twice daily to achieve a steady-state trough concentration of 10 ng/mL. A 50-kg patient with an AST level greater than 70 U/L would require a dose of 2.6 mg. Marked interindividual variability (43% to 93%) and residual random error (3.3 ng/mL) were observed. Predictions made using the final model were reasonably nonbiased (0.56 ng/mL), but imprecise (4.8 ng/mL). Pharmacokinetic information obtained will assist in tacrolimus dosing; however, further investigation into reasons for the pharmacokinetic variability of tacrolimus is required.

Population pharmacokinetics and association between A77 1726 plasma concentrations and disease activity measures following administration of leflunomide to people with rheumatoid arthritis

Aims To investigate the concentration-effect relationship and pharmacokinetics of leflunomide in patients with rheumatoid arthritis (RA). Methods Data were collected from 23 RA patients on leflunomide therapy (as sole disease modifying antirheumatic drug (DMARD)) for at least 3 months. Main measures were A77 1726 (active metabolite of leflunomide) plasma concentrations and disease activity measures including pain, duration/intensity of morning stiffness, and SF-36 survey. A population estimate was sought for apparent clearance (CL/F ) and volume of distribution was fixed (0.155 l kg(-1)). Factors screened for influence on CL/F were weight, age, gender and estimated creatinine clearance. Results Significantly higher A77 1726 concentrations were seen in patients with less swollen joints and with higher SF-36 mental summary scores than in those with measures indicating more active disease (P < 0.05); concentration-effect trends were seen with five other disease activity measures. Statistical analysis of all disease activity measures showed that mean A77 1726 concentrations in groups with greater control of disease activity were significantly higher than those in whom the measures indicated less desirable control (P < 0.05). There was large between subject variability in the dose-concentration relationship. A steady-state infusion model best described the pharmacokinetic data. Inclusion of age as a covariate decreased interindividual variability (P < 0.01), but this would not be clinically important in terms of dosage changes. Final parameter estimate (% CV interindividual variability) for CL/F was 0.0184 l h(-1) (50%) (95% CI 0.0146, 0.0222). Residual (unexplained) variability (% CV) was 8.5%. Conclusions This study of leflunomide in patients using the drug clinically indicated a concentration-effect relationship. From our data, a plasma A77 1726 concentration of 50 mg l(-1) is more likely to indicate someone with less active disease than is a concentration around 30 mg l(-1). The marked variability in pharmacokinetics suggests a place for individualized dosing of leflunomide in RA therapy.

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.

The population pharmacokinetics of citalopram after deliberate self-poisoning: A Bayesian approach

Defining the pharmacokinetics of drugs in overdose is complicated. Deliberate self-poisoning is generally impulsive and associated with poor accuracy in dose history. In addition, early blood samples are rarely collected to characterize the whole plasma-concentration time profile and the effect of decontamination on the pharmacokinetics is uncertain. The aim of this study was to explore a fully Bayesian methodology for population pharmacokinetic analysis of data that arose from deliberate self-poisoning with citalopram. Prior information on the pharmacokinetic parameters was elicited from 14 published studies on citalopram when taken in therapeutic doses. The data set included concentration-time data from 53 patients studied after 63 citalopram overdose events (dose range: 20-1700 mg). Activated charcoal was administered between 0.5 and 4 h after 17 overdose events. The clinical investigator graded the veracity of the patients' dosing history on a 5-point ordinal scale. Inclusion of informative priors stabilised the pharmacokinetic model and the population mean values could be estimated well. There were no indications of non-linear clearance after excessive doses. The final model included an estimated uncertainty of the dose amount which in a simulation study was shown to not affect the model's ability to characterise the effects of activated charcoal. The effect of activated charcoal on clearance and bioavailability was pronounced and resulted in a 72% increase and 22% decrease, respectively. These findings suggest charcoal administration is potentially beneficial after citalopram overdose. The methodology explored seems promising for exploring the dose-exposure relationship in the toxicological settings.

Population pharmacokinetics of metformin in late pregnancy

The pharmacokinetic disposition of metformin in late pregnancy was studied together with the level of fetal exposure at birth. Blood samples were obtained in the third trimester of pregnancy from women with gestational diabetes or type 2 diabetes, 5 had a previous diagnosis of polycystic ovary syndrome. A cord blood sample also was obtained at the delivery of some of these women, and also at delivery of others who had been taking metformin during pregnancy but from whom no blood had been taken. Plasma metformin concentrations were assayed by a new, validated, reverse-phase HPLC method, A 2-compartment, extravascular maternal model with transplacental partitioning of drug to a fetal compartment was fitted to the data. Nonlinear mixed-effects modeling was performed in'NONMEM using FOCE with INTERACTION. Variability was estimated using logarithmic interindividual and additive residual variance models; the covariance between clearance and volume was modeled simultaneously. Mean (range) metformin concentrations in cord plasma and in maternal plasma were 0.81 (range, 0.1-2.6) mg/L and 1.2 (range, 0. 1-2.9) mg/L, respectively. Typical population values (interindividual variability, CV%) for allometrically scaled maternal clearance and volume of distribution were 28 L/h/70 kg (17.1%) and 190 L/70 ka (46.3%), giving a derived population-wide half-life of 5.1 hours. The placental partition coefficient for metformin was 1.07 (36.3%). Neither maternal age nor weight significantly influenced the pharmacokinetics. The variability (SD) of observed concentrations about model-predicted concentrations was 0.32 mg/L. The pharmacokinetics were similar to those in nonpregnant patients and, therefore, no dosage adjustment is warranted. Metformin readily crosses the placenta, exposing the fetus to concentrations approaching those in the maternal circulation. The sequelae to such exposure, ea, effects on neonatal obesity and insulin resistance, remain unknown.

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

Objective: The objective of the study was to characterise the population pharmacokinetic properties of itraconazole and its active metabolite hydroxyitraconazole in a representative paediatric population of cystic fibrosis and bone marrow transplant (BMT) patients and to identify patient characteristics influencing the pharmacokinetics of itraconazole. The ultimate goals were to determine the relative bioavailability between the two oral formulations (capsules vs oral solution) and to optimise dosing regimens in these patients. Methods: All paediatric patients with cystic fibrosis or patients undergoing BMT at The Royal Children's Hospital, Brisbane, QLD, Australia, who were prescribed oral itraconazole for the treatment of allergic bronchopulmonary aspergillosis (cystic fibrosis patients) or for prophylaxis of any fungal infection (BMT patients) were eligible for the study. Blood samples were taken from the recruited patients as per an empirical sampling design either during hospitalisation or during outpatient clinic visits. ltraconazole and hydroxy-itraconazole plasma concentrations were determined by a validated high-performance liquid chromatography assay with fluorometric detection. A nonlinear mixed-effect modelling approach using the NONMEM software to simultaneously describe the pharmacokinetics of itraconazole and its metabolite. Results: A one-compartment model with first-order absorption described the itraconazole data, and the metabolism of the parent drug to hydroxy-itraconazole was described by a first-order rate constant. The metabolite data also showed one-compartment characteristics with linear elimination. For itraconazole the apparent clearance (CLitraconazole) was 35.5 L/hour, the apparent volume of distribution (V-d(itraconazole)) was 672L, the absorption rate constant for the capsule formulation was 0.0901 h(-1) and for the oral solution formulation was 0.96 h-1. The lag time was estimated to be 19.1 minutes and the relative bioavailability between capsules and oral solution (F-rel) was 0.55. For the metabolite, volume of distribution, V-m/(F (.) f(m)), and clearance, CL/(F (.) fm), were 10.6L and 5.28 L/h, respectively. The influence of total bodyweight was significant, added as a covariate on CLitraconazoie/F and V-d(itraconazole)/F (standardised to a 70kg person) using allometric three-quarter power scaling on CLitraconazole/F, which therefore reflected adult values. The unexplained between-subject variability (coefficient of variation %) was 68.7%, 75.8%, 73.4% and 61.1% for CLitraconazoie/F, Vd(itraconazole)/F, CLm/(F (.) fm) and F-rel, respectively. The correlation between random effects of CLitraconazole and Vd((itraconazole)) was 0.69. Conclusion: The developed population pharmacokinetic model adequately described the pharmacokinetics of itraconazole and its active metabolite, hydroxy-itraconazole, in paediatric patients with either cystic fibrosis or undergoing BMT. More appropriate dosing schedules have been developed for the oral solution and the capsules to secure a minimum therapeutic trough plasma concentration of 0.5 mg/L for these patients.