Pharmacokinetic-pharmacodynamic modeling of CNS drug effects: an overview

Impact of pharmacokinetic-pharmacodynamic modelling in early clinical drug development

Early clinical pharmacology studies in healthy subjects are often dissociated from patient studies. In this review we encourage the use of modelling and simulation techniques to generate valuable information at an early stage of clinical development. We illustrate these principles by presenting 5 different case studies from a spectrum of therapeutic drug classes. Their application leads to a better understanding of drug characteristics early on, thereby facilitating the design of dose-finding studies in the target patient population and saving resources.

Pharmacokinetics and pharmacodynamics of metamizol in co-administration with morphine under acute and chronic treatments in arthritic rats

Objective

To investigate the relationship between metamizol pharmacokinetics and the antinociceptive effect produced after subcutaneous administration of metamizol (177.8 mg/kg) alone or in combination with morphine (3.2 mg/kg), under acute and chronic treatments.

Methods

Antinociception was assessed using the pain-induced functional impairment model in rat (PIFIR). Serial blood samples were collected from the same animals to study the pharmacokinetics of metamizol.

Key findings

The co-administration of the drugs in single dose, confirmed the potentiation of their individual antinociceptive effects. When the drugs were administered alone following the chronic schedule, a pronounced tolerance development to their antinociceptive effects was found, whereas it was significantly attenuated when they were administered together. Metamizol pharmacokinetics was unaltered by the presence of morphine. Plasma concentrations of 4-methylaminoantipyrine, an active metabolite markedly decreased under chronic administration.

Conclusions

The mechanism involved in the potentiation of the antinociceptive effect produced by the combination, cannot be explained by the interaction of morphine on metamizol pharmacokinetics. Other pharmacokinetic interactions along with known pharmacodynamic interactions in which metamizol active metabolites contribute, should be considered. The frequency of administration enhances tolerance development and induces metamizol elimination process.

Kinetics of drug action in disease states: towards physiology-based pharmacodynamic (PBPD) models

Gerhard Levy started his investigations on the "Kinetics of Drug Action in Disease States" in the fall of 1980. The objective of his research was to study inter-individual variation in pharmacodynamics. To this end, theoretical concepts and experimental approaches were introduced, which enabled assessment of the changes in pharmacodynamics per se, while excluding or accounting for the cofounding effects of concomitant changes in pharmacokinetics. These concepts were applied in several studies. The results, which were published in 45 papers in the years 1984-1994, showed considerable variation in pharmacodynamics. These initial studies on kinetics of drug action in disease states triggered further experimental research on the relations between pharmacokinetics and pharmacodynamics. Together with the concepts in Levy's earlier publications "Kinetics of Pharmacologic Effects" (Clin Pharmacol Ther 7(3): 362-372, 1966) and "Kinetics of pharmacologic effects in man: the anticoagulant action of warfarin" (Clin Pharmacol Ther 10(1): 22-35, 1969), they form a significant impulse to the development of physiology-based pharmacodynamic (PBPD) modeling as novel discipline in the pharmaceutical sciences. This paper reviews Levy's research on the "Kinetics of Drug Action in Disease States". Next it addresses the significance of his research for the evolution of PBPD modeling as a scientific discipline. PBPD models contain specific expressions to characterize in a strictly quantitative manner processes on the causal path between exposure (in terms of concentration at the target site) and the drug effect (in terms of the change in biological function). Pertinent processes on the causal path are: (1) target site distribution, (2) target binding and activation and (3) transduction and homeostatic feedback.

Electroencephalogram effects of armodafinil: comparison with behavioral alertness

Development of central nervous system-acting drugs would be enhanced by suitable biomarkers that reflect the targeted pathophysiologic brain state. The electroencephalogram (EEG) has several characteristics of an ideal biomarker and can be promptly adapted to pre-clinical and clinical testing. The aim of this study was to evaluate EEG as a measure of the wakefulness-promoting effect of armodafinil in sleep deprived healthy subjects. Armodafinil pharmacodynamics were simultaneously assessed by EEG- and behavioral-based measures including a well-established measure of alertness. Using two quantitative EEG-based measures-power spectral and event-related brain activity analyses-we observed that armodafinil mitigated the slowing of brain activity and the decrease of the event-related brain activity caused by sleep deprivation. Armodafinil-induced changes in EEG are in agreement and explain up to 73.1% of the armodafinil-induced changes in alertness. Our findings suggest that EEG can serve as a marker of the wakefulness-promoting drug effect.

Characterization of bladder selectivity of antimuscarinic agents on the basis of in vivo drug-receptor binding

The in vivo muscarinic receptor binding of antimuscarinic agents (oxybutynin, solifenacin, tolterodine, and imidafenacin) used to treat urinary dysfunction in patients with overactive bladder is reviewed. Transdermal administration of oxybutynin in rats leads to significant binding of muscarinic receptors in the bladder without long-term binding in the submaxillary gland and the abolishment of salivation evoked by oral oxybutynin. Oral solifenacin shows significant and long-lasting binding to muscarinic receptors in mouse tissues expressing the M₃ subtype. Oral tolterodine binds more selectively to muscarinic receptors in the bladder than in the submaxillary gland in mice. The muscarinic receptor binding of oral imidafenacin in rats is more selective and longer-lasting in the bladder than in other tissues such as the submaxillary gland, heart, colon, lung, and brain, suggesting preferential muscarinic receptor binding in the bladder. In vivo quantitative autoradiography with (+)N-[¹¹C]methyl-3-piperidyl benzilate in rats shows significant occupancy of brain muscarinic receptors with the intravenous injection of oxybutynin, solifenacin, and tolterodine. The estimated in vivo selectivity in brain is significantly greater for solifenacin and tolterodine than for oxybutynin. Imidafenacin occupies few brain muscarinic receptors. Similar findings for oral oxybutynin were observed with positron emission tomography in conscious rhesus monkeys with a significant disturbance of short-term memory. The newer generation of antimuscarinic agents may be advantageous in terms of bladder selectivity after systemic administration.

Significance of protein binding in pharmacokinetics and pharmacodynamics

The significance of plasma protein binding on drug efficacy and, subsequently, the clinical relevance of changes in protein binding has been controversially discussed for decades. The uncertainty concerning the impact of plasma protein binding on a drug's pharmacological activity is, in part, related to the approach used when investigating and interpreting protein binding effects in vitro and in vivo. Frequently, a generalized one-size-fits-all approach, such as "protein binding does matter/does not matter," may not be applicable. An appropriate analysis requires careful consideration of both pharmacokinetic and pharmacodynamic processes, as they both contribute to the safety and efficacy of drugs. Therefore, the aim of this article is to provide a concise review of the theoretical concepts of protein binding, and to discuss relevant examples where applicable.

The forefront for novel therapeutic agents based on the pathophysiology of lower urinary tract dysfunction: bladder selectivity based on in vivo drug-receptor binding characteristics of antimuscarinic agents for treatment of overactive bladder

We have reviewed the binding of antimuscarinic agents, used to treat urinary dysfunction in patients with overactive bladder, to muscarinic receptors in target and non-target tissues in vivo. Transdermal administration of oxybutynin in rats led to significant binding in the bladder without long-term binding in the submaxillary gland and the abolishment of salivation evoked by oral oxybutynin. Oral solifenacin showed significant and long-lasting binding to muscarinic receptors in mouse tissues expressing the M(3) subtype. Oral tolterodine bound more selectively to muscarinic receptors in the bladder than in the submaxillary gland in mice. The muscarinic receptor binding activity of oral darifenacin in mice was shown to be pronounced and long-lasting in the bladder, submaxillary gland, and lung. In vivo quantitative autoradiography using (+)N-[(11)C]methyl-3-piperidyl benzilate in rats showed significant occupancy of brain muscarinic receptors on intravenous injection of oxybutynin, propiverine, solifenacin, and tolterodine. The estimated in vivo bladder selectivity compared to brain was significantly greater for solifenacin and tolterodine than oxybutynin. Darifenacin occupied few brain muscarinic receptors. Similar findings were also observed with positron emission tomography in conscious rhesus monkeys. The newer generation of antimuscarinic agents may be advantageous in the bladder selectivity after systemic administration.

