Pharmacokinetics and pharmacodynamics of orally administered clonidine: a model-based approach.

BACKGROUND/AIMS: The oral clonidine test is a diagnostic procedure performed in children with suspected growth hormone (GH) deficiency. It is associated with untoward effects, including bradycardia, hypotension and sedation. Serum clonidine levels have not previously been assessed during this test. METHODS: In 40 children referred for an oral clonidine test, blood samples were drawn for clonidine and GH. Vital statistics and sedation scores were recorded until 210 min post-dose. We explored the relationship between clonidine concentrations and effects such as GH peak and blood pressure. RESULTS: Of 40 participants, 5 children were GH deficient. Peak clonidine concentrations of 0.846 ± 0.288 ng/ml were reached after 1 h. Serum levels declined slowly, with concentrations of 0.701 ± 0.189 ng/ml 210 min post-dose. A large interindividual variation of serum levels was observed. During the procedure, systolic blood pressure dropped by 12.8%, diastolic blood pressure by 19.7% and heart rate by 8.4%. Moderate sedation levels were observed. Concentration-effect modeling showed that the amount of GH available for secretion as determined by previous bursts was an important factor influencing GH response. CONCLUSION: Clonidine concentrations during the test were higher than necessary according to model-based predictions. A lower clonidine dose may be sufficient and may produce fewer side effects.

Effect of altered AGP plasma binding on heart rate changes by S(-)-propranolol in rats using mechanism-based estimations of in vivo receptor affinity (K(B,vivo)).

In contrast to the impact of plasma protein binding on pharmacokinetics, no quantitative in vivo information is available on its impact on pharmacodynamics. The pharmacokinetic-pharmacodynamic relationship of the model drug S(-)-propranolol was evaluated using mechanism-based estimations of in vivo receptor affinity (K(B,vivo)), under conditions of altered plasma protein binding resulting from different levels of alpha-1-acid glycoprotein (AGP). Male Wistar Kyoto rats with isoprenaline-induced tachycardia received an intravenous infusion of S(-)-propranolol, on postsurgery day 2 (n = 7) and day 7 (n = 8) with elevated and normal plasma protein binding, respectively. Serial blood samples were taken in parallel to heart rate measurements. AGP concentrations at 2 and 7 days postsurgery were 708 +/- 274 and 176 +/- 111 microg/mL (mean +/- SE), respectively. Using nonlinear mixed effects modeling, AGP concentration was a covariate for intercompartmental clearance for the third compartment of the pharmacokinetic model of S(-)-propranolol. Individual values of AGP concentrations ranged between 110 and 1150 microg/mL, and were associated with K(B,vivo) values of S(-)-propranolol from 7.0 to 30 nM. Using the K(B,vivo) for S(-)-propranolol with correction for average values for normal and elevated plasma protein binding, nearly identical values were found. This confirms, strictly quantitative, earlier indications that plasma protein binding restricts the pharmacodynamics of S(-)-propranolol.

Influence of plasma protein binding on pharmacodynamics: Estimation of in vivo receptor affinities of beta blockers using a new mechanism-based PK-PD modelling approach.

The objective of this investigation was to examine in a systematic manner the influence of plasma protein binding on in vivo pharmacodynamics. Comparative pharmacokinetic-pharmacodynamic studies with four beta blockers were performed in conscious rats, using heart rate under isoprenaline-induced tachycardia as a pharmacodynamic endpoint. A recently proposed mechanism-based agonist-antagonist interaction model was used to obtain in vivo estimates of receptor affinities (K(B,vivo)). These values were compared with in vitro affinities (K(B,vitro)) on the basis of both total and free drug concentrations. For the total drug concentrations, the K(B,vivo) estimates were 26, 13, 6.5 and 0.89 nM for S(-)-atenolol, S(-)-propranolol, S(-)-metoprolol and timolol. The K(B,vivo) estimates on the basis of the free concentrations were 25, 2.0, 5.2 and 0.56 nM, respectively. The K(B,vivo)-K(B,vitro) correlation for total drug concentrations clearly deviated from the line of identity, especially for the most highly bound drug S(-)-propranolol (ratio K(B,vivo)/K(B,vitro) approximately 6.8). For the free drug, the correlation approximated the line of identity. Using this model, for beta-blockers the free plasma concentration appears to be the best predictor of in vivo pharmacodynamics.

Pharmacokinetic-pharmacodynamic modelling of fluvoxamine 5-HT transporter occupancy in rat frontal cortex.

