Digoxin population pharmacokinetics from routine clinical data: role of patient characteristics for estimating dosing regimens.

Routine clinical pharmacokinetic data collected from patients receiving digoxin have been analysed to evaluate the role of patient characteristics for estimating dosing regimens. The data were analysed using NONMEM, a computer program designed for population pharmacokinetic analysis that allows pooling of data. The pharmacokinetic model of digoxin was described using a one-compartment steady-state model. The effect of a variety of developmental and demographic factors on clearance was investigated. NONMEM estimates indicate that digoxin clearance was influenced by the demographic variables of age, total body weight, serum creatinine and sex. The interindividual variability in digoxin clearance was modelled with additive error with an estimated standard deviation of 46.15 L day-1 and the intraindividual variability, or residual error was 0.209 ng mL-1. The dosing method based on clearance values obtained by NONMEM analysis allowed the prediction of the steady-state concentration as a function of maintenance dose with acceptable error for therapeutic drug monitoring.     

Population pharmacokinetics of free carbamazepine in adult Omani epileptic patients

OBJECTIVE: The aim of this study was to define the pharmacokinetic profile of free carbamazepine (F-CBZ) in adult Omani epileptic patients in order to improve on dosing schedules through population pharmacokinetic analysis using the NONMEM program. METHOD: Steady-state trough F-CBZ serum concentrations, carbamazepine (CBZ) dosing history and associated information were collected prospectively. RESULTS: Forty-eight patients with two or more available F-CBZ serum concentrations (total of 149 dose/serum concentration pairs) met our inclusion criteria. Patients were taking CBZ (200-1200 mg/day) in monotherapy. The analysis assumed a one-compartmental open model with first-order absorption and elimination. The apparent clearance (CL/F) and apparent volume of distribution (V/F) and their interindividual variabilities were estimated using the program. The population estimates for clearance (CL; modelled independently of dose) and volume of distribution were 13.2 +/- 0.6 l/h and 525 +/- 44 1, respectively. However, CL increased as a function of dosing rate and consequently was modelled as a linear function of steady-state concentration. In order to validate these results, the predictions of the population model were tested against data from 13 further patients subjected to the same inclusion criteria but who were not included in the original analysis. The predictions were good, being unbiased (P=0.31), and had an average deviation from the observed values of 18%. CONCLUSION: In order to establish steady-state dosage regimens, a population pharmacokinetic model is proposed, based on the patient's dose, to estimate the individual CL for an Omani epileptic patient receiving CBZ in monotherapy.     

Population pharmacokinetics of digoxin in pediatric patients

Digoxin pharmacokinetics were studied in a pediatric population with an age range of 6 days to 1 year using the population pharmacokinetic approach. Digoxin data were analyzed by mixed-effects modeling according to a one-compartment steady-state pharmacokinetic model using NONMEM software. The final model selected for the population prediction of digoxin clearance in pediatric patients was as follows: [equation: see text] Individual empirical Bayesian estimates were generated on the basis of the population estimates and were used to correlate the optimum dose of digoxin and patient age according to the following equation: [equation: see text] This equation and its derived nomogram may be used for the initial dosing of digoxin in children aged between 0 and 1 year. The use of this nomogram in routine monitoring requires further pharmacokinetic and clinical validation.     

Pharmacokinetics of dexamethasone in a rat model of rheumatoid arthritis

Dexamethasone (DEX) is often given for the treatment of rheumatoid arthritis and clinical dosing regimens of DEX have often been based empirically. This study tests whether the inflammation processes in a rat model of rheumatoid arthritis alters the clearance and volume of distribution of DEX when compared with healthy controls. Groups of healthy and arthritic male Lewis rats received either a low (0.225 mg/kg) or high (2.25 mg/kg) intramuscular dose of DEX. Arthritis was induced by intradermal injection of type II porcine collagen in incomplete Freund's adjuvant emulsion at the base of the tail. DEX was dosed in the arthritic animals 22 days post arthritis induction. Plasma DEX concentrations were determined by HPLC. Plasma concentration versus time data were analysed by non-compartmental analysis and pharmacokinetic model fitting using the population pharmacokinetic software NONMEM V. A linear bi-exponential pharmacokinetic model with extravascular input described the data for both healthy and arthritic animals. Clearance was the only parameter determined statistically different between both groups (healthy=1.05 l/h/kg, arthritic=1.19 l/h/kg). The steady-state volume of distribution for both groups was 4.85 l/kg. The slight difference in clearance was visibly undetectable and unlikely to produce meaningful changes in DEX disposition in arthritic rats. Copyright (c) 2008 John Wiley & Sons, Ltd.     

