ABCB1 and cytochrome P450 genotypes and phenotypes: influence on methadone plasma levels and response to treatment

BACKGROUND AND OBJECTIVE: The in vivo implication of various cytochrome P450 (CYP) isoforms and of P-glycoprotein on methadone kinetics is unclear. We aimed to thoroughly examine the genetic factors influencing methadone kinetics and response to treatment. METHODS: Genotyping for CYP1A2, CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, ABCB1, and UGT2B7 polymorphisms was performed in 245 patients undergoing methadone maintenance treatment. To assess CYP3A activity, the patients were phenotyped with midazolam. RESULTS: The patients with lower CYP3A activity presented higher steady-state trough (R,S)-methadone plasma levels (4.3, 3.0, and 2.3 ng/mL x mg for low, medium, and high activity, respectively; P = .0002). As previously reported, CYP2B6*6/*6 carriers had significantly higher trough (S)-methadone plasma levels (P = .0001) and a trend toward higher (R)-methadone plasma levels (P = .07). CYP2D6 ultrarapid metabolizers presented lower trough (R,S)-methadone plasma levels compared with the extensive or intermediate metabolizers (2.4 and 3.3 ng/mL x mg, respectively; P = .04), whereas CYP2D6 poor metabolizer status showed no influence. ABCB1 3435TT carriers presented lower trough (R,S)-methadone plasma levels (2.7 and 3.4 ng/mL . mg for 3435TT and 3435CC carriers, respectively; P = .01). The CYP1A2, CYP2C9, CYP2C19, CYP3A5, and UGT2B7 genotypes did not influence methadone plasma levels. Only CYP2B6 displayed a stereoselectivity in its activity. CONCLUSION: In vivo, CYP3A4 and CYP2B6 are the major CYP isoforms involved in methadone metabolism, with CYP2D6 contributing to a minor extent. ABCB1 genetic polymorphisms also contribute slightly to the interindividual variability of methadone kinetics. The genetic polymorphisms of these 4 proteins had no influence on the response to treatment and only a small influence on the dose requirement of methadone.

Functional imaging of the parotid glands using blood oxygenation level dependent (BOLD)-MRI at 1.5T and 3T

PURPOSE: To evaluate the function of the parotid glands before and during gustatory stimulation, using an intrinsic susceptibility-weighted MRI method (blood oxygenation level dependent, BOLD-MRI) at 1.5T and 3T. MATERIALS AND METHODS: A total of 10 and 13 volunteers were investigated at 1.5T and 3T, respectively. Measurements were performed before and during gustatory stimulation using ascorbate. Circular regions of interest (ROIs) were delineated in the left and right parotid glands, and in the masseter muscle for comparison. The effects of stimulation were evaluated by calculating the difference between the relaxation rates, DeltaR(2)*. Baseline and stimulation were statistically compared (Student's t-tests), merging both parotid glands. RESULTS: The averaged DeltaR(2)* values prestimulation obtained in all parotid glands were stable (-0.61 to 0.38 x 10(-3) seconds(-1)). At 3T, these values were characterized by an initial drop (to -2.7 x 10(-3) seconds(-1)) followed by a progressive increase toward the baseline. No significant difference was observed between baseline and parotid gland stimulation at 1.5T, neither for the masseter muscle at both field strengths. A considerable interindividual variability (over 76%) was noticed at both magnetic fields. CONCLUSION: BOLD-MRI at 3T was able to detect DeltaR(2)* changes in the parotid glands during gustatory stimulation, consistent with an increase in oxygen consumption during saliva production.

Modeling of transfer kinetics at the serum-cerebrospinal fluid barrier in rabbits with experimental meningitis: application to grepafloxacin

The goals of the present study were to model the population kinetics of in vivo influx and efflux processes of grepafloxacin at the serum-cerebrospinal fluid (CSF) barrier and to propose a simulation-based approach to optimize the design of dose-finding trials in the meningitis rabbit model. Twenty-nine rabbits with pneumococcal meningitis receiving grepafloxacin at 15 mg/kg of body weight (intravenous administration at 0 h), 30 mg/kg (at 0 h), or 50 mg/kg twice (at 0 and 4 h) were studied. A three-compartment population pharmacokinetic model was fit to the data with the program NONMEM (Nonlinear Mixed Effects Modeling). Passive diffusion clearance (CL(diff)) and active efflux clearance (CL(active)) are transfer kinetic modeling parameters. Influx clearance is assumed to be equal to CL(diff), and efflux clearance is the sum of CL(diff), CL(active), and bulk flow clearance (CL(bulk)). The average influx clearance for the population was 0.0055 ml/min (interindividual variability, 17%). Passive diffusion clearance was greater in rabbits receiving grepafloxacin at 15 mg/kg than in those treated with higher doses (0.0088 versus 0.0034 ml/min). Assuming a CL(bulk) of 0.01 ml/min, CL(active) was estimated to be 0.017 ml/min (11%), and clearance by total efflux was estimated to be 0.032 ml/min. The population kinetic model allows not only to quantify in vivo efflux and influx mechanisms at the serum-CSF barrier but also to analyze the effects of different dose regimens on transfer kinetic parameters in the rabbit meningitis model. The modeling-based approach also provides a tool for the simulation and prediction of various outcomes in which researchers might be interested, which is of great potential in designing dose-finding trials.

Prognosis in pediatric hematologic malignancies is associated with serum concentration of mannose-binding lectin-associated serine protease-2 (MASP-2)

BACKGROUND: Mannose-binding lectin (MBL) and MBL-associated serine protease-2 (MASP-2) are key components of the lectin pathway of complement activation. Their serum concentrations show a wide interindividual variability. This study investigated whether the concentration of MBL and MASP-2 is associated with prognosis in pediatric patients with cancer. METHODS: In this retrospective multicenter study, MBL and MASP-2 were measured by commercially available ELISA in frozen remnants of serum taken at diagnosis. Associations of overall survival (OS) and event-free survival (EFS) with MBL and MASP-2 were assessed by multivariate Cox regression accounting for prognostically relevant clinical variables. RESULTS: In the 372 patients studied, median serum concentration of MBL was 2,808 microg/L (range, 2-10,060) and 391 microg/L (46-2,771) for MASP-2. The estimated 4-year EFS was 0.60 (OS, 0.78). In the entire, heterogeneous sample, MBL and MASP-2 were not significantly associated with OS or EFS. In patients with hematologic malignancies, however, higher MASP-2 was associated with better EFS in a significant and clinically relevant way (hazard ratio per tenfold increase (HR), 0.22; 95% CI, 0.09-0.54; P = 0.001). This was due to patients with lymphoma (HR, 0.11; 95% CI, 0.03-0.47; P = 0.003), but less for those with acute leukemia (HR, 0.35; 95% CI, 0.11-1.15; P = 0.083). CONCLUSION: In this study, higher MASP-2 was associated with better EFS in pediatric patients with hematologic malignancies, especially lymphoma. Whether MASP-2 is an independent prognostic factor affecting risk stratification and anticancer therapy needs to be assessed in prospective, disease-specific studies.

The pharmacokinetic profile of fesoterodine: similarities and differences to tolterodine

BACKGROUND: Fesoterodine is a new antimuscarinic agent developed for the treatment of overactive bladder. Fesoterodine itself is inactive and is rapidly and extensively converted by ubiquitous esterases to its principal active moiety, 5-hydroxymethyl tolterodine (5-HMT). 5-HMT is formed via biotransformation of both fesoterodine and tolterodine, albeit by different metabolising enzymes, viz. esterases and CYP2D6 respectively. Tolterodine is a potent muscarinic receptor antagonist and has been used for the treatment of overactive bladder for over ten years. The objective of this study was to establish the pharmacokinetic profile of fesoterodine and to highlight ist potential pharmacokinetic advantages over tolterodine. DESIGN: Single-centre, open-label, randomised, 4-way crossover study in a total of 24 healthy male volunteers. Single oral doses of 4, 8, or 12 mg fesoterodine were administered after an overnight fast. In addition, the 8 mg dose was also administered after a standard high-fat and high-calorie breakfast. Blood and urine samples for the analysis of 5-HMT were collected before and multiple times after drug administration for pharmacokinetic analysis. RESULTS: The mean peak plasma concentration (Cmax) of 5-HMT and the mean area under the time versus concentration curve (AUC) increased proportionally with the fesoterodine dose. These two parameters were some 2-fold higher in CYP2D6 poor metabolisers, whereas the time to peak plasma concentration (tmax) and half life (t1/2) were not influenced by the dose or the CYP2D6 metaboliser status. If fesoterodine was taken following a high-fat breakfast, we observed small increases in Cmax and AUC. In spite of these modest genetic influences and food effects on the pharmacokinetics of fesoterodine, the overall interindividual variability in Cmax levels was relatively little compared to previously published reports using tolterodine. CONCLUSIONS: Due to the esterase-mediated cytochrome P450-independent formation of 5-HMT and involvement of multiple metabolic and renal excretion pathways in the elimination of 5-HMT, the effects of patient-intrinsic and -extrinsic factors on the pharmacokinetics of fesoterodine are only modest, with some 2-fold higher 5-HMT exposure. Therefore, in contrast to tolterodine, no reduction of fesoterodine dosage is required under conditions of reduced elimination. In most cases of drug interaction or renal/hepatic impairment, the fesoterodine dose may be increased to 8 mg/day based on individual patients' response, or patients may be required to remain at the initial recommended dose of 4 mg/day.

Anesthesia or sedation for gastroenterologic endoscopies

PURPOSE OF REVIEW: Because propofol is the sedative preferred by gastroenterologists, we focus this review on gastroenterologist-directed propofol sedation, provide simulations of the respiratory depressant effect of different dosing protocols and give a perspective on future developments in computer-assisted sedation techniques. RECENT FINDINGS: Propofol use by nonanesthesiologists remains a contraindication in the package insert of propofol in most countries. Sedation guidelines produced by the American Society of Gastroenterology partially contradict those produced by the American Society of Anesthesiologists for sedation by nonanesthesiologists, whereas the German guidelines were developed with anesthesiologists involved. The use of fospropofol, recently approved by the US Food and Drug Administration for sedation, is considered an alternative to propofol by some gastroenterologists. Methodological errors in earlier pharmacological studies have to be solved before widespread use of fospropofol is justified, however. Our simulations show that dosing protocols with small boluses administered at reasonable intervals induce less respiratory depression than large boluses. Interindividual variability of propofol-induced respiratory depression is illustrated by different pharmacokinetic and dynamic parameter sets used in the simulation. Two computer-assisted propofol infusion systems are currently being investigated. They not only incorporate the target effect but also the side effects, which may limit respiratory depression. SUMMARY: Propofol use by gastroenterologists may be well tolerated if appropriate patient selection, staff training, monitoring and low-dose sedation protocols are applied.

In vivo analysis of efavirenz metabolism in individuals with impaired CYP2A6 function

INTRODUCTION: The antiretroviral drug efavirenz (EFV) is extensively metabolized into three primary metabolites: 8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV. There is a wide interindividual variability in EFV plasma exposure, explained to a great extent by cytochrome P450 2B6 (CYP2B6), the main isoenzyme responsible for EFV metabolism and involved in the major metabolic pathway (8-hydroxylation) and to a lesser extent in 7-hydroxylation. When CYP2B6 function is impaired, the relevance of CYP2A6, the main isoenzyme responsible for 7-hydroxylation may increase. We hypothesize that genetic variability in this gene may contribute to the particularly high, unexplained variability in EFV exposure in individuals with limited CYP2B6 function. METHODS: This study characterized CYP2A6 variation (14 alleles) in individuals (N=169) previously characterized for functional variants in CYP2B6 (18 alleles). Plasma concentrations of EFV and its primary metabolites (8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV) were measured in different genetic backgrounds in vivo. RESULTS: The accessory metabolic pathway CYP2A6 has a critical role in limiting drug accumulation in individuals characterized as CYP2B6 slow metabolizers. CONCLUSION: Dual CYP2B6 and CYP2A6 slow metabolism occurs at significant frequency in various human populations, leading to extremely high EFV exposure.

Pharmacogenetics-based population pharmacokinetic analysis of etravirine in HIV-1 infected individuals

Objectives: Etravirine (ETV) is metabolized by cytochrome P450 (CYP) 3A, 2C9, and 2C19. Metabolites are glucuronidated by uridine diphosphate glucuronosyltransferases (UGT). To identify the potential impact of genetic and non-genetic factors involved in ETV metabolism, we carried out a two-step pharmacogenetics-based population pharmacokinetic study in HIV-1 infected individuals. Materials and methods: The study population included 144 individuals contributing 289 ETV plasma concentrations and four individuals contributing 23 ETV plasma concentrations collected in a rich sampling design. Genetic variants [n=125 single-nucleotide polymorphisms (SNPs)] in 34 genes with a predicted role in ETV metabolism were selected. A first step population pharmacokinetic model included non-genetic and known genetic factors (seven SNPs in CYP2C, one SNP in CYP3A5) as covariates. Post-hoc individual ETV clearance (CL) was used in a second (discovery) step, in which the effect of the remaining 98 SNPs in CYP3A, P450 cytochrome oxidoreductase (POR), nuclear receptor genes, and UGTs was investigated. Results: A one-compartment model with zero-order absorption best characterized ETV pharmacokinetics. The average ETV CL was 41 (l/h) (CV 51.1%), the volume of distribution was 1325 l, and the mean absorption time was 1.2 h. The administration of darunavir/ritonavir or tenofovir was the only non-genetic covariate influencing ETV CL significantly, resulting in a 40% [95% confidence interval (CI): 13–69%] and a 42% (95% CI: 17–68%) increase in ETV CL, respectively. Carriers of rs4244285 (CYP2C19*2) had 23% (8–38%) lower ETV CL. Co-administered antiretroviral agents and genetic factors explained 16% of the variance in ETV concentrations. None of the SNPs in the discovery step influenced ETV CL. Conclusion: ETV concentrations are highly variable, and co-administered antiretroviral agents and genetic factors explained only a modest part of the interindividual variability in ETV elimination. Opposing effects of interacting drugs effectively abrogate genetic influences on ETV CL, and vice-versa.

Vitamin D time profile based on the contribution of non-genetic and genetic factors in HIV-infected individuals of European ancestry

BACKGROUND Vitamin D deficiency is prevalent in HIV-infected individuals and vitamin D supplementation is proposed according to standard care. This study aimed at characterizing the kinetics of 25(OH)D in a cohort of HIV-infected individuals of European ancestry to better define the influence of genetic and non-genetic factors on 25(OH)D levels. These data were used for the optimization of vitamin D supplementation in order to reach therapeutic targets. METHODS 1,397 25(OH)D plasma levels and relevant clinical information were collected in 664 participants during medical routine follow up visits. They were genotyped for 7 SNPs in 4 genes known to be associated with 25(OH)D levels. 25(OH)D concentrations were analyzed using a population pharmacokinetic approach. The percentage of individuals with 25(OH)D concentrations within the recommended range of 20-40ng/ml during 12 months of follow up and several dosage regimens were evaluated by simulation. RESULTS A one-compartment model with linear absorption and elimination was used to describe 25(OH)D pharmacokinetics, while integrating endogenous baseline plasma concentrations. Covariate analyses confirmed the effect of seasonality, body mass index, smoking habits, the analytical method, darunavir/r and the genetic variant in GC (rs2282679) on 25(OH)D concentrations. 11% of the interindividual variability in 25(OH)D levels was explained by seasonality and other non-genetic covariates and 1% by genetics. The optimal supplementation for severe vitamin D deficient patients was 300000 IU two times per year. CONCLUSIONS This analysis allowed identifying factors associated with 25(OH)D plasma levels in HIV-infected individuals. Improvement of dosage regimen and timing of vitamin D supplementation is proposed based on those results.

Predicting 5-fluorouracil toxicity: DPD genotype and 5,6-dihydrouracil:uracil ratio

AIM Decreased DPD activity is a major cause of 5-fluorouracil (5-FU) toxicity, but known reduced-function variants in the DPD gene (DPYD) explain only a part of DPD-related 5-FU toxicities. Here, we evaluated the baseline (pretherapeutic) plasma 5,6-dihydrouracil:uracil (UH2:U) ratio as a marker of DPD activity in the context of DPYD genotypes. MATERIALS & METHODS DPYD variants were genotyped and plasma U, UH2 and 5-FU concentrations were determined by liquid chromatography-tandem mass spectrometry in 320 healthy blood donors and 28 cancer patients receiving 5-FU-based chemotherapy. RESULTS Baseline UH2:U ratios were strongly correlated with generally low and highly variable U concentrations. Reduced-function DPYD variants were only weakly associated with lower baseline UH2:U ratios. However, the interindividual variability in the UH2:U ratio was reduced and a stronger correlation between ratios and 5-FU exposure was observed in cancer patients during 5-FU administration. CONCLUSION These results suggest that the baseline UH2:U plasma ratio in most individuals reflects the nonsaturated state of DPD and is not predictive of decreased DPD activity. It may, however, be highly predictive at increased substrate concentrations, as observed during 5-FU administration.

Exaggerated pulmonary hypertension and right ventricular dysfunction in high-altitude dwellers with patent foramen ovale.

BACKGROUND There is considerable interindividual variability in pulmonary artery pressure among high-altitude (HA) dwellers, but the underlying mechanism is not known. At low altitude, a patent foramen ovale (PFO) is present in about 25% of the general population. Its prevalence is increased in clinical conditions associated with pulmonary hypertension and arterial hypoxemia, and it is thought to aggravate these problems. METHODS We searched for a PFO (transesophageal echocardiography) in healthy HA dwellers (n = 22) and patients with chronic mountain sickness (n = 35) at 3,600 m above sea level and studied its effects (transthoracic echocardiography) on right ventricular (RV) function, pulmonary artery pressure, and vascular resistance at rest and during mild exercise (50 W), an intervention designed to further increase pulmonary artery pressure. RESULTS The prevalence of PFO (32%) was similar to that reported in low-altitude populations and was not different in participants with and without chronic mountain sickness. Its presence was associated with RV enlargement at rest and an exaggerated increase in right-ventricular-to-right-atrial pressure gradient (25 ± 7 mm Hg vs 15 ± 9 mm Hg, P < .001) and a blunted increase in fractional area change of the right ventricle (3% [-1%, 5%] vs 7% [3%, 16%], P = .008) during mild exercise. CONCLUSIONS These findings show, we believe for the first time, that although the prevalence of PFO is not increased in HA dwellers, its presence appears to facilitate pulmonary vasoconstriction and RV dysfunction during a mild physical effort frequently associated with daily activity. TRIAL REGISTRY ClinicalTrials.gov; No.: NCT01182792; URL: www.clinicaltrials.gov.

The genetic mediation of individual differences in sensitivity to pain and its inhibition

The underlying bases of the considerable interindividual variability in pain-related traits are starting to be revealed. Although the relative importance of genes versus experience in human pain perception remains unclear, rodent populations display large and heritable differences in both nociceptive and analgesic sensitivity. The identification and characterization of particularly divergent populations provides a powerful initial step in the genetic analysis of pain, because these models can be exploited to identify genes contributing to the behavior-level variability. Ultimately, DNA sequence differences representing the differential alleles at pain-relevant genes can be identified. Thus, by using a combination of “top-down” and “bottom-up” strategies, we are now able to genetically dissect even complex biological traits like pain. The present review summarizes the current progress toward these ends in both humans and rodents.

Avaliação de microtécnicas de extração para análise de lamotrigina em plasma de pacientes epilépticos por eletroforese capilar

A lamotrigina (LTG) é um fármaco pertencente à classe das feniltriazinas utilizado no tratamento de crises epilépticas generalizadas e focais e no tratamento adjunto da epilepsia refratária. Devido à alta variabilidade interindividual, às interações medicamentosas e aos efeitos adversos apresentados durante a administração da LTG, a monitorização terapêutica nos pacientes que fazem uso deste fármaco é necessária para ajuste de dose individual e evitar os efeitos adversos. Assim, o objetivo deste trabalho foi a avaliação de duas técnicas de microextração: a microextração em fase líquida com fibras ocas (HF-LPME) e a microextração líquido-líquido dispersiva (DLLME) para análise da lamotrigina em amostras de plasma de pacientes epilépticos. Primeiramente foram definidas as condições eletroforéticas: foi utilizado um capilar de sílica fundida de 75 ?m de diâmetro interno e 50 cm de comprimento efetivo. O eletrólito de corrida (BGE) foi composto por ácido 2-morfolinoetanosulfônico (MES), na concentração de 130 mmol L-1 e pH 5,0. As análises foram realizadas à temperatura de 20°C e tensão de 15 kV. A amostra foi injetada hidrodinamicamente (0,5 psi por 10 s) e a detecção foi feita em 214 nm. Nestas condições a LTG e o padrão interno (PI), lidocaína, puderam ser analisados em menos de 7 minutos. A HF-LPME foi avaliada no modo de 3 fases, usando 500 ?L de plasma e 3,5 mL de solução fosfato de sódio 50 mmol L-1 pH 9,0 como fase doadora. O solvente utilizado para impregnar a fibra foi o 1-octanol. Como fase aceptora foram utilizados 60 ?L de solução de ácido clorídrico pH 4,0. Para avaliação da DLLME, foi necessária uma etapa de pré-tratamento da amostra (500 ?L de plasma) com 1 mL de acetonitrila. Após isto, 1,3 mL do sobrenadante foram adicionados a 4 mL de solução fosfato de sódio 50 mmol L-1 pH 9,0 e 120 ?L de clorofórmio (solvente extrator) foram injetados nesta amostra aquosa e 165 ?L de fase sedimentada foram recuperados. As características de desempenho analítico para ambos os métodos foram avaliadas, sendo obtida linearidade na faixa de concentração plasmática de 1-20 ?g/mL e limite inferior de quantificação (LIQ) de 1 ?g mL-1. Os ensaios de precisão e exatidão apresentaram valores de acordo com os guias oficiais. Além disso, os métodos foram seletivos, não apresentaram efeito residual e as amostras foram estáveis. Os valores de recuperação foram de 54,3 e 23% para HF-LPME e DLLME, respectivamente. Os métodos validados foram aplicados com sucesso em amostras de plasma de pacientes epilépticos em tratamento com a LTG. Além disso, as duas técnicas foram comparadas e a HF-LPME apresentou vantagens em relação à DLLME, mostrando ser uma técnica promissora para análise de matrizes complexas, com reduzido consumo de solvente orgânico e possibilidade de automação.

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.