BACKGROUND AND PURPOSE: The pharmacokinetic-pharmacodynamic (PK-PD) correlation of fluvoxamine 5-HT transporter (SERT) occupancy was determined in rat frontal cortex ex vivo. EXPERIMENTAL APPROACH: Rats (n=47) with permanent arterial and venous cannulas received a 30 min intravenous infusion of fluvoxamine (1 or 7.3 mg kg(-1)). At various time points after dosing, brains were collected for determination of fluvoxamine concentration and SERT occupancy. In addition, the time course of fluvoxamine concentration in plasma was determined up to the time of brain collection. In a separate study (n=26), the time course of fluvoxamine concentration in brain extracellular fluid (ECF) and plasma was determined. The results of the investigations were interpreted by nonlinear mixed effects modeling. KEY RESULTS: Highest SERT occupancy was reached at the first time point (10 or 15 min) and maintained for 1.5 and 7 h after 1 and 7.3 mg kg(-1), respectively. Thereafter, SERT occupancy decreased linearly at a rate of 8% h(-1). SERT occupancy could be directly related to plasma, brain ECF and brain tissue concentrations by a hyperbolic function (Bmax model). Maximal SERT occupancy (Bmax) was 95%. Estimated concentrations at half-maximal SERT occupancy (EC50) in plasma, ECF and brain tissue were 0.48, 0.22 and 14.8 ng mL(-1) respectively. The minimum value of the objective function decreased 12 points for ECF and brain tissue concentrations relative to plasma (P<0.01), presumably as a result of nonlinear brain distribution. CONCLUSIONS AND IMPLICATIONS: The proposed PK-PD model constitutes a useful basis for prediction of the time course of ex vivo SERT occupancy in behavioural studies with selective serotonin reuptake inhibitors.

Influence of biophase distribution and P-glycoprotein interaction on pharmacokinetic-pharmacodynamic modelling of the effects of morphine on the EEG.

BACKGROUND AND PURPOSE: The aim was to investigate the influence of biophase distribution including P-glycoprotein (Pgp) function on the pharmacokinetic-pharmacodynamic correlations of morphine's actions in rat brain. EXPERIMENTAL APPROACH: Male rats received a 10-min infusion of morphine as 4 mg kg(-1), combined with a continuous infusion of the Pgp inhibitor GF120918 or vehicle, 10 or 40 mg kg(-1). EEG signals were recorded continuously and blood samples were collected. KEY RESULTS: Profound hysteresis was observed between morphine blood concentrations and effects on the EEG. Only the termination of the EEG effect was influenced by GF120918. Biophase distribution was best described with an extended catenary biophase distribution model, with a sequential transfer and effect compartment. The rate constant for transport through the transfer compartment (k(1e)) was 0.038 min(-1), being unaffected by GF120918. In contrast, the rate constant for the loss from the effect compartment (k(eo)) decreased 60% after GF120918. The EEG effect was directly related to concentrations in the effect compartment using the sigmoidal E(max) model. The values of the pharmacodynamic parameters E(0), E(max), EC(50) and Hill factor were 45.0 microV, 44.5 microV, 451 ng ml(-1) and 2.3, respectively. CONCLUSIONS AND IMPLICATIONS: The effects of GF120918 on the distribution kinetics of morphine in the effect compartment were consistent with the distribution in brain extracellular fluid (ECF) as estimated by intracerebral microdialysis. However, the time-course of morphine concentrations at the site of action in the brain, as deduced from the biophase model, is distinctly different from the brain ECF concentrations.

Population pharmacokinetic modelling of non-linear brain distribution of morphine: influence of active saturable influx and P-glycoprotein mediated efflux.

BACKGROUND AND PURPOSE: Biophase equilibration must be considered to gain insight into the mechanisms underlying the pharmacokinetic-pharmacodynamic (PK-PD) correlations of opioids. The objective was to characterise in a quantitative manner the non-linear distribution kinetics of morphine in brain. EXPERIMENTAL APPROACH: Male rats received a 10-min infusion of 4 mg kg(-1) of morphine, combined with a continuous infusion of the P-glycoprotein (Pgp) inhibitor GF120918 or vehicle, or 40 mg kg(-1) morphine alone. Unbound extracellular fluid (ECF) concentrations obtained by intracerebral microdialysis and total blood concentrations were analysed using a population modelling approach. KEY RESULTS: Blood pharmacokinetics of morphine was best described with a three-compartment model and was not influenced by GF120918. Non-linear distribution kinetics in brain ECF was observed with increasing dose. A one compartment distribution model was developed, with separate expressions for passive diffusion, active saturable influx and active efflux by Pgp. The passive diffusion rate constant was 0.0014 min(-1). The active efflux rate constant decreased from 0.0195 min(-1) to 0.0113 min(-1) in the presence of GF120918. The active influx was insensitive to GF120918 and had a maximum transport (N(max)/V(ecf)) of 0.66 ng min(-1) ml(-1) and was saturated at low concentrations of morphine (C(50)=9.9 ng ml(-1)). CONCLUSIONS AND IMPLICATIONS: Brain distribution of morphine is determined by three factors: limited passive diffusion; active efflux, reduced by 42% by Pgp inhibition; low capacity active uptake. This implies blood concentration-dependency and sensitivity to drug-drug interactions. These factors should be taken into account in further investigations on PK-PD correlations of morphine.