Pharmacokinetic variability of routinely administered bisoprolol in middle-aged and elderly Japanese patients

The nonlinear mixed effects model (NONMEM) was used to analyze the pharmacokinetics of routinely administered bisoprolol in middle-aged and elderly Japanese patients. The subjects consisted of 29 males and 11 females with a mean age of 63.5+/-10.1. Data on the plasma concentration of bisoprolol from 94 blood samples obtained at steady-state following repetitive oral administration were analyzed using the NONMEM program, where a one-compartment model with repetitive bolus dosing was parameterized in terms of oral clearance (CL/F) and apparent volume of distribution (V/F). Individual CL/F values were correlated with body weight (WT) and creatinine clearance (CLcr). The relation between CLcr and the CL/F of bisoprolol was not altered by the CYP2D6 and CYP2C19 genotypes, gender, or age. The mean CL/F value estimated with NONMEM was 0.0612.WT+1.15.CLcr (l/h), and the mean V/F value was 2.61.WT (l). The residual interindividual variability of CL/F and V/F were 22.0% and 12.6%, respectively. The pharmacokinetic variability of bisoprolol is small even in routinely treated Japanese patients, provided that both body weight and renal function are taken into account for the prediction of oral clearance of the drug.     

Dichloroacetate: population pharmacokinetics with a pharmacodynamic sequential link model

Dichloroacetate (DCA) is a small molecule that reduces ambient concentrations of lactate in man. It was the purpose of this study to develop pharmacokinetic and pharmacodynamic models for determination of a dose for a pivotal Phase III clinical trial of DCA in patients with traumatic brain injury (TBI). Population pharmacokinetic and pharmacodynamic models were developed for DCA using NONMEM software. The pharmacokinetic data were fit to a physiologic two-compartment model, and the pharmacodynamic data were fit to an indirect physiologic response model. Simulations were employed to evaluate various dosing strategies for consideration in a pivotal Phase III clinical trial of DCA. For the pharmacokinetic model, it was discovered that the clearance of DCA decreased on multiple dosing from 4.82 L/h to 1.07 L/h and that the pharmacokinetics and pharmacodynamics in TBI patients could not be predicted from normal volunteers. Population pharmacokinetic modeling and simulation of the expected effects of several dosing strategies were useful procedures for designing a Phase III trial.     

Detection of Carbamazepine-induced Changes in Valproic Acid Relative Clearance in Man by Simple Pharmacokinetic Screening

Selecting the optimum dose of valproic acid is difficult because the pharmacokinetics are complicated by inter-patient variability and by effects arising as a result of co-administration with other antiepileptic drugs. The multiple peak approach has been used to evaluate the effect of age, total body weight, dose, gender and co-medication (carbamazepine-induced change) on population estimates of valproic acid relative clearance. Routine clinical pharmacokinetic data (n = 479) were collected from 207 epilepsy patients on combination therapy. The data were analysed by a simple steady-state pharmacokinetic model with the use of NONMEM, a computer program designed for population pharmacokinetic analysis that enables pooling of data. NONMEM estimates suggested that the rate of valproic acid clearance in patients receiving concomitant administration of valproic acid and carbamazepine decreased non-linearly with increasing total body weight in the maturation process, and increased non-linearly with increasing valproic acid dose. The clearance in females was 5.7% less than in males. NONMEM estimates also suggested that the rate of valproic acid clearance increased non-linearly with increasing carbamazepine dose. Concomitant administration of valproic acid and carbamazepine with other antiepileptic drugs resulted in an increase in valproic acid clearance of 10%. The final regression model of valproic acid relative clearance was CL = 606TBW^0.168 ×DOSE^0.414 × CBZDOSE^0.095 × 0.943^GEN × 1.10^CO, where CL is the clearance (mL kg^? h^?), TBW is the total body weight (kg), DOSE is the dose of valproic acid, CBZDOSE is the dose of carbamazepine, GEN = 0 for males and 1 for females and CO = 0 for concomitant administration of valproic acid and carbamazepine and 1 for concomitant administration of valproic acid and carbamazepine with other antiepileptic drugs. This technique can be used to estimate the pharmacokinetic parameters of a population from sparse data collected during routine clinical care and to determine the extent to which patient characteristics influence drug pharmacokinetics.     

A pharmacokinetic model for alpha interferon administered subcutaneously