Decreased Efficacy of GABAA-receptor Modulation by Midazolam in the Kainate Model of Temporal Lobe Epilepsy

PURPOSE

The objective of this investigation was to characterize quantitatively the time-dependent changes in midazolam (MDL) efficacy in the silent period after induction of status epilepticus (SE) in rats. The changes in MDL efficacy were correlated to changes in ex vivo GABA(A)-receptor expression.

METHODS

MDL efficacy was quantified by pharmacokinetic-pharmacodynamic (PK-PD) modeling by using the beta-frequency of the EEG as PD end point. Two PK-PD experiments were performed in each animal: the first experiment before and the second experiment at either day 4 or day 14 after SE. SE was induced by repetitive intraperitoneal injections with kainate. GABA(A)-receptor expression was determined by ex vivo autoradiography with [(3)H]flumazenil.

RESULTS

The concentration versus EEG effect relation of midazolam was successfully described by the sigmoidal E(max) model. The maximal effect on the beta-frequency of the EEG (E(max)) was reduced to 51.6 +/- 35.6% and 25.8 +/- 33.7% of the original value at 4 and 14 days after induction of SE. The ex vivo study with [(3)H]flumazenil showed that the observed reductions in E(max) were paralleled by a reduction in GABA(A)-receptor density.

CONCLUSIONS

The efficacy of MDL is decreased in the silent period after SE, which can be partly accounted for by a reduction in GABA(A)-receptor density.

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.

Pharmacokinetic/pharmacodynamic modelling of the anti-hyperalgesic and anti-nociceptive effect of adenosine A1 receptor partial agonists in neuropathic pain

The objective of this investigation was to characterise the pharmacokinetic-pharmacodynamic correlation of adenosine A1 receptor partial agonists in the chronic constriction injury model of neuropathic pain. Following intravenous administration of 8-methylamino-N6-cyclopentyl-adenosine (MCPA; 10 mg/kg) and 2'deoxyribose-N6-cyclopentyl-adenosine (2'dCPA; 20 mg/kg), the time course of the effect on the mechanical paw pressure threshold was determined in conjunction with plasma concentrations. Population pharmacokinetic/pharmacodynamic analysis was applied to derive individual concentration-effect relationships. A composite model consisting of an E(max) model for the anti-hyperalgesic effect in combination with a linear model for the anti-nociceptive effect accurately described the concentration-effect relationship. For both compounds, a full anti-hyperalgesic effect was observed. The values of the EC50 for the anti-hyperalgesic effect were (mean+/-S.D.): 3170+/-1460 and 2660+/-1200 ng/ml for MCPA and 2'dCPA versus 178+/-51 ng/ml for the reference full agonist 5'deoxyribose-N6-cyclopentyl-adenosine (5'dCPA). The values of the slope for the anti-nociceptive effect were 1.9+/-0.30 and 1.2+/-0.20 g.microl/ng, respectively, versus 55+/-8 g microl/ng for 5'dCPA. Adenosine A1 receptor partial agonists behave as full agonists with regard to the anti-hyperalgesic effect in neuropathic pain, but the anti-nociceptive effect is diminished.

Pharmacodynamic assessment of the benztropine analogues AHN-1055 and AHN-2005 using intracerebral microdialysis to evaluate brain dopamine levels and pharmacokinetic/pharmacodynamic modeling

PURPOSE

The benztropine (BZT) analogues bind with high affinity to the dopamine transporter (DAT) and demonstrate a behavioral and pharmacokinetic profile unlike that of cocaine. The development of a predictive pharmacokinetic/pharmacodynamic (PK/PD) model to characterize the concentration-effect relationship between the BZT analogues and brain dopamine (DA) levels is an important step in the evaluation of these compounds as potential cocaine abuse pharmacotherapies. Hence, the objective of this study was to mathematically characterize the PD of BZT analogues and cocaine, using appropriate PK/PD models.

METHODS

Dialysis probes were stereotaxically implanted into the nucleus accumbens of Sprague-Dawley rats (275-300 g). Extracellular fluid (ECF) DA levels were measured after intravenous administration of the BZT analogues AHN-1055 and AHN-2005, as well as cocaine using high performance liquid chromatography-electrochemical detection (HPLC-ECD). PD models were used to describe the relationship between the BZT analogues or cocaine and brain microdialysate DA, and suitability was based on standard goodness-of-fit criteria.

RESULTS

The BZT analogues produced a sustained increase in brain microdialysate DA levels in comparison to cocaine. The time of maximum concentration (T(max)) for brain microdialysate DA was 2 h for AHN-1055 and 1 h for AHN-2005 compared to a T(max) of 10 min for cocaine. The duration of brain microdialysate DA elevation was approximately 12-24 h for the BZTs in comparison to 1 h for cocaine. An indirect model with inhibition of loss of response and a sigmoid E(max) model best described the PK/PD for the BZT analogues and cocaine, respectively. The 50% of maximum inhibition (IC(50)) of the loss of DA was lower for AHN-2005 (226 +/- 27.5 ng/ml) compared to AHN-1055 (321 +/- 19.7 ng/ml). In addition, the EC(50) for cocaine was 215 +/- 11.2 ng/ml.

CONCLUSIONS

The slow onset and long duration of BZT analogue-induced DA elevation may avoid the reinforcing effects and craving of cocaine. Further, the developed models will be useful in characterizing the PK/PD of other analogues and aid in the assessment of the therapeutic efficacy of the BZT analogues as substitute medications for cocaine abuse.

Population pharmacokinetic–pharmacodynamic modelling of the anti-hyperalgesic effect of 5′deoxy-N6-cylopentyladenosine in the mononeuropathic rat

The objective of this investigation was to characterise the pharmacokinetic-pharmacodynamic correlation of 5'-deoxy-N6-cyclopentyl-adenosine (5'dCPA) in the chronic constriction injury model of neuropathic pain. Following intravenous administration of 5'dCPA (0.30 or 0.75 mg kg(-1)), the time course of the drug concentration in plasma was determined in conjunction with the effect on (1) the mechanical paw pressure and (2) the Von Frey Hair monofilament withdrawal threshold. Population pharmacokinetic-pharmacodynamic analysis was applied to derive individual concentration-effect relationships. For mechanical paw pressure a composite model consisting of an Emax model for the anti-hyperalgesic effect in combination with a linear model for the anti-nociceptive effect accurately described the data. The EC50 for the anti-hyperalgesic effect was 178+/-51 ng ml(-1) and the slope of the anti-nociceptive effect 0.055+/-0.008 g ml ng(-1). For the Von Frey Hair monofilament withdrawal threshold responders and non-responders were observed. Typically, in responders, full pain relief was observed at concentrations exceeding 100 ng ml(-1). The high plasma concentrations required for the anti-hyperalgesic effect relative to the receptor affinity are consistent with restricted transport of 5'dCPA to the site of action in the spinal cord and/or the brain.

Pharmacokinetic-pharmacodynamic modelling in the early development phase of anti-psychotics: a comparison of the effects of clozapine, S 16924 and S 18327 in the EEG model in rats