Pharmacokinetic-pharmacodynamic modelling of S(-)-atenolol in rats: reduction of isoprenaline-induced tachycardia as a continuous pharmacodynamic endpoint.

BACKGROUND AND PURPOSE: For development of mechanism-based pharmacokinetic-pharmacodynamic (PK-PD) models, continuous recording of drug effects is essential. We therefore explored the use of isoprenaline in the continuous measurement of the cardiovascular effects of antagonists of beta-adrenoceptors (beta-blockers). The aim was to validate heart rate as a pharmacodynamic endpoint under continuous isoprenaline-induced tachycardia by means of PK-PD modelling of S(-)-atenolol. EXPERIMENTAL APPROACH: Groups of WKY rats received a 15 min i.v. infusion of 5 mg kg(-1) S(-)-atenolol, with or without i.v. infusion of 5 microg kg(-1) h(-1) isoprenaline. Heart rate was continuously monitored and blood samples were taken. KEY RESULTS: A three-compartment model best described the pharmacokinetics of S(-)-atenolol. The PK-PD relationship was described by a sigmoid Emax model and an effect compartment was used to resolve the observed hysteresis. In the group without isoprenaline, the variability in heart rate (30 b.p.m.) approximated the maximal effect (Emax=43+/-18 b.p.m.), leaving the parameter estimate of potency (EC50=28+/-27 ng ml(-1)) unreliable. Both precise and reliable parameter estimates were obtained during isoprenaline-induced tachycardia: 517+/-13 b.p.m. (E0), 168+/-15 b.p.m. (Emax), 49+/-14 ng ml(-1) (EC50), 0.042+/-0.012 min(-1) (k(eo)) and 0.95+/-0.34 (n). CONCLUSIONS AND IMPLICATIONS: Reduction of heart rate during isoprenaline-induced tachycardia is a reliable pharmacodynamic endpoint for beta-blockers in vivo in rats. Consequently this experimental approach will be used to investigate the relationship between drug characteristics and in vivo effects of different beta-blockers.

Dioxins, dioxin-like PCBs and non-dioxin-like PCBs in foodstuffs: occurrence and dietary intake in The Netherlands.

Data on occurrence of dioxins (polychlorinated dibenzo-p-dioxins [PCDDs] and dibenzofurans [PCDFs]), dioxin-like PCBs (polychlorinated non-ortho and mono-ortho biphenyls) and non-dioxin-like PCBs (as represented by the so-called indicator-PCBs: congeners 28, 52, 101, 118, 138, 153 and 180) in food products consumed in The Netherlands that were collected in measurement programs carried out during 1998 and 1999, and combined with food consumption data to assess the dietary intake of these persistent food contaminants. The estimated median life-long-averaged intake of the sum of dioxins and dioxin-like PCBs in the population is 1.2 pg WHO-TEQ (toxic equivalents) per kg body weight (bw) per day, while the estimated median life-long-averaged intake of indicator-PCBs is 5.6 ng per kg bw per day. The contribution of different food groups to the total intake of both dioxins + dioxin-like PCBs and non-dioxin-like PCBs is fairly uniformly distributed over the foods consumed: meat products (23% and 27%, respectively), dairy products (27% and 17%, respectively), fish (16% and 26%, respectively), eggs (4% and 5%, respectively), vegetable products (13% and 7%, respectively), and industrial oils and fats (17% and 18%, respectively). Compared with earlier intake estimations the present estimation shows a continued reduction in the intake of dioxins as well as PCBs. This reduction is related to the decrease in the concentration of these substances in the majority of foodstuffs. Nevertheless, a small part of the population still has a rather high life-long averaged intake: 8% of the population is exposed to intake levels above the tolerable weekly intake for dioxins and dioxin-like PCBs of 14 pg WHO-TEQ per kg bw per week, as recently derived by the Scientific Committee on Food of the European Commission. For the non-dioxin-like PCBs an internationally accepted maximum intake level is still lacking. However, to provide risk managers with a health-based guideline to prevent health effects of exposure to non-dioxin-like PCBs, the (international) derivation of a tolerable daily intake is recommended. Monitoring the dietary intake of PCBs is just as important as monitoring the intake of dioxins and dioxin-like PCBs, and attempts to decrease the exposure to both compound classes need continuous attention.

Modeling relationships between indoor and outdoor air quality.

Information about the ratio between indoor and outdoor concentrations (IO ratios) of air pollutants is a crucial component in human exposure assessment. The present study examines the relationship between indoor and outdoor concentrations as influenced by the combined effect of time patterns in outdoor concentrations, ventilation rate, and indoor emissions. Two different mathematical approaches are used to evaluate IO ratios. The first approach involves a dynamic mass balance model that calculates distributions of transient IO ratios. The second approach assumes a linear relationship between indoor and outdoor concentrations. We use ozone and benzene as examples in various modeling exercises. The modeled IO ratio distributions are compared with the results obtained from linear fits through plots of indoor versus outdoor concentrations.