AIMS: To model the pharmacokinetic profiles of alpha interferon (alphaIFN) after a single subcutaneous (s.c.) injection of 3 million units of alpha 2b interferon, to correlate the pharmacokinetic parameters with patient demographic covariates, and to develop a limiting sampling strategy for determining the alphaIFN plasma area under the curve of concentration vs time (AUC). METHODS: The plasma alphaIFN pharmacokinetics were determined in 27 patients with chronic hepatitis C virus infection after the first s.c. injection of the drug. Ten patients had normal renal function and 17 were chronic haemodialysis patients. Plasma samples were assayed by an Elisa method. Concentration-time data was analysed by a population approach using NONMEM. RESULTS: The pharmacokinetic model which better described the concentration vs time data was a one-compartment model with two processes of absorption: a zero-order followed by a first-order process. The mean clearance of dialysis patients represented 37% (with 95% confidence interval: 30% -44%) of the mean value of the patients with normal renal function. The volume of distribution was significantly correlated to the body surface area. Bayesian analysis using NONMEM allowed determination of the individual plasma AUC from three samples within the 24 h period post s.c. injection. CONCLUSIONS: The present pharmacokinetic model will allow one to obtain individual parameters such as, the area under the curve of concentration vs time from a limited-sampling strategy, and to perform pharmacokinetic-pharmacodynamic analysis of combined alphaIFN plasma concentrations and viraemic data.     

Dosing of gentamicin in patients with end-stage renal disease receiving hemodialysis

The aim of this study was to evaluate dosing schedules of gentamicin in patients with end-stage renal disease and receiving hemodialysis. Forty-six patients were recruited who received gentamicin while on hemodialysis. Each patient provided approximately 4 blood samples at various times before and after dialysis for analysis of plasma gentamicin concentrations. A population pharmacokinetic model was constructed using NONMEM (version 5). The clearance of gentamicin during dialysis was 4.69 L/h and between dialysis was 0.453 L/h. The clearance between dialysis was best described by residual creatinine clearance (as calculated using the Cockcroft and Gault equation), which probably reflects both lean mass and residual clearance mechanisms. Simulation from the final population model showed that predialysis dosing has a higher probability of achieving target maximum concentration (Cmax) concentrations (> 8 mg/L) within acceptable exposure limits (area under the concentration-time curve [AUC] values > 70 and < 120 mg x h/L per 24 hours) than postdialysis dosing.     

The prediction of steady-state plasma phenobarbitone concentrations (following low-dose phenobarbitone) to refine its use as an indicator of compliance.

1. A model for predicting the steady-state plasma concentration of phenobarbitone following low-dose phenobarbitone used as an indicator of compliance was derived using data for 10 healthy volunteers. 2. Each volunteer was given a single 30 mg oral dose of phenobarbitone and the pharmacokinetics were described. Subsequently, volunteers were given phenobarbitone 2 mg daily for 28 days and a further pharmacokinetic profile determined during and after this period. 3. An initial predicted estimate of steady-state plasma drug concentration was made using each volunteer's demographic details. This estimate was revised by Bayesian analysis using single timed samples (24, 48, 72 or 96 h) following the single dose. 4. The model was tested on a further 10 healthy volunteers given a single 8 mg dose and who were subsequently given 2 mg daily for 28 days. 5. The revised estimate of peak steady-state plasma phenobarbitone concentration utilising the 96 h post-single dose concentration (356 ng ml-1) was least biased (mean prediction error +/- 95% CI = 10.6 +/- 19.8 ng ml-1) and most precise (root mean square error +/- 95% CI = 28.3 +/- 19.0 ng ml-1). In all cases the peak or trough steady-state drug concentration was within 13% of the predicted value. 6. The model reflected compliance accurately in a further eight volunteers with simulated partial (two-thirds) compliance. 7. The use of a predictive model using Bayesian analysis to estimate expected steady-state plasma phenobarbitone concentrations could increase further the usefulness of low-dose phenobarbitone as an indicator of compliance.     

Population-based investigation of relative clearance of digoxin in Japanese neonates and infants by multiple-trough screen analysis