1. The use of pharmacokinetic/pharmacodynamic (PK/PD) analysis in early compound development was investigated in the rat for two developmental anti-psychotic compounds with clozapine as a positive control. 2. Three plasma samples were collected from each of eight animals according to a pre-defined sampling matrix allowing a total of 12 time points for PK analysis. Quantitative electroencephalography (QEEG), particularly the theta and beta frequencies, was used as a measurement of pharmacological effect. 3. PK/KD modelling of the sparse PK data available relative to a rich set of PD data was achieved using a population approach in NONMEM (IV). Individual PK parameter estimates were incorporated into a PK/PD model. 4. Qualitative EEG changes in rat and human were similar for clozapine, but different for the two developmental compounds, suggesting that changes in these PD parameters may not be specifically related to the anti-psychotic activity. 5. Although no definitive data are available concerning the signal specificity of EEG frequency bands with respect to dopaminergic or serotonergic receptor activity, qualitative and quantitative differences seen in EEG parameters are likely to result from the multiple receptor occupancy for these compounds. 6. The results confirm the value of population PK/PD modelling in conjunction with sparse sampling to enable determination of concentration effect relationships in the pre-clinical development programme of CNS-active drugs.

Clinical pharmacology of psychotropic drugs

No Abstract

A population pharmacokinetic-pharmacodynamic analysis of repeated measures time-to-event pharmacodynamic responses: the antiemetic effect of ondansetron

This paper presents and illustrates methodology for specifying, estimating, and evaluating a predictive model for repeated measures time-to-event responses. The illustrative example specifies a model of the antiemetic effect vs. concentration relationship for the 5-HT3 antagonist ondansetron in the human ipecac model for emesis. A key part of this model is a time-dependent log hazard function for emesis that is increased by ipecac administration and decreased by ondansetron concentration. The model is fit using an approximate maximum likelihood method. The data consist of the time free of emeses and, for those individuals with emetic episodes, the time(s) of the episode(s). Model evaluation is accomplished using residual plots adapted to time-to-event data and a "posterior predictive check" wherein observed data statistics are compared to those obtained from data simulated from the fitted model. The ondansetron concentration required to obtain a 50% reduction in the hazard of emesis is estimated to be 1.4 +/- 0.2 ng/ml, and the rate constant for elimination of ipecac-induced hazard is 1.5 +/- 0.2 hr-1.

Rate of change of blood concentrations is a major determinant of the pharmacodynamics of midazolam in rats

The objective of this investigation was to characterize quantitatively the influence of the rate of increase in blood concentrations on the pharmacodynamics of midazolam in rats. The pharmacodynamics of midazolam were quantified by an integrated pharmacokinetic-pharmacodynamic modelling approach. Using a computer controlled infusion technique, a linear increase in blood concentrations up to 80 ng ml(-1) was obtained over different time intervals of 16 h, resulting in rates of rise of the blood concentrations of respectively, 1.25, 1.00, 0.87, 0.46, 0.34 and 0.20 ng ml(-1) min(-1). In one group of rats the midazolam concentration was immediately brought to 80 ng ml(-1) and maintained at that level for 4 h. Immediately after the pretreatment an intravenous bolus dose was given to determine the time course of the EEG effect in conjunction with the decline of midazolam concentrations. The increase in beta activity (11.5-30 Hz) of the EEG was used as pharmacodynamic endpoint. For each individual animal the relationship between blood concentration and the EEG effect could be described by the sigmoidal Emax model. After placebo, the values of the pharmacodynamic parameter estimates were Emax = 82+/-5 microV, EC50,u = 6.4+/-0.8 ng ml(-1) and Hill factor = 1.4+/-0.1. A bell-shaped relationship between the rate of change of midazolam concentration and the value of EC50,u was observed with a maximum of 21+/-5.0 ng ml(-1) at a rate of change of 0.46 ng ml(-1) min(-1); lower values of EC50,u were observed at both higher and lower rates. The findings of this study show that the rate of change in plasma concentrations is an important determinant of the pharmacodynamics of midazolam in rats.

Phamacokinetic-pharmacodynamic modelling of behavioural responses

Drug concentrations at the site of action in studies on behavioural pharmacology, are seldom constant. Therefore, observed changes in behaviour can be due to the natural time course of behavioural processes, but equally to changes in drug concentration, and it is therefore crucial to separate the former from the latter. One solution is keeping drug concentrations constant. However, one can also exploit the variation in drug concentration caused by absorption, distribution and elimination of a drug. This is done by simultaneous measurement of drug effect and concentration, while the drug enters and leaves a biologically relevant compartment, such as blood or cerebrospinal fluid. The concept of determining concentration-effect curves in individual animals, by monitoring in parallel drug effect and changes in concentration in one single experiment, has not yet found wide application in behavioural studies. The fact that behavioural processes, like any other physiological process, change over time, may have contributed to the scarcity of pharmacokinetic-pharmacodynamic (PK/PD) studies in behavioural pharmacology. However, there are now mathematical techniques that allow PK/PD modelling even if the effect parameter changes over time or cannot be properly assessed in every instance. Here we use PK/PD modelling to characterize fear-induced ultrasonic vocalizations and the anxiolytic effect of buspirone. This approach reduces the number of animals required to assess concentration-effect relationships. More importantly, it allows the identification of differences in individual drug response over a wide range of concentrations. Consequently, we suggest that PK/PD modelling can be used as a tool to study drug-induced changes in behavioural response. An introduction in PK/PD modelling is presented.

Pharmacokinetic-pharmacodynamic modeling of the coexistence of stimulatory and sedative components for midazolam

Midazolam increased the shorter-response rate and decreased the reinforcement rate of a contingency-controlled timing behavior--a differential-reinforcement-of-low-rate 45-s schedule. The responding rate changes observed were immediately interpretable as functions of midazolam concentration during a 3-h session--a period for investigating the onset, peak, and disappearance of midazolam effect--in rats. That the midazolam pharmacokinetic-pharmacodynamic model was a direct application of our alprazolam pharmacokinetic-pharmacodynamic model implies that both drugs exhibit similar pharmacological effects. The two peaks of the shorter-response rate increases produced by midazolam were modeled as a stimulation-sedation model that consisted of two opposing effect-link sigmoidal Emax functions. The stimulation-sedation model suggested that midazolam possesses both stimulatory and sedative effects in a continuous but sequential fashion, and hypothesizes the coexistence of stimulation and sedation components for midazolam; this model may help delineate possible mechanisms for rebound side effects and of tolerance in humans. The reinforcement rate was, then, an index for evaluating the deficit in timing performance.

Influence of different fat emulsion-based intravenous formulations on the pharmacokinetics and pharmacodynamics of propofol

PURPOSE

The influence of different intravenous formulations on the pharmacokinetics and pharmacodynamics of propofol was investigated using the effect on the EEG (11.5-30 Hz) as pharmacodynamic endpoint.

METHODS

Propofol was administered as an intravenous bolus infusion (30 mg/kg in 5 min) or as a continuous infusion (150 mg/kg in 5 hours) in chronically instrumented male rats. Propofol was formulated as a 1% emulsion in an Intralipid 10%-like fat emulsion (Diprivan-10, D) or as a 1%- or 6% emulsion in Lipofundin MCT/LCT-10% (P1% and P6%, respectively). EEG was recorded continuously and arterial blood samples were collected serially for the determination of propofol concentrations using HPLC.

RESULTS