OBJECTIVE: The steady-state concentrations of digoxin at trough levels were studied to establish the role of patient characteristics in estimating doses for digoxin using routine therapeutic drug monitoring data. METHOD: The data (n = 448) showing steady state after repetitive oral administration in 172 hospitalized neonates and infants were analyzed using Nonlinear Mixed Effect Model (NONMEM), a computer program designed to analyze pharmacokinetics in study populations by allowing pooling of data. Analysis of the pharmacokinetics of digoxin was accomplished using a simple steady-state pharmacokinetic model. The effects of a variety of developmental and demographic factors on the clearance of digoxin were investigated. RESULTS: Estimates generated using NONMEM indicated that clearance of digoxin (l.h-1) was influenced by the demographic variables of age, total body weight, serum creatinine, the coadministration of spironolactone, and the presence or absence of congestive heart failure. The interindividual variability in digoxin clearance was modeled with proportional errors with an estimated coefficient of variation of 32.1%, and the residual variability was 28.9%. In the validation set of 66 patients, the performance (bias, precision) of the final population model was good (mean prediction error -0.04 ng.ml-1; mean absolute prediction error 0.20 ng.ml-1).     

Population pharmacokinetic analysis of micafungin in Japanese patients with fungal infections

The object of this analysis was to develop a population pharmacokinetic model of micafungin, a new anti-fungal agent of the echinocandin class, to optimize dosing in Japanese patients with fungal infections. Population pharmacokinetics parameters were determined using NONMEM based on pharmacokinetic data from 198 subjects in seven clinical studies, comprising four phase I, two phase II and one pediatric phase III study. The healthy subjects received intravenous infusion of 2.5-150 mg micafungin. Adult and pediatric patients, age range of 8 month to 15 yeras old, were received 25-150 mg and 1-6 mg/kg daily, respectively. A total of 1825 micafungin plasma samples were available for this analysis. Two-compartment pharmacokinetic model was adopted. The clearance of micafungin was influenced by body weight in children and platelet counts (PLT). However the PLT accounted for less than 20% of the variation of micafungin clearance in Japanese subjects. In conclusions, body weight is the primary covariate factor in pediatric patients. The dose adjustment by body weight would be required only pediatric patients for the micafungin therapy in Japanese patients with fungal infection.     

Population pharmacokinetic modelling and estimation of dosing strategy for NXY-059, a nitrone being developed for stroke

BACKGROUND AND OBJECTIVES: NXY-059 (disufenton sodium, Cerovive, a nitrone with neuroprotective and free radical trapping properties (in experimental stroke) is under development for the treatment of acute stroke. The objectives of this study were to develop a population pharmacokinetic model for NXY-059 in acute stroke patients and to estimate individualised dosing strategies for NXY-059 using preclinical pharmacological and clinical pharmacokinetic information and knowledge of characteristics of the patient population. METHODS: NXY-059 was given as a continuous intravenous infusion for 72 hours, including a 1-hour loading infusion. Maintenance infusion rates were individualised based on creatinine clearance (CL(CR)). Population pharmacokinetic models were derived using NONMEM software. Optimal dosing strategies, individualised based on CL(CR) or bodyweight, were estimated using the population pharmacokinetic models, empirical covariate distributions relevant for the target population, and a target definition. Dosing strategies were selected based on target fulfillment criteria and parsimony. PATIENTS: Pharmacokinetic data from 179 patients with acute ischaemic or haemorrhagic stroke, included in two clinical studies, were used for the analyses. Patients were aged 34-92 years with varying degrees of renal impairment (estimated CL(CR) 20-143 mL/min). MAIN OUTCOME MEASURES AND RESULTS: The final population model based on data from both studies comprised a two-compartment model with unexplained interpatient variability for clearance (23% coefficient of variation [CV]) and central volume of distribution (40% CV). Part of the variability in clearance and volume of distribution was explained by CL(CR) and bodyweight, respectively. Typical clearance was estimated to 4.54 L/h in a patient with CL(CR) of 70 mL/min. The preferred dosing strategy for NXY-059 comprised an initial loading infusion (the same for all patients) followed by an individualised maintenance infusion on the basis of CL(CR) (three dosing categories) with cut-off values (at which infusion rates are incremented or decremented) of 50 and 80 mL/min. CONCLUSION: The results illustrate how an individualised dosing strategy, given a pharmacokinetic target, for NXY-059 was successfully optimised through estimation using the increasing pharmacokinetic and pharmacodynamic knowledge during a clinical drug development programme. The chosen dosing strategy of NXY-059 provides an easily adapted treatment regimen for acute stroke, resulting in early achievement of target plasma concentrations.     

Rationale for Fixed Dosing of Pertuzumab in Cancer Patients Based on Population Pharmacokinetic Analysis

Objective Our objectives were to develop the population pharmacokinetic (PK) for pertuzumab and examine the variability of steady-state trough serum concentrations (C_SS,trough) and exposure after fixed, body-weight-based, or body-surface area (BSA)-based dosing methods in cancer patients. Methods Pertuzumab was administered by IV infusion (every 3 weeks) either as a weight-based dose (0.5–15 mg/kg) or a fixed dose (420 or 1050 mg). Data from three clinical studies, comprising 153 patients and 1458 concentration-time points, were pooled for this analysis using NONMEM. Results A linear two-compartment model best described the data. Body weight and BSA were significant covariates affecting clearance (CL) and distribution volume (Vc), respectively. However, weight and BSA only explained small percentage of interpatient variability for CL and Vc, respectively. Simulation results indicated that PK profiles were very similar after the three dosing methods. Compared to fixed dosing, weight- and BSA-based dosing only reduced the population variability of C_SS,trough moderately. Conclusion A population PK model was developed for pertuzumab, the first monoclonal IgG1 antibody in a new class of agents known as HER dimerization inhibitors. In addition, our analyses demonstrate the feasibility of administering pertuzumab using a fixed dose in women with ovarian and breast cancers.     

Population pharmacokinetics of lamotrigine in elderly patients

Lamotrigine is being used more frequently in elderly patients. Dosing of lamotrigine in elderly patients is based largely on studies from younger adults and not evidence-based data from elderly patients. The goal of this study is to determine the pharmacokinetic parameters, such as clearance, and the factors that have a significant effect on these parameters to provide evidence-based information that can be used to dose elderly patients taking lamotrigine. Lamotrigine plasma concentrations from 148 elderly patients (aged 59-92 years) were used to develop a population pharmacokinetic model. Data were analyzed using NONMEM. Model evaluation was performed using the bootstrap approach and predictive check. The results showed that the blood urea nitrogen/serum creatinine ratio, weight, and phenytoin use significantly affect apparent clearance of lamotrigine. These results show that clinicians may need to take into account these covariates when dosing lamotrigine in this population.     

Population pharmacokinetic analysis of bisoprolol

The technique of population pharmacokinetic analysis was employed to study the variability in the dose concentration relationship of bisoprolol during its clinical development. The influence of demographic factors on the variability of clearance was investigated in 3 different populations: group I, patients (including an elderly group) with essential hypertension receiving multiple oral doses of bisoprolol 10 or 20mg for 3 months; group II, patients with different degrees of renal impairment and healthy controls; and group III, patients with different types of hepatic impairment and healthy controls. Patients and controls in groups II and III received only a single oral dose of bisoprolol 10mg. The 3 data sets were analysed separately, using a non-linear mixed effects model (the NONMEM program). A 2-compartment pharmacokinetic model with first-order absorption described the data adequately. The typical values of volume of central compartment, volume of distribution at steady-state and the absorption rate constant for the 3 populations were: for group I, 68L, 235L, and 0.7h-1; for group II, 28L, 179L, and 0.3h-1; and for group III, 55L, 256L, and 0.4h-1, respectively. Plasma clearance was related to age in group I, to serum creatinine in group II and to aspartate transaminase activity in group III. The 68% confidence limits for clearance and elimination half-life were 8.2 to 21.5 L/h and 7.6 to 19.7h, respectively, for 50-year-old patients in group I. The analysis predicted that progressive increases in serum creatinine or aspartate transaminase activity will result in only a 50% reduction of clearance.     

A three-step approach combining Bayesian regression and NONMEM population analysis: application to midazolam

NONMEM, the only available supported program for population pharmacokinetic analysis, does not provide the analyst with individual subject parameter estimates. As a result, the relationship between pharmacokinetic parameters and demographic factors such as age, gender, and body weight cannot be sought by plotting demographic factors vs. kinetic parameters. To overcome this problem, we devised a three-step approach. In step 1, an initial NONMEM analysis provides the population pharmacokinetic parameters without taking into account the demographic factors. Step 2 consists of individual bayesian regressions using the measured drug concentrations for each subject and the population pharmacokinetic parameters obtained in step 1. The bayesian parameter estimates of the individual subject can be plotted against the demographic factors of interest. From the scatter plots, it can be seen which are the demographic factors that appear to affect the pharmacokinetic parameters. In step 3, the NONMEM analysis is resumed, and the demographic factors found in step 2 are entered into the NONMEM regression model in a stepwise manner. This method was used to analyze the pharmacokinetics of midazolam in 64 subjects from 714 plasma concentrations and 11 demographic factors. CL (elimination clearance) and V1 were found to be a function of body weight. Age and liver disease were found to decrease CL. Of the 11 demographic factors recorded for each patient, none was found to influence VSS or intercompartmental clearance.     

The utility of mixed-effects covariate analysis in rapid selection of doses in pediatric subjects: a case study with fexofenadine hydrochloride

Fexofenadine hydrochloride is a non-sedating antihistamine that is used in the treatment of symptoms associated with seasonal allergic rhinitis and chronic idiopathic urticaria. A pooled analysis of pharmacokinetic data from children 6 months to 12 years of age and adults was conducted to identify the dose(s) in children that produce exposures comparable to those in adults for the treatment of seasonal allergic rhinitis. The pharmacokinetic parameter database included peak and overall exposure data from 269 treatment exposures from 136 adult subjects, and 90 treatment exposures from 77 pediatric allergic rhinitis patients. The data were pooled and analysed using NONMEM software, version 5.0. A covariate model based on body weight and age and a power function model based on body weight were identified as appropriate models to describe the variability in fexofenadine oral clearance and peak concentration, respectively. Individual oral clearance estimates were on average 44%, 36% and 61% lower in children 6 to 12 years (n=14), 2 to 5 years (n=21), and 6 months to 2 years (n=42), respectively, compared with adults. Trial simulations (n=100) were carried out based on the final pharmacostatistical models and parameter estimates to identify the appropriate dose(s) in children relative to the marketed dose of 60 mg fexofenadine hydrochloride in adults. The trials were designed as crossover studies in 18 subjects comprising various potential dosing regimens with and without weight stratification. Pharmacokinetic parameter variability was assumed to have a log-normal distribution. Individual weights and ages were simulated using mean (SD) estimates derived from the studies used in this analysis and proportional measurement/model mis-specification errors derived from the analysis were incorporated into the simulation. The results indicated that a 30 mg dose of fexofenadine hydrochloride administered to children 1 to 12 years of age and weighing >10.5 kg and a 15 mg dose administered to children 6 months and older and weighing 相似文献   

Postinduction carbamazepine clearance in an adult psychiatric population

The objective of this study was to describe the postinduction clearance of carbamazepine (CBZ) in adult psychiatric patients by population pharmacokinetic analysis using the NONMEM program. Specifically, an estimate of CBZ clearance and insight into the effect of common patient characteristics on clearance were sought. Steady-state trough CBZ serum concentrations, CBZ dosing history, concomitant drug administration, and other data from 45 patients were collected retrospectively. A one-compartment model with first-order absorption and first-order elimination was used, with absorption rate, bioavailability, and volume of distribution fixed to literature values. No evidence was found that race, sex, age, ethanol use, smoking, and concomitant lithium significantly affected CBZ clearance. In the final model, clearance was based on lean body weight. The coefficient of variation for clearance estimates was 16.5%. Residual variability was modest. Estimates for volume of distribution, rates of absorption and elimination, and bioavailability could not be pursued rigorously. Although these results may assist in understanding CBZ disposition in this population, their general clinical application should be approached with caution.     

Population pharmacokinetics of procainamide from routine clinical data

Routine clinical pharmacokinetic data collected from patients receiving procainamide were analysed to estimate population pharmacokinetic parameters. 116 plasma concentration determinations for procainamide and 14 timed urine collections for the drug and its major metabolite N-acetylprocainamide (NAPA) were obtained from 39 patients, mostly males. The data were analysed using NONMEM, a computer program designed for population pharmacokinetic analysis that allows pooling of data from many individuals. Estimates of the influence of weight, height, renal function, and the presence of congestive heart failure (CHF) on the renal clearance (CLR), acetylation clearance (CLA), miscellaneous metabolic clearance (CLO), and volume of distribution (Vd) of procainamide were obtained. The mean (SE) CLR, CLA, CLO and Vd for procainamide in a 70kg patient with normal renal function were estimated to be 14.4 (2.3) L/h, 10.1 (1.7) L/h, 1.2 (1.3) L/h, and 136.0 (20.0) L, respectively. These pharmacokinetic parameters vary linearly with bodyweight; height adds no information if weight is known. The presence of CHF has no significant effect on either CLO or Vd, but reduces CLA and CLR by 11% (p less than 0.01). Even after adjustments for CHF, renal function and weight, the total clearance and Vd of procainamide vary unpredictably among individuals, with a coefficient of variation between 30 and 40%, and less than 50%, respectively.