Following bolus infusion, the pharmacokinetics of the various propofol emulsions could adequately be described by a two-compartmental pharmacokinetic model. The average values for clearance (Cl), volume of distribution at steady-state (Vd,ss) and terminal half-life (t1/2, lambda 2) were 107 +/- 4 ml/min/kg, 1.38 +/- 0.06 l/kg and 16 +/- 1 min, respectively (mean +/- S.E. n = 22). No significant differences were observed between the three propofol formulations. After continuous infusion these values were 112 +/- 11 ml/min/kg, 5.19 +/- 0.41 l/kg and 45 +/- 3 min, respectively (mean +/- S.E., n = 20) with again no statistically significant differences between the three propofol formulations. Comparison between the bolus- and the continuous infusion revealed a statistically significant difference for both Vd,ss and t1/2, lambda 2 (p < 0.05), whereas Cl remained unchanged. In all treatment groups infusion of propofol resulted in a burst-suppression type of EEG. A profound hysteresis loop was observed between blood concentrations and EEG effect for all formulations. The hysteresis was minimized by a semi-parametric method and resulted in a biphasic concentration-effect relationship of propofol that was described non-parametrically. For P6% a larger rate constant onset of drug effect (t1/2,keo) was observed compared to the other propofol formulations (p < 0.05).

CONCLUSIONS

The pharmacokinetics and pharmacodynamics of propofol are not affected by to a large extent the type of emulsion nor by the concentration of propofol in the intravenous formulation.

Pharmacodynamics of temazepam in primary insomnia: assessment of the value of quantitative electroencephalography and saccadic eye movements in predicting improvement of sleep

BACKGROUND AND OBJECTIVE

Quantitative electroencephalographic parameters and saccadic eye movements are frequently used as pharmacodynamic measures of benzodiazepine effect. We investigated the relationship between these measures and the hypnotic effect.

METHODS

The correlation between the pharmacodynamic measures and sleep quality was determined in 21 patients with primary insomnia. The pharmacokinetic-pharmacodynamic relationships were characterized after administration of 20 mg oral temazepam. The hypnotic effect was determined on the basis of polysomnographic sleep recordings and a subjective sleep evaluation questionnaire. Correlations between pharmacodynamic measures and the improvement of sleep were investigated.

RESULTS

The pharmacokinetic-pharmacodynamic relationships for the parameters derived from electroencephalography and saccadic eye movements showed considerable interindividual variability. Administration of temazepam led to a significant improvement in the objective parameters sleep period efficiency, wake time after sleep onset, and sleep efficiency and in the subjective assessment of sleep quality. No significant correlations were observed between the pharmacokinetic-pharmacodynamic-derived parameters and the improvement in objective or subjective sleep parameters.

CONCLUSION

In subjects with primary insomnia the administration of 20 mg oral temazepam results in changes in both the pharmacodynamic measures and in quality of sleep. No individual correlations between the pharmacodynamic measures and quality of sleep were observed. We concluded that the investigated pharmacodynamic measures are of value in the first assessment of clinical efficacy and for the selection of the dose(s) to be investigated in subsequent trials that aim at showing clinical efficacy. However, the conclusive quantification of clinical efficacy should be performed only on the basis of the clinical end point itself.

Pharmacokinetic-pharmacodynamic modeling of mivacurium in rats

The pharmacokinetic (PK) and pharmacodynamic (PD) characteristics of the neuromuscular blocking agent mivacurium were evaluated separately in two groups of rats receiving 0.6 mg kg-1 of mivacurium in a 2.5-min intravenous continuous (iv) infusion. The PK parameters for mivacurium were determined in the first group. A two-compartment model describes the kinetics of mivacurium in plasma. The estimates of the apparent volume of distribution at steady-state and plasma clearance [mean(SE)] were 650 (123) mL kg-1 and 9.9 (0.75) mL min-1 kg-1, respectively. In the second group, the evoked tibialis anterior muscle tension was monitored. The PK parameters derived from the first group were used to compute mivacurium plasma concentrations (C) at the times the PD measurements were recorded in the second group. The concentration-neuromuscular effect [% depression of initial twitch tension (E)] relationship was analyzed by two approaches. (1) The relationship of estimated effect site concentrations versus E; a sigmoidal Emax model described the effect compartment concentrations versus E relationship. The estimate [mean(SE)] of Cess50 (steady-state plasma concentration eliciting half of maximum E) was 0.65 (0.01) microgram mL-1. The value [mean-(SE)] of Keo (rate constant of equilibration between plasma and effect site) was estimated at 0.32 (0.03) min-1. (2) The relationship of descending limb C versus E; a sigmoidal Emax model described such relationship. The estimate [mean(SE)] of C50 (post-infusion C eliciting half of maximum E) was 0.57(0.03) microgram mL-1. The PD properties of mivacurium were also evaluated in another two groups of animals receiving either 5- or 10-min continuous iv infusion; PK and PD parameters obtained from the 2.5-min infusion experiments were used to predict the time course of E in the groups receiving 0.6 mg kg-1 of mivacurium in 5- and 10-min infusions; simulations using the estimated parameters adequately describe the time course of E in those groups. The effect of mivacurium on the mean arterial blood pressure (MAP) was also investigated; a 10% nonsignificant decrease (p > 0.05) in MAP was found in all groups.

Comparison of methods for the assessment of central nervous system stimulant response after dextroamphetamine administration to healthy male volunteers

The objective of this investigation was to evaluate a series of potential pharmacodynamic measures of central nervous system stimulation, including quantitative electroencephalography (EEG) and neuroendocrine, mood, and psychomotor performance measures. The reproducibility and sensitivity of the measures were compared. The study was conducted in two parts. The first part investigated the interindividual and intraindividual variability associated with a series of potential pharmacologic response measures under baseline (i.e., drug-free) conditions. It was an open-label, three-period pilot study in which healthy male volunteers underwent a series of tests (EEG, a visual continuous performance task, a finger tapping task, and self-rated mood scales) repeatedly during each study period. The second part evaluated the sensitivity of a series of potential response measures to detect the effects of dextroamphetamine, and was a double-blind, placebo-controlled, four-period crossover study in nine healthy male volunteers. Subjects received 5 mg, 10 mg, or 20 mg of dextroamphetamine or placebo orally and underwent the same series of tests as in Part I in addition to blood collection for determination of serum prolactin and dextroamphetamine levels. Peripheral response to dextroamphetamine was assessed by heart rate and blood pressure measurement. The greatest variability among days, within days, and among participants was associated with the quantitative electroencephalographic parameters studied. First-session effects were apparent for several of the tests, including EEG. Consistent response on EEG (increased alpha power) to dextroamphetamine was observed only in the three subjects who had a baseline alpha activity greater than 35%. Serum prolactin levels were inversely associated with the amount of dextroamphetamine administered, with the largest decrease in serum prolactin levels observed after the 5-mg dose, and this finding was statistically significant. Mood scales showed that three of nine participants experienced dysphoria after at least one dose level of dextroamphetamine. The effect on mood was generally greater as the dose increased. Doses could not be distinguished based on the results of the performance tests. Serum prolactin concentration was the most sensitive measure of central nervous system stimulation on EEG produced by dextroamphetamine under these study conditions. Cardiovascular measures were more sensitive measures of dextroamphetamine effects than the central nervous system measures.

Acute tolerance to subjective but not cardiovascular effects of d-amphetamine in normal, healthy men

This is a descriptive report on the relationship between the pharmacokinetics and pharmacodynamics of d-amphetamine in healthy, normal volunteers. Six men, aged 22 to 31, attended two experimental sessions during which they received single oral doses of 20 mg of d-amphetamine. Plasma levels of drug and measures of drug effect were collected predrug and at regular intervals for 24 hours after drug administration. Plasma drug levels peaked at 4 hours and remained at detectable levels for 24 hours after drug administration. Subjective ratings, including "feel drug" and "feel high" peaked at 1 1/2 to 2 hours and returned to baseline levels by 3 to 4 hours. Evaluation of phase plots (i.e., drug effect vs. drug concentration) indicated that acute tolerance developed to the subjective but not to the cardiopressor effects of d-amphetamine. This finding implies that individuals who repeatedly administer the drug to maintain certain levels of subjective effects may increase plasma drug levels and physiologic effects to toxic levels.

Deoxyribose analogues of N6-cyclopentyladenosine (CPA): partial agonists at the adenosine A1 receptor in vivo

1. The purpose of the present study was to quantify the cardiovascular effects of the 2'-, 3'-, 5'-deoxyribose analogues of the selective adenosine A1 receptor agonist, N6-cyclopentyladenosine (CPA) in vivo. The blood concentration-effect relationships of the compounds were assessed in individual rats and correlated to their receptor binding characteristics. 2. The pharmacokinetics and pharmacodynamics of the compounds were determined after a single intravenous infusion of 0.80 mg kg (-1) (63 micromol kg(-1)of 2' dCPA. The heart rate (HR) and mean arterial blood pressure (MAP) were monitored continuously during the experiment and serial arterial blood samples were taken for analysis of drug concentration. 3. The relationship between blood concentrations and the reductions in both heart rate and blood pressure were described according to the sigmoidal Emax model. For the bradycardiac effect, the potencies based on free drug concentrations (EC50,u) of 5'dCPA, 3'dCPA, 2'dCPAin blood were 5.9 +/- 1.7, 18 +/- 4 and 260 +/- 70 ng ml (-1) (19 +/- 6, 56 +/- 11 and 830 +/- 210 nM), respectively, and correlated well with the adenosine A1 receptor affinity in vitro. The Emax value of 2'dCPA was significantly less than those of the other compounds, suggesting that this compound may be regarded as a partial agonist when compared to the other analogues. The rank order of the maximal reduction in heart rate of the compounds corresponded well with the order of the GTP-shifts, as determined in vitro. 4. It is concluded that deoxyribose derivatives of CPA may be partial agonists for the adenosine A1 receptor and may serve as tools for further investigation of adenosine receptor partial agonism in vivo.

Pharmacokinetics and pharmacodynamics of Ro 41-3696, a novel nonbenzodiazepine hypnotic

This report describes the first evaluation in humans of Ro 41-3696. Based on its preclinical profile, Ro 41-3696, a nonbenzodiazepine partial agonist at the benzodiazepine receptor, offers promising perspectives as an innovative hypnotic drug in that it does not exhibit most of the disadvantages associated with full agonists. Single oral doses of 0.1, 0.3, 1.0, 3.0, 10, and 30 mg were administered sequentially to six groups of six healthy male volunteers in a placebo-controlled, double-blind design. Tolerability was assessed and pharmacokinetic and pharmacodynamic measurements were conducted during a period of 28 hours after drug intake. Ro 41-3696 was well tolerated at all doses, causing no clinically relevant changes in vital signs or laboratory parameters. At doses of 10 and 30 mg there were signs of unsteady gait, indicating a central nervous system depressant effect. Pharmacokinetic analyses revealed that Ro 41-3696 was absorbed and eliminated rapidly (tmax = approximately 1 hour; t1/2 = approximately 4 hours). At all times plasma levels of Ro 41-3290, the desethylated derivative of Ro 41-3696, were higher than those of the parent drug (tmax and t1/2 values = approximately 2 and 8 hours, respectively). Area under the curve (AUC) data indicated dose-proportional pharmacokinetics for both Ro 41-3696 and Ro 41-3290. Performance in both a tracking and a memory search test was significantly affected by doses of 10 and 30 mg, and long-term memory, as assessed by a word learning and recall test, was slightly impaired at these doses. The results of this study support the initiation of therapeutic efficacy studies with Ro 41-3696 in doses up to approximately 5 mg and further exploration of the characteristics of Ro 41-3290.

Modelling of the pharmacodynamic interaction of an A1 adenosine receptor agonist and antagonist in vivo: N6-cyclopentyladenosine and 8-cyclopentyltheophylline

1. The purpose of this investigation was to develop a pharmacokinetic-pharmacodynamic model for the interaction between an adenosine A1 receptor agonist and antagonist in vivo. The adenosine A1 receptor agonist, N6-cyclopentyladenosine (CPA) and the antagonist, 8-cyclopentyltheophylline (CPT) were used as model drugs. The CPA-induced reduction in mean arterial pressure and heart rate were used as measurements of effect. 2. Four groups of eight rats each received 200 micrograms kg-1 of CPA i.v. in 5 min during a steady-state infusion of CPT at a rate of 0, 57, 114 or 228 micrograms kg-1 h-1. The haemodynamic parameters were continuously measured and frequent blood samples were taken to determine the pharmacokinetics of the drugs. 3. CPT had no influence on the pharmacokinetics of CPA and the baseline values of the haemodynamic variables. Furthermore, no clear antagonism by CPT was observed of the CPA-induced reduction in mean arterial pressure. However, CPT antagonized the effect on heart rate, and with increasing CPT concentrations, a parallel shift of the CPA concentration-effect relationship to the right was observed. 4. An agonist-antagonist interaction model was used to characterize the interaction quantitatively. On the basis of this model, the pharmacodynamic parameters of both CPA and CPT could be estimated. For CPA the values were (mean +/- s.e.): Emax = 198 +/- 11 b.p.m., EC50 = 2.1 +/- 0.7 ng ml-1, Hill factor = 2.3 +/- 0.6 and for CPT: EC50 = 3.7 +/- 0.3 ng ml-1 and Hill factor = 3.1 +/- 0.1. 5. It is concluded that the competitive agonist-antagonist interaction model may be of value to characterize quantitatively the pharmacodynamic interactions between adenosine A1 receptor ligands in vivo.

Pharmacotherapy of insomnia: practice and prospects

Insomnia is a complex complaint which is often multifactorial in origin. Pharmacotherapy can only be an adjunct in the treatment of insomnia and hypnotics should be given on an intermittent basis for short periods of time. An overview is presented of the currently available hypnotics, of which benzodiazepines are still the most widely prescribed. New drugs which bind to specific receptor subtypes or which are partial benzodiazepine receptor agonists might overcome the disadvantages associated with chronic benzodiazepine use, but more long-term investigations are needed.

Integrated pharmacokinetics and pharmacodynamics of the novel catechol-O-methyltransferase inhibitor tolcapone during first administration to humans

OBJECTIVES

To assess the tolerability, pharmacokinetics and pharmacodynamics of single oral doses of the novel catechol-O-methyltransferase (COMT) inhibitor tolcapone in healthy volunteers.

METHODS

In this double-blind, placebo-controlled, ascending-single-dose study, doses of 5 to 800 mg tolcapone were administered orally to eight sequential groups of six young healthy male volunteers. Adverse events, vital signs, and clinical laboratory variables were recorded. Pharmacokinetic parameters of tolcapone and its 3-O-methylmetabolite were determined. Pharmacodynamics were assessed by determination of COMT activity in erythrocytes.

RESULTS

Tolcapone was well tolerated at all dose levels and did not exert a detectable influence on vital sign measurements. The drug was rapidly absorbed and showed dose-proportional pharmacokinetics. Its mean elimination half-life was 2.0 +/- 0.8 hours (n = 42). Plasma levels of the 3-O-methylmetabolite of tolcapone were not proportional to dose, and its formation was delayed at higher doses. Its elimination half-life was 32 +/- 7 hours (n = 29). Tolcapone caused a rapid and reversible inhibition of COMT activity in erythrocytes. At doses of 200 mg and higher, COMT activity was inhibited by more than 80%. The pharmacokinetic-pharmacodynamic relationship could be described by an inhibitory Emax model and suggested that metabolites of tolcapone did not substantially contribute to its inhibitory activity.

CONCLUSIONS

The novel COMT inhibitor tolcapone was well tolerated at oral doses of 5 to 800 mg. Tolcapone concentration-dependently inhibited COMT activity in erythrocytes and exhibited dose-proportional kinetics. Further investigations into its applicability in the treatment of Parkinson's disease are warranted.

Pharmacokinetic modelling of the haemodynamic effects of the A2a adenosine receptor agonist CGS 21680C in conscious normotensive rats

1. The aim of the present investigation was to determine the relationship between the blood concentration and haemodynamic effects of the adenosine A2a receptor agonist, CGS 21680C (the sodium salt of 2-p-(2-carboxyethyl)phenylethylamino-5'-N-ethylcarboxamidoadeno sin e) in conscious normotensive rats. 2. Chronically cannulated rats were randomly assigned to three groups which received 300, 1000 or 3000 micrograms kg-1 (0.56, 1.9 or 5.6 mumol kg-1) of CGS 21680C intravenously over 15 min. The mean arterial blood pressure (MAP) and heart rate (HR) were monitored continuously during the experiment and serial arterial blood samples were taken for analysis of drug concentration. The ratio MAP/HR was also calculated, which may reflect changes in total peripheral resistance on the assumption that no changes in stroke volume occur. 3. For each individual rat the reduction in mean arterial pressure was related to the blood concentration according to the sigmoidal Emax model. The concentration-effect relationships were consistent for the different treatment groups. The potency based on free drug concentrations (EC50,u) was 5.8 ng ml-1 (11 nM) (mean +/- s.e.; n = 19) and correlated well with the reported adenosine A2a receptor affinity (Ki 19 nM). In comparison with the reduction in blood pressure, CGS 21680C exhibited a greater potency for the reduction of the ratio MAP/HR. 4. It is concluded that estimates can be obtained for the potency and intrinsic activity of adenosine A2a receptor agonists in vivo by pharmacokinetic-pharmacodynamic analysis of mean arterial pressure data in a rat model. In future studies, total peripheral resistance may also be useful as a pharmacodynamic parameter for A24 activation, provided that possible changes of the stroke volume are also assessed.

Pharmacokinetic-pharmacodynamic modelling of the EEG effects of midazolam in individual rats: influence of rate and route of administration

1. The purpose of the present investigation was to quantify the concentration-pharmacological effect relationship of midazolam in individual rats by use of effect parameters derived from aperiodic EEG analysis. By varying the rate and route of administration the role of (inter)active metabolites and development of acute tolerance was evaluated. 2. The pharmacokinetics and pharmacodynamics of midazolam were determined after intravenous administration of 10 mg kg-1 during 5, 30 and 60 min and oral administration of 15 mg kg-1. Following intravenous administration the pharmacokinetics were most adequately described by a bi-exponential equation. The values (mean +/- s.e. mean, n = 20) of clearance, volume of distribution at steady-state and terminal half-life were 67 +/- 2 ml min-1 kg-1, 1.61 +/- 0.071 kg-1 and 27 +/- 1 min, respectively. Following oral administration midazolam was rapidly absorbed with a systemic availability of 45 +/- 9%. 3. The averaged amplitudes in the 11.5-30 Hz (beta) frequency band of the fronto-central lead on the left-hemisphere, as derived by aperiodic EEG analysis, was selected as a measure of the pharmacological effect of midazolam. By pharmacokinetic-pharmacodynamic modelling the individual concentration-EEG effect relationships of midazolam were derived, which were successfully quantified by the sigmoidal Emax model. No marked and systematic differences in pharmacodynamic parameters were found between the rates and routes of administration. The averaged pharmacodynamic parameters of midazolam obtained after combining the results of all rates and routes of administration were (mean + s.e.mean, n = 27): Eo = 61 + 3puV s 1, Emax = 85 + 3 Vs 1, EC50 = 40 + 3 ngmlP-1 and N = 0.84 + 0.04. 4. The results of the present study show that the concentration-EEG effect relationship of midazolam can be characterized in individual animals using the amplitudes in the 11.5-30 (beta) frequency band as a measure of pharmacological response. Acute tolerance did not develop and (inter)active metabolites did not contribute to this effect parameter within the time span of the experiments.

Personality and benzodiazepine sensitivity in anxious patients and control subjects

Cloninger has recently proposed a model of personality variability that is based on three independent heritable traits of harm avoidance, novelty seeking, and reward dependence, each of which is thought to be mediated by a separate neurochemical and neuroanatomic mechanism. The current study tested hypotheses generated on the basis of this theory in anxious patients and control subjects. Eighteen patients with panic disorder, 12 patients with generalized anxiety disorder, and 21 control subjects underwent both personality testing and assessment of their sensitivity to diazepam, as measured by slowing of saccadic eye movement velocity. As expected, anxious patients displayed higher harm avoidance scores than controls. Although an inverse correlation between harm avoidance and benzodiazepine sensitivity was predicted, no relationship between these variables was found in any diagnostic group. However, a significant correlation was found between novelty-seeking scores and sensitivity to diazepam. This finding, although not predicted by Cloninger's theory, is consistent with prior preclinical and human studies.

Benzodiazepine pharmacodynamics: utility of eye movement measures

The utility of several measures of saccadic and smooth pursuit eye movements as benzodiazepine pharmacodynamic measures was explored in 24 psychiatrically and medically health control subjects. Measures of sedation and memory impairment were also included. Subjects received four logarithmically increasing doses of intravenous diazepam at 15-min intervals on 1 day resulting in monotonically increasing plasma diazepam levels, and placebo on another day in random order 1 week apart. Measures were collected twice at baseline, once after each dose of diazepam/placebo and twice more, 15 and 30 min after the last dose. Peak saccadic velocity and smooth pursuit gain showed the greatest overall and dose-dependent drug effect among eye movement measures. Although effect sizes at the highest dose for memory impairment and self-rated sedation were comparable to these two measures, reliability (i.e., placebo-day fluctuation) with these measures was considerably poorer. Log-linear pharmacodynamic modeling was used to calculate the effective dose (ED30) or concentration (EC30) required to reduce saccadic velocity or pursuit gain by 30%. Almost all (23/24) subjects had linear and easily interpretable plots for saccadic velocity, while a majority (19/24) of subjects had interpretable plots for smooth pursuit gain. The distribution of ED30 and EC30 values showed a wide range of sensitivities to diazepam. These findings suggest that saccadic velocity and smooth pursuit gain are sensitive, reliable, quantitative benzodiazepine pharmacodynamic measures.

Pharmacokinetic-pharmacodynamic relationships in phase I/phase II of drug development

Early investigation of pharmacokinetic-pharmacodynamic relationships in Phase I/II may facilitate the further clinical development of a new drug. Although some pharmacology assessments in Phase I are often only surrogates for the therapeutic effect, PK-PD modelling of those effects provides in general crucial information on the drug's potency in vivo. A mathematical PK-PD expression allows explorative simulations on the rate of onset of drug action, on the intensity and duration of the effects for doses in future clinical trials, or in situations of altered drug kinetics. Furthermore, understanding of the PK-PD relationship early on in drug development may anticipate unnecessary exposure of human subjects to inappropriate drug doses or trials.

Electroencephalogram effect measures and relationships between pharmacokinetics and pharmacodynamics of centrally acting drugs

Electroencephalogram (EEG) effect parameters may be useful in pharmacokinetic-pharmacodynamic modelling studies of drug effects on the central nervous system (CNS). Effect parameters derived from a quantitative analysis of the EEG appear to be perfectly suited to characterise the relationships between pharmacokinetics and pharmacodynamics of benzodiazepines and intravenous anaesthetics. EEG parameters represent many of the characteristics of ideal pharmacodynamic measures, being continuous, objective, sensitive and reproducible. These features provide the opportunity to derive concentration-effect relationships for these drugs in individuals, which yield important quantitative information on the potency and intrinsic efficacy of these drugs. The EEG techniques presented can be used to study the influences of factors such as age, disease, chronic drug use and drug interactions on the concentration-effect relationships of psychotropic drugs. An important issue is the choice of the EEG parameter to characterise the CNS effects of the compounds. More attention must be paid to evaluating the relevance of EEG parameters to the pharmacological effects of the drugs. Knowledge of the relationship between EEG effect parameters and clinical effects of drugs under different physiological and pathophysiological conditions is crucial to determining the value of EEG parameters in drug effect monitoring. Pharmacodynamic parameters derived from the concentration-EEG effect relationship may be correlated to pharmacodynamic parameters obtained from other in vitro and in vivo effect measurements. These comparisons revealed that changes in the amplitudes in the beta frequency band of EEG signals is a relevant measure of pharmacological effect intensity of benzodiazepines, which reflects their affinity and intrinsic efficacy at the central gamma-aminobutyric acid (GABA) benzodiazepine receptor complex. The exact EEG correlates of the anxiolytic, anticonvulsant, sedative and hypnotic actions of benzodiazepines have not yet clearly been elucidated. For intravenous anaesthetics, close correlations between the potency determined with EEG measurements and clinical measures of anaesthetic depth have been established, suggesting that, in principle, EEG parameters can adequately reflect depth of anaesthesia. However, more study is required to further substantiate these findings.

Differences in intrinsic efficacy of benzodiazepines are reflected in their concentration-EEG effect relationship

1. The relevance of EEG effect parameters as a measure of the central nervous system effects of benzodiazepines was evaluated. The concentration-EEG effect relationships of the benzodiazepine agonist midazolam, partial agonist bretazenil, antagonist flumazenil and inverse agonist Ro 19-4603 were quantified and compared with the intrinsic efficacy and affinity of these compounds at the gamma-aminobutyric acid (GABA)-benzodiazepine receptor complex. 2. The pharmacokinetics and pharmacodynamics of the compounds were determined after a single intravenous bolus administration of 5 mg kg-1 midazolam, 2.5 mg kg-1 bretazenil, 10 mg kg-1 flumazenil or 2.5 mg kg-1 Ro 19-4603 to male Wistar derived rats. In a separate experiment the distribution between blood, cerebrospinal fluid and brain concentrations of these compounds was determined. A sensitive assay was developed to measure bretazenil and Ro 19-4603 concentrations in small samples of biological fluids. 3. The benzodiazepine-induced changes in amplitudes in the 11.5-30 Hz frequency band, as determined by aperiodic analysis, was used as EEG effect measure. Concentration-EEG effect relationships were derived by a pharmacokinetic-pharmacodynamic modelling procedure and in the case of midazolam, bretazenil and Ro 19-4603 successfully quantified by the sigmoidal Emax model. Large differences in maximal effect of midazolam (Emax = 73 +/- 2 microVs-1), bretazenil (Emax = 19 +/- 1 microVs-1) and Ro 19-4603 (Emax = -6.5 +/- 0.4 microVs-1) were observed, reflecting their differences in intrinsic efficacy. A close correlation was found between the EC50 values based on free drug concentration and receptor affinity as determined by displacement of [3H]-flumazenil in a washed brain homogenate at 37 degrees C. In the concentration range of receptor saturation flumazenil did not produce any changes in the EEG effect measure.4. The study demonstrated that the change in amplitudes in the 11.5-30 Hz frequency band of the EEG is a relevant measure of the pharmacological effect intensity of benzodiazepines, because it seems to reflect their affinity and intrinsic efficacy at the central GABA-benzodiazepine receptor complex.

The relationship between the concentration of temazepam in cerebrospinal fluid and sedation in man

Twenty-six patients received oral temazepam and subsequently spinal anaesthesia. Blood and lumbar cerebrospinal fluid temazepam levels were measured together with the degree of sedation. The plasma and cerebrospinal fluid concentrations correlated well with the temazepam dose but even better with the weight standardised dose (r = 0.65, p = 0.0003 and r = 0.75, p = 0.00001 respectively). Both the plasma and cerebrospinal fluid concentrations of temazepam were correlated with the patient's sedation (r = 0.42 p = 0.037, and r = 0.46 p = 0.021 respectively), but neither was strong. Thus, although the drug concentration at the receptor may be a major factor in producing sedation, other factors, possibly the receptor population or their responsiveness, are also important contributors.

Biophase equilibration times

Various methods for describing how quickly a drug equilibrates at the biophase are proposed. The biophase equilibration time (BET) is the time it takes the biophase drug level to reach a given percentage (p) of its predicted steady state in a drug administration that leads to a steady-state condition. The time to reach biophase equilibrium may be defined as the BET value for p = 95, and the 50% biophase equilibration time is obtained when p = 50. Biophase equilibration profiles (BEPs), obtained by plotting p versus BET, give a dynamic representation of the approach to equilibrium and may serve as an indicator of the rate of drug delivery to the biophase. A pharmacodynamic system analysis method is proposed to determine BETs and BEPs from the biophase conduction function. The approach is demonstrated using pharmacodynamic data from the CNS effect of amobarbital evaluated by an aperiodic analysis of EEG recordings. The relevance of the BET and/or BEP principles in optimal computer-controlled drug infusion, drug design, and evaluation of targeted drug delivery is discussed. Both vascular and extravascular drug administrations are considered in the analysis.

A system approach to pharmacodynamics. III: An algorithm and computer program, COLAPS, for pharmacodynamic modeling

Many pharmacodynamic (PD) models may be generalized in the form E(t) = N(L[c(t)]), where E(t) is a recorded effect response, c(t) is a sampled drug level, N is a nonlinear autonomic function, and L is a linear operator that commonly is a convolution operation. The NL class of PD models includes the traditional effect compartment PD models as a subclass, but is not limited to such models. An algorithm and computer program named COLAPS, based on system analysis principles and hysteresis minimization, that enable N and L to be empirically determined for the NL class of models without addressing specific kinetic structure aspects ("model independence") are presented. The kinetic concepts of biophase conduction and transduction functions are used by COLAPS. Such an approach is more general than the effect compartment approaches because it does not assume first-order transport principles. The pitfalls of hysteresis minimization in PD modeling are discussed and the procedures taken by COLAPS to avoid these pitfalls are outlined. A transformation technique prevents improper convergence to a point. A novel reparameterization scheme is introduced that maximizes the flexibility of the kinetic functions and extends the generality of the analysis. Inequality function constraints are maintained without the need for troublesome constrained nonlinear optimization procedures. Usage of the COLAPS program is illustrated in the analysis of the PD of amobarbital. The COLAPS program resulted in an excellent minimization of the effect versus biophase level hysteresis. The biophase conduction function, the biophase drug level (normalized), and the transduction curve were determined. The transduction curve showed clear biphasic behavior.

Intra-operative patient-controlled sedation

Patient-controlled sedation using propofol is described in 23 ASA 1 patients undergoing surgical extraction of third molars. Cardiovascular and respiratory stability was maintained and sedation was no deeper than full eyelid closure with prompt response to verbal command. Amnesia for the extractions was common (70%), but there was no significant amnesia for recovery room events. Intra-operative patient cooperation was good, postoperative patient satisfaction high and all patients stated they liked the self-administration and would use the same technique again. Propofol dose was correlated with both procedure length (r = 0.587; p = 0.003) and surgical difficulty (rs = 0.550; p = 0.010).

The extracellular concentration of the anti-epileptic drug valproate in the rat brain as determined with microdialysis and an automated HPLC procedure

The cellular and extracellular brain concentration of valproic acid in freely moving rats has been estimated after intravenous injection of sodium valproate. Some rats were provided with a stereotaxically implanted dialysis probe in the striatum and a cannula in the heart through which the drug was injected and which allowed regular removal of blood. In other rats tissue levels of valproic acid were determined 5 and 90 min after drug injection. Valproic acid was determined by an automated precolumn derivatization procedure followed by HPLC separation and fluorimetric detection. Extracellular concentration was proportional to the blood concentration at every time interval, indicating rapid exchange of the drug between the two compartments. About 50% of the striatal content of valproate was in extracellular space. The experiments demonstrated the usefulness of microdialysis to estimate both the extracellular concentration and the average cellular drug levels, provided a sensitive analysis procedure is available.

The place of simultaneous pharmacokinetic pharmacodynamic modeling in new drug development: trends and perspectives

The potential applications in drug development of pharmacokinetic-pharmacodynamic modeling are numerous: optimal medication regimen, design of galenic forms, identification of specific effects of metabolites and enantiomers. Nevertheless, this methodology is presently under-used, as appears from an analysis of the literature. We examine this point as well as progress in both non-invasive pharmacodynamic measurement and specific experimental design that could lead to a future extension of PK-PD in toxicology and phase I drug development.

Pharmacokinetic-pharmacodynamic modeling of the central nervous system effects of heptabarbital using aperiodic EEG analysis

The concentration EEG effect relationship of heptabarbital was modeled using effect parameters derived from aperiodic EEG analysis. Male Wistar rats (n = 10) received an intravenous infusion of heptabarbital at a rate of 6-9 mg/kg per min until burst suppression with isoelectric periods of 5 sec or longer. Arterial blood samples were obtained and EEG was measured continuously until recovery of baseline EEG and subjected to aperiodic analysis for quantification. Two EEG parameters, the amplitudes per second (AMP) and the total number of waves per second (TNW), in five discrete frequency ranges and for two EEG leads were used as descriptors of the drug effect on the brain. The EEG parameters responded both qualitatively and quantitatively different to increasing concentrations of heptabarbital. Monophasic concentration effect curves (decrease) were found for the frequency ranges greater than 2.5 Hz and successfully quantified with a sigmoidal Emax model after collapsing the hysteresis by a nonparametric modeling approach. For the parameter TNW in the 2.5-30 Hz frequency range the value of the pharmacodynamic parameters EC50, Emax, and n (means +/- SD) were 78 +/- 7 mg/L, 11.4 +/- 1.7 waves/sec and 5.0 +/- 1.5, respectively. For other discrete frequency ranges, differences in EC50 were observed, indicating differences in sensitivity to the effect of heptabarbital. In the 0.5 +/- 2.5 Hz frequency range biphasic concentration effect relationships (increase followed by decrease) were observed. To fully account for the hysteresis in these concentration effect relationships, postulation of two effect compartments was necessary. To characterize these biphasic effect curves two different pharmacodynamic models were evaluated. Model 1 characterized the biphasic concentration effect relationship as the summation of two sigmoidal Emax models, whereas Model 2 assumed the biphasic effect to be the result of only one inhibitory mechanism of action. With Model 1 however realistic parameter estimation was difficult because the maximal increase could not be measured, resulting in high correlations between parameter estimates. This seriously limits the value of Model 1. Model 2 involves besides estimation of the classical pharmacodynamic parameters Emax, EC50, and n also estimation of the maximal disinhibition Amax. This model is a new approach to characterize biphasic drug effects and allows, in principle, reliable estimation of all relevant pharmacodynamic parameters.

Pharmacodynamics of tolerance development to the anesthetic and anticonvulsant effects of phenobarbital in rats

In this study, the potential development of tolerance towards the anesthetic and anticonvulsant effects of phenobarbital after chronic administration to rats was investigated, differentiating between dispositional and functional tolerance. Chronic exposure to phenobarbital by daily ip injections of 20 or 100 mg/kg for 14 days caused a 30% increase in clearance within 1 week. No changes occurred in the total amounts of phenobarbital and para-hydroxyphenobarbital excreted into urine and feces. Pharmacodynamic effects were quantitated after 2 weeks of treatment. The anesthetic effect was measured by slow iv infusion of phenobarbital until onset of loss of righting reflex (LRR), followed by measurement of drug concentrations in serum (both total and free), brain, and cerebrospinal fluid (CSF). With a phenobarbital dose of 100 mg/kg daily, CSF concentrations at the onset of LRR significantly increased from 136 +/- 12 to 176 +/- 21 mg/L (p less than 0.001), indicating that functional tolerance developed for the anesthetic effect. Protection of phenobarbital against convulsions induced by pentylenetetrazol (PTZ), as measured by the elevation of the PTZ plasma threshold concentration necessary to elicit seizures at the EC50 of phenobarbital, was not altered. The discrepancy observed in the development of functional adaptation of the CNS demonstrates that different mechanisms of action are reflected in the different measures of the anesthetic and anticonvulsant effects of phenobarbital that were utilized.

Pharmacokinetic-pharmacodynamic modelling of the anticonvulsant effect of oxazepam in individual rats

1. The purpose of this investigation was to examine in vivo drug-concentration anticonvulsant effect relationships of oxazepam in individual rats following administration of a single dose. 2. Whole blood concentration vs time profiles of oxazepam were determined following administration of doses of 4, 8 and 12 mg kg-1. The pharmacokinetics could be described by an open 2-compartment pharmacokinetic model. Following 12 mg kg-1 the values (mean +/- s.e., n = 11) of clearance and volume of distribution were 28 +/- 2 ml min-1 kg-1 and 2.6 +/- 0.31 kg-1, respectively, and were not significantly different from the values obtained at the other doses. 3. The anticonvulsant effect was quantitated by a new technique which allows repetitive determination of the convulsive threshold by direct cortical stimulation within one rat. Significant dose-dependent elevations of the seizure threshold were observed. 4. By pharmacokinetic-pharmacodynamic modelling, a log-linear relationship was found between concentration and anticonvulsant effect. Following 12 mg kg-1 the values (mean +/- s.e., n = 11) of the pharmacodynamic parameters slope and minimal effective concentration (Cmin) were 243 +/- 27 microA and 0.11 +/- 0.02 mg l-1, respectively and not significantly different from the values obtained at the other doses. 5. In a repeatability study the pharmacodynamic parameters were determined twice on two different occasions with an interval of two weeks in the same group of 11 rats. The inter-animal variability in the pharmacodynamic parameter slope was 46%, whereas the intra-animal variability was 24 +/- 18%. The value of the minimal effective concentration was in each animal and on each occasion close to zero within the relatively narrow range of 0.01-0.30mgI. 6. The results of this study showed that it is possible to determine in vivo concentration-anticonvulsant effect relationships of oxazepam under non-steady-state conditions in individual rats. The anti-convulsant effect of oxazepam appeared to be a rapidly reversible direct effect and acute tolerance did not develop within the time frame of the experiments.