![]() |
|
|
Vol. 303, Issue 3, 1130-1137, December 2002
Leiden/Amsterdam Center for Drug Research, Division of Pharmacology, Gorlaeus Laboratories, Leiden, The Netherlands (K.P.Z., J.R.-P., H.J.M., P.H.V., M.D.); and Mathematical Institute, Leiden University, Niels Bohrweg, Leiden, The Netherlands (L.A.P.)
| |
Abstract |
|---|
|
|
|---|
The objective of this investigation was to compare the in vivo potency and intrinsic activity of buspirone and its metabolite 1-(2-pyrimidinyl)-piperazine (1-PP) in rats by pharmacokinetic-pharmacodynamic modeling. Following intravenous administration of buspirone (5 or 15 mg/kg in 15 min) or 1-PP (10 mg/kg in 15 min), the time course of the concentrations in blood were determined in conjunction with the effect on body temperature. The pharmacokinetics of buspirone and 1-PP were analyzed based on a two-compartment model with metabolite formation. Differences in the pharmacokinetics of buspirone and 1-PP were observed with values for clearance of 13.1 and 8.2 ml/min and for terminal elimination half-life of 25 and 79 min, respectively. At least 26% of the administered dose of buspirone was converted into 1-PP. Complex hypothermic effects versus time profiles were observed, which were successfully analyzed on the basis of a physiological indirect response model with set-point control. Both buspirone and 1-PP behaved as partial agonists relative to R-(+)-8-hydroxy-2-(di-n-propylamino)tetralin (R-8-OH-DPAT) with values of the intrinsic activity of 0.465 and 0.312, respectively. Differences in the potency were observed with values of 17.6 and 304 ng/ml for buspirone and 1-PP, respectively. The results of this analysis show that buspirone and 1-PP behave as partial 5-hydroxytryptamine1A agonists in vivo and that following intravenous administration the amount of 1-PP formed is too small to contribute to the hypothermic effect.
| |
Introduction |
|---|
|
|
|---|
Buspirone
is a selective 5-HT1A agonist that is used clinically as an
anxiolytic drug (for reviews on its pharmacological properties, see
New, 1990
and Fulton and Brogdon, 1997
).
In vivo buspirone is metabolized into 1-(2-pyrimidinyl)-piperazine
(1-PP), which also possesses affinity to the 5-HT1A
receptor. This metabolite could therefore contribute to the effect of
buspirone (Bianchi et al., 1988
; Manahan-Vaughan
et al., 1995
; Cao and Rodgers, 1997
). At present
there is limited quantitative information on the magnitude of this
effect. In theory the contribution of the metabolite to the effect is
determined by the relative concentrations of buspirone and 1-PP, their
relative in vivo potency and intrinsic activity, and the nature of the
pharmacodynamic interaction between the two.
Recently we have developed a physiological
pharmacokinetic-pharmacodynamic (PK-PD) model to quantitatively
characterize the pharmacodynamics of 5-HT1A receptor
agonists in vivo using the hypothermic effect as a pharmacodynamic
endpoint (Zuideveld et al., 2001
,
2002a
). The hypothermic
response is considered a robust marker of 5-HT1A activity
(Millan et al., 1993
; Cryan et al., 1999
). It has been shown that the PK-PD model allows for the
estimation of the in vivo potency and intrinsic activity of a wide
array of different 5-HT1A agonists such as R-8-OH-DPAT,
S-(
)-8-hydroxy-2-(di-n-propylamino)tetralin (S-8-OH-DPAT),
N-[2-[4-(2-methoxyphenyl)-1-piperzinyl]ethyl]-n-2-pyridinyl-cy-clohexanecarboxamide (WAY-100,635), and Flesinoxan (Zuideveld et al., 2001
,
2002a
, 2002b
). In the present investigation we apply this model to
determine the relative in vivo potency and intrinsic activity of
buspirone and its metabolite 1-PP, taking into consideration the
potential pharmacodynamic interaction between the two. In this respect
it is important to consider the mechanism of the hypothermic response of buspirone and 1-PP. The hypothermic effect of buspirone has been
widely studied (Higgins et al., 1988
; Koenig et
al., 1988
; Young et al., 1993
) and is generally
considered a specific 5-HT1A receptor-mediated response
(Koenig et al., 1988
; Kommisarev et al.,
1990
; Cryan et al., 1999
). 1-PP however has
affinity for both the 5-HT1A receptor
(KD = 1035 nM = 94 ng/ml
1) and the
2-adrenoceptor
(
2-AR) (KD = 40 nM = 3.6 ng/ml
1) (Gobert et al., 1995
) and it
is therefore possible that this effect is not mediated by the
5-HT1A receptor (Kommisarev et al., 1990
).
| |
Materials and Methods |
|---|
|
|
|---|
Experiments that were approved by the Ethics Committee of the University of Leiden were performed on male Wistar rats (Broekman BV, Someren, The Netherlands) weighing 315 ± 3 g (mean ± S.E.M., n = 27). They were housed in standard plastic cages (six per cage before surgery and individually after surgery). They were housed under normal 12-h light/dark cycle (lights on at 7 AM and lights off at 7 PM) and a temperature of 21°C. During the light period, a radio was turned on for background noise. Acidified water and food (laboratory chow; Hope Farms, Woerden, The Netherlands) were provided ad libitum before the experiment.
Surgical Procedure
Eight days before the experiment, the rats were operated upon. The animals were anesthetized with an intramuscular injection of 0.1 ml/kg Domitor (1 mg/ml medetomidine hydrochloride; Pfizer, Capelle a/d IJsel, The Netherlands) and 1 ml/kg Ketalar (50 mg/ml ketamine base; Pfizer, Hoofddorp, The Netherlands). Indwelling pyrogen-free cannulae for drug administration were implanted into the right jugular vein (Polythene, 14 cm, 0.52-mm i.d., 0.96-mm o.d.) for infusions of buspirone and 1-PP. For blood sampling, the left femoral artery (Polythene, 4 cm, 0.28-mm i.d., 0.61-mm o.d. + 20 cm, 0.58-mm i.d., 0.96-mm o.d.) was cannulated. Cannulae were tunneled subcutaneously to the back of the neck and exteriorized. To prevent coagulation of blood, the cannulae were filled with a 25% (w/v) solution of polyvinylpyrrolidone (PVP) (Brocacef, Maarssen, The Netherlands) in a 0.9% (w/v) pyrogen-free sodium chloride solution (NPBI, Emmer-Compascuum, The Netherlands) that contained 50 IU/ml heparin (Leiden University Medical Center, Leiden, The Netherlands). Just before the experiment, the PVP solution was removed and the cannulae were flushed with saline containing 20 IU/ml heparin. The skin in the neck was stitched with normal sutures, and the skin in the groin was closed with wound clips. A telemetric transmitter (Physiotel implant TA10TA-F40 system; Data Sciences International, St. Paul, MN) with a weight of approximately 7 g, which had been made pyrogen free with CIDEX (22 g/l glutaraldehyde; Johnson and Johnson Medical Ltd., Gargrave, Skipton, U.K.) for at least 2 h, was implanted into the abdominal cavity for the measurement of core body temperature. After surgery, an injection of the antibiotic ampicillin (0.6 ml/kg of a 200 mg/ml, A.U.V., Cuijk, The Netherlands) was administered.
Experimental Protocol
Dosage Regimen. The experiments were performed 8 days after surgery. Rats received infusions of buspirone and 1-PP. For buspirone, infusions of 5 and 15 mg/kg in 15 min (12.96 and 38.91 µmol/kg for n = 6 and n = 7, respectively) were given. For 1-PP, an infusion of 10 mg/kg (60.90 µmol/kg for n = 6) in 15 min was administered. Six rats received vehicle treatments in which an equivalent amount of saline was infused. For the infusions, an external cannula with a specific volume was filled with a solution of the drug in an amount of saline, which was calculated according to the weight of the rat and connected to the infusion pump (BAS BeeHive, BAS Bioanalytical Systems Inc., West Lafayette, IN). All the experiments started between 9:00 and 9:30 AM.
Blood Sampling.
Approximately 15 to 18 serial blood samples
of 50 µl were taken according to a fixed time schedule to determine
the pharmacokinetics of the drug. The exact amount was measured with a
capillary (Servoprax, Wesel, Germany) and transferred into a glass
centrifuge tube containing 400 µl of Millipore water for hemolysis.
During the experiment the samples were kept on ice. After the
experiment, samples were stored at
20°C pending analysis.
Data Acquisition
Temperature Measurements. To measure the body temperature of the rat a telemetric system (Physiotel Telemetry system; Data Sciences International) was used. A telemetric transmitter (Physiotel implant TA10TA-F40; Data Sciences International) was implanted in the abdominal cavity of the rat. The transmitter measured the body temperature every 30 s for a 2-s period and signaled the data to a receiver (Physiotel receiver, model RPC-1; Data Sciences International). The receiver was connected to the computer through a BCM 100 consolidation matrix (Data Sciences International). The temperature profiles and the room temperature data (C10T temperature adapter; Data Sciences International) were processed and visualized using the Dataquest LabPro software (Data Sciences International, running under OS/2 Warp, IBM).
HPLC Analysis of Buspirone and 1-PP.
The blood
concentrations of buspirone and 1-PP were determined simultaneously by
HPLC using UV detection. The HPLC system consisted of a Shimadzu LC-10
AD pump (flow rate: 1 ml/min; Shimadzu Corporation, Kyoto, Japan), a
Waters 712WISP autosampler (Waters, Milford, MA), a Waters Sperisorb S5
OD/CN column (mixed phase 4.6 mm × 250 mm) together with a
Sperisorb S5CN guard column (all Waters), a Spectroflow 757 UV
absorbance detector set at 240 nm (Kratos Analytical, Ramsey, NJ), and
a C-R3A Chromatopac integrator (Shimadzu Corporation). The mobile phase
consisted of a mixture of a 0.05 M potassium phosphate buffer (pH 5)
and acetonitrile (60:40, v/v), to which 0.005 M KCl was added. The
analytes were extracted from blood using an extensive liquid-liquid
extraction to allow for the simultaneous analysis of both buspirone and
1-PP. To 50 µl of blood hemolyzed in 400 µl of water, 100 µl of
the internal standard solution (5000 ng/ml
4-[3-(benzotriazol-1-yl)propyl]-1-(2-methoxyphenyl)-piperazine) was
added. The mixture was deproteinized by adding 2 ml of acetonitrile and
centrifuged. The supernatant was transferred to a clean tube and
volumes of 50 µl of hydrochloric acid (1 M), 1 ml of water and 3 ml
of distilled ethyl acetate were added. After separation, 2 ml of borate
buffer (0.2 M, pH 10) and 3 ml of dichloromethane was added to the
aqueous solution. The organic phase was transferred to a clean tube and
evaporated to dryness under vacuum at 37°C. The residue was
re-dissolved in 100 µl of 30% acetonitrile, of which 50 µl was
injected into the HPLC system. On the day of analysis, a 9-point
calibration curve was prepared for both buspirone and 1-PP by spiking
50 µl of blood hemolyzed in 400 µl of water with 50 µl of
buspirone and 1-PP. This resulted in a blood concentration range of 10 to 10,000 ng/ml for buspirone and 1 to 10,000 ng/ml for 1-PP. Samples
were processed as described, and the peak-area ratio of either
buspirone or 1-PP over the internal standard were calculated.
Calibration curves were constructed by weighted linear regression
([weight factor = 1/(peak-area ratio)]2). Recovery,
intra-, and inter-assay variation were determined using spiked blood of
500 and 1500 ng/ml for buspirone and 50, 500, and 1500 ng/ml for 1-PP.
For buspirone, recovery was (n = 3, mean ± S.E.M., corrected for volume loss) 41 ± 3.6% and 45 ± 8.6% for 500 and 1500 ng/ml, respectively. For 1-PP, recovery was
(n = 3, mean ± S.E.M., corrected for volume loss)
88 ± 9.8%, 82 ± 3.8%, and 88 ± 16.5% for 50, 500, and 1500 ng/ml, respectively. Inter-assay variation (n = 10) was 9.9% and 19.2% (accuracy:
8.9% and 1.2%) for 500 and
1500 ng/ml buspirone and 9.3, 14.4, and 16.3% (accuracy:
12.6,
0.4, and
0.8%) for 50, 500, and 1500 ng/ml 1-PP. Intra-day assay
variation (n = 4) was 4.4 and 3.4% (accuracy:
2.5
and 1.4%) for 500 and 1500 ng/ml buspirone and 7.4, 9.8, and 7.9%
(accuracy: 12.6,
1.7, and
17.3%) for 50, 500, and 1500 ng/ml 1-PP.
Detection limit (signal-to-noise ratio 3) using 50 µl of blood was
22.2 and 5.2 ng/ml for buspirone and 1-PP, respectively.
Chemicals. Buspirone and 1-PP were generously donated by Bristol-Myers Squibb (Princeton, NJ). 4-[3-(Benzotriazol-1-yl)propyl]-1-(2-methoxyphenyl)-piperazine maleate was purchased from Tocris (Lanford, Bristol, UK). Dichloromethane was purchased from Riedel-de Haën (Seelze, Germany). All other chemicals used were of analytical grade (Baker, Deventer, The Netherlands).
Data Analysis
A population approach was used to quantify both the pharmacokinetics and pharmacodynamics of buspirone and 1-PP. Using this approach, the population is taken as the unit of analysis while taking into account both intra-individual variability in the model parameters as well as inter-individual residual error. Modeling was performed using the nonlinear mixed effects modeling software NONMEM developed by Sheiner & Beal (version V 1.1, NONMEM project Group, University of California, San Francisco, CA). Individual predictions were obtained in a Bayesian post hoc step.
Pharmacokinetic Analysis.
The concentration-time profiles of
the buspirone administration, the 1-PP formed during the buspirone
administration, and the 1-PP following direct administration were best
described using a two-compartment model with metabolite formation. In
this model, the input function for the metabolite following
administration of the parent drug is the fraction of buspirone that is
converted into 1-PP. This model is mathematically described as follows
(Houston, 1985
)
|
(1) |
|
(2) |
|
1-PP in eq. 1 represents the clearance for the
conversion of buspirone into 1-PP. The value of this clearance can be
used to determine the fraction of buspirone that is converted into
1-PP, according to
|
(3) |
|
(4) |
is the population value for parameter P,
Pi is the individual value and
i is the
random deviation of Pi from P. The
values of
i are assumed to be independently normally
distributed with mean zero and variance
2. Residual
error was characterized by a proportional error model,
|
(5) |
accounts for the residual deviance of the model predicted value from
the observed concentration. Different
values were assigned for the
observations of buspirone, 1-PP during buspirone infusions, and 1-PP
when modeled simultaneously. The values for
are assumed to be
independently normally distributed with mean zero and variance
2. The values for the population
,
2,
and
2 are estimated using the first order method in
NONMEM.
|
Pharmacodynamics Analysis.
Recently we have proposed a
physiological model, to describe the complex time course of the
hypothermic response in vivo following the administration of
5-HT1A receptor agonists (Zuideveld et al., 2001
). The model contains a set-point control that can be
attenuated by the drug receptor interaction and utilizes the concept of
an indirect physiological response model (Dayneka et al.,
1993
). This model considers a 0th-order rate constant
associated with the warming of the body (kin)
and a first-order rate constant associated with the cooling of the body
(kout). The model is further characterized by
the phenomenological system parameter A and a parameter
, which
determines the amplification in the system (for details, see
Zuideveld et al., 2001
). In this model, the stimulus
generated by the drug-receptor interaction, which drives the
physiological processes that lower temperature and characterize the
drug-receptor interaction in terms of the drugs' potency and intrinsic
activity, is described by the function f(C),
|
(6) |
|
(7) |
|
(8) |
|
(9) |
|
(10) |
accounts for the residual deviance of the model predicted value from the observed value. The values for the population
,
2, and
2 are estimated using the centering first-order
conditional estimation method with the first-order model in NONMEM. A
conditional estimation method was used due to the high degree of
nonlinearity of the model and the high density of the data. The
centering option gives the average estimate of each element of
together with a P value that can be used to assess whether
this value is sufficiently close to zero. The occurrence of an average
that is significantly different from zero indicates an uncentered
or a biased fit. This method was not chosen because the average
estimates of each element of
were expected to be different from
zero, but rather to greatly decrease computing time as required with
just the conditional estimation method. To further decrease computing
time only 1/16th of the temperature data set was used for modeling,
reducing the temperature measurements from over 900 measurements per
individual to approximately 60. The implication of this reduction is
that there is a data point every 8 min, as opposed to every 0.5 min. This reduction did not void the integrity of the data profiles.
Statistical Analysis.
Goodness-of-fit was analyzed using the
objective function and various diagnostic methods. Model selection was
based on the Akaike Information Criterion (Akaike, 1974
)
and assessment of the accuracy and precision of parameter estimates and
correlations between the parameter estimates. Comparisons of maximal
decreases in body temperature were performed using an unpaired
t test, where a P value <0.05 was considered significant.
| |
Results |
|---|
|
|
|---|
Pharmacokinetics.
During the experiments blood samples were
taken to construct individual concentration-time profiles for both
buspirone and 1-PP. Individual concentration profiles were predicted by
fitting a population pharmacokinetic model to the data (eq. 1). On the basis of the concentration curves, the goodness-of-fit plots and the
minimum value of the objective function, two-compartment models were
selected for buspirone, 1-PP, and 1-PP following administration of
buspirone. The concentration-time profiles of the 1-PP administration were fitted both separately and simultaneously with the "1-PP during
the buspirone administration" profiles. Ultimately the values of the
pharmacokinetic parameters of 1-PP were fixed to those obtained upon
the separate administration. These values were subsequently entered as
constants in the analysis of the concentration profiles of buspirone
and 1-PP on the basis of the two-compartment metabolite model (eq. 1),
which resulted in a significant improvement of the objective function.
The obtained pharmacokinetic parameter estimates were used as input in
the pharmacodynamic analysis because it is believed that these
concentrations are closest to the "true" concentrations at each
observation of the body temperature. Both the individually measured
concentrations and the population predicted profiles following
administration of buspirone and 1-PP are depicted in Fig.
2. All parameters were estimated in their
mixed effect form, the random effect being incorporated in exponential
error models. The results of the analysis are represented in Table
1. The clearance found for buspirone was
13 ml/min resulting in a half-life of 25 min. For 1-PP, a clearance of
8.2 ml/min was found, which corresponds to a half-life of 79 min. Based
on the pharmacokinetic analysis (eq. 3), 26% of the administered dose
of buspirone is metabolized to 1-PP. Pharmacokinetic parameter
estimates were independent of administered dose and no statistically
significant correlations between the parameter estimates and between
the parameters from the different dosing groups were detected. Table 1
further displays the inter-individual variation, and the
intra-individual variation, which are both found to be reasonable for
both buspirone and 1-PP. For the different drugs and administrations,
different intra-individual variation was predicted, which significantly
improved the fit. The precision of the parameters prediction, as
represented by their 95% confidence intervals, is good.
|
|
Pharmacodynamics.
The average time-effect profiles for the
effect on body temperature following administration of vehicle,
buspirone, and 1-PP are depicted in Fig.
3. The average baseline temperature
(S.E.M., n = 33) was 38.0 ± 0.005°C. Following
administration of 5 and 15 mg/kg buspirone in a 15-min infusion, a
maximal decrease of 2.8 ± 0.3°C was observed after
approximately 55 min. Administration of 10 mg/kg 1-PP in a 15-min
infusion resulted in a maximum temperature decrease of 1.6 ± 0.2°C after approximately 36 min. The temperature decreases were all
statistically significantly (P < 0.05) different from
the vehicle treatment, consisting of saline. After reaching a maximal
decrease, a rapid recovery was observed followed by a plateau phase
before the body temperature returned to baseline.
|
). These parameters
were estimated using the centering first-order conditional estimation
method with a first-order model. The average values for all
values
were not significantly different from 0. Population parameter estimates are depicted in Table 2. Individual post
hoc predictions of the parameters were not biased with respect to the
different treatment groups. Interestingly, no difference was found in
the pharmacodynamic parameter estimates for buspirone regardless of
whether the sigmoid pharmacodynamic model (eq. 6) or the noncompetitive
interaction model (eq. 8) were used. The parameters obtained with the
competitive interaction model (eq. 9) differ substantially, in
particular with respect to the values of SC50
and Smax. No difference was found in the
system-specific parameters between the various models or various
compounds. The inter-individual variations of the pharmacodynamic parameter estimates for 1-PP were larger than for buspirone. Figure 4 depicts the individual observations,
predictions, and the population predictions for each model and two
typical individual rats.
|
|
| |
Discussion |
|---|
|
|
|---|
The present study provides novel information on the pharmacokinetic and pharmacodynamic relationship of buspirone and its metabolite 1-PP in vivo. Specifically, the present investigation has also focused on the relative potency and intrinsic activity of the metabolite 1-PP and the extent to which 1-PP contributes to the effect of the buspirone. This study shows that with values of the relative intrinsic activity of 0.465 and 0.312 both buspirone and 1-PP behave as partial agonists relative to R-8-OH-DPAT. Furthermore it is shown that the potency of 1-PP is 20-fold lower compared with buspirone and that following intravenous administration of the parent drug, formation of 1-PP does not contribute to the observed effect.
In comparative pharmacodynamic investigations in vivo, it is important
to take differences in pharmacokinetics between the compounds into
account. Therefore in the present investigation, the effects of
buspirone and 1-PP were determined by an integrated PK-PD approach. In
the present study, the value of the clearance was 13 ml/min resulting
in a half-life of 25 min for buspirone and a clearance of 8.2 ml/min
with a corresponding half-life of 79 min for the metabolite 1-PP. Both
values are in close agreement with the values reported by Caccia
et al. (1985)
. In the present study, it was estimated on the
basis of the two-compartment model with metabolite formation that at
least 26% of the administered dose of buspirone is metabolized into
1-PP. Again, this appears to be consistent with findings of
Caccia et al. (1986)
and Jajoo et al.
(1989a
, 1989b
) who
showed that 25% of the administered dose is excreted into urine as
1-PP. Overall the pharmacokinetics of both buspirone and 1-PP could be
characterized successfully using the metabolite model.
To quantify the pharmacodynamic response of both buspirone and 1-PP the
hypothermic response was chosen as a pharmacodynamic endpoint.
Recently, we have developed an integrated
pharmacokinetic-pharmacodynamic model to characterize
5-HT1A receptor-mediated hypothermia, in terms of potency
and intrinsic activity of the agonists used (Zuideveld et al.,
2001
). The 5-HT1A agonist-induced hypothermic
response was chosen, because it is considered to be a robust marker for 5-HT1A activity (Millan et al., 1993
;
Cryan et al., 1999
) and behaves as an ideal biomarker,
which is continuous, reproducible, and sensitive enough to discriminate
between a full and a partial agonist (Hadrava et al.,
1996
) and is selective (Millan et al., 1993
;
Cryan et al., 1999
). The PK-PD model itself combines an indirect physiological response model (Dayneka et al.,
1993
) with a set-point control. It is believed that the
5-HT1A receptors play a role in maintaining the body's
set-point temperature (Lin et al., 1983
; Schwartz
et al., 1998
; Zeisberger, 1998
). In addition numerous reports suggest that this set-point is regulated through an
interplay between the 5-HT1A (hypothermia) and the
5-HT2A/C (hyperthermia) receptor system (Gudelsky et
al., 1986
; Salmi and Ahlenius, 1998
). Therefore
the proposed model is considered to reflect physiology. This is
confirmed by the observation that a competitive interaction model
accurately characterizes the interaction between R-8-OH-DPAT and
WAY-100,635 (Zuideveld et al., 2002a
). In addition, this
model has been applied successfully to characterize the in vivo
pharmacodynamics of a range of 5-HT1A ligands including R-
and S-8-OH-DPAT, WAY-100,635 and flesinoxan (Zuideveld et al., 2001
, , 2002b
) showing
differences in both in vivo potency as well as intrinsic activity.
Regarding the hypothermic response induced by buspirone, much evidence
has been presented indicating that it is mediated through the
5-HT1A receptor (Higgins et al., 1988
;
Young et al., 1993
). Furthermore its effect can be
antagonized by selective 5-HT1A receptor antagonists
(Koenig et al., 1988
; Cryan et al., 1999
) and its selectivity over other receptors that may also mediate a
hypothermic response, e.g., the
2-AR (Dilsaver et
al., 1989
; MacDonald et al., 1989
), µ-opioid
(OP3) (Spencer et al., 1988
), dopamine D3
(Barik and de Beaurepaire, 1998
) and adenosine
A1 receptors (Malhotra and Gupta, 1997
) is
high (New, 1990
; Gobert et al., 1995
).
1-PP on the other hand has affinity for both the 5-HT1A
(KD = 1035 nM = 94 ng/ml) and the
2-AR (KD = 40 nM = 3.6 ng/ml) (Gobert et al., 1995
). The mechanism of the
hypothermic response of 1-PP is therefore subject to debate. Despite
the fact that 1-PP behaves as an antagonist (rather than an agonist) at the
2-AR, evidence has been presented that suggests that
its hypothermic response may not be mediated by 5-HT1A
receptors (Komissarev et al., 1990
). In this respect it
is of interest that in a separate set of experiments it has been shown
that the selective 5-HT1A receptor agonist WAY-100,635,
when administered in doses that inhibit the effect of 8-OH-DPAT in a
competitive manner, does not influence the hypothermic response of 1-PP
(Zuideveld et al. unpublished observations). This shows
that 1-PP apparently induces hypothermia via a mechanism that is
different from that of buspirone and other selective 5-HT1A agonists.
In the present study, the PK-PD relationships of the hypothermic effect
for both buspirone and 1-PP are characterized in a quantitative manner
using the set-point model developed for characterization of
5-HT1A receptor agonists. Although the mechanism of the
1-PP-induced hypothermic response is unclear, the set-point model could
be applied successfully to describe the 1-PP-induced hypothermic response. This seems justified since the observed complex time-effect profile bears a great similarity to that following administration of a
selective 5-HT1A agonist. Furthermore, it is well
established that the mechanism of
2-AR-mediated
hypothermia is very similar to that of 5-HT1A-mediated
hypothermia. Both receptors are connected to the regulation of the
body's set-point and utilize the same effector mechanisms
(Frank et al., 1997
; Sallinen et al.,
1997
). Pharmacodynamic analysis revealed that the hypothermic
response of 1-PP could be well characterized using this model, showing that 1-PP behaves as a partial agonist with a value of the parameter Smax of 0.312 relative to R-8-OH-DPAT.
Furthermore the potency is roughly 20 times lower than that of
buspirone itself.
In the pharmacodynamic analysis, we have utilized both a non- and a
competitive interaction model to describe the possible interaction
between buspirone and 1-PP. The noncompetitive interaction model as
proposed by Ariens and Simonis (1964)
was chosen because of its simplicity, and since the maximal stimulus of the system is
known, as defined by the maximal response of R-8-OH-DPAT, the parameters are identifiable. A specific feature of this model is that
it assumes the interaction to be additive, with the consequence that
neither synergism nor antagonism can be estimated. Therefore, the
pharmacodynamics of buspirone has also been characterized with a
regular sigmoidal model, disregarding 1-PP. As can be concluded from
the pharmacodynamic analysis, the parameter estimates for potency,
intrinsic activity, and the slope factor are very similar between the
interaction model and the regular sigmoidal model, suggesting that 1-PP
hardly plays a significant role in buspirone's hypothermic response,
essentially adding zero to the effect. The data was also analyzed using
a competitive interaction model as proposed by Holford and
Sheiner (1981)
. The major disadvantage of this model is that
essentially the slope factor n determines whether either
synergism or antagonism is predicted. Despite the fact that the
parameters obtained with the competitive interaction model deviate the
most from the ones found with the noncompetitive and sigmoidal model,
the slope factors are relatively close to 1, and it is concluded that
the competitive interaction model is more sensitive to a relatively
small amount of 1-PP in the body. The peak concentration of 1-PP during
the 15 mg/kg dose of buspirone is 210 ± 34 ng/ml after 60 min,
which only represents a mere 69% of its SC50,
and therefore too low to contribute significantly. However, as
buspirone is subjected to large scale presystemic metabolism
(Caccia et al., 1985
), which is expected given its high
clearance (present investigation), it should be anticipated that in
humans higher concentrations of 1-PP will be found after an oral
administration, which could ultimately interfere with buspirone effect
in a significant manner. In this respect it is interesting to note that
based on the potency and its ratio with its affinity for the
5-HT1A receptor [potency, 304 ng/ml]/[affinity, 94 ng/ml] = 3.2 for 1-PP versus [potency, 17.6 ng/ml]/[affinity, 2.7 ng/ml] = 6.5 for buspirone, it is conceivable that 1-PP mediates its
hypothermic effect through the 5-HT1A receptor.
In the present study, the PK-PD relationship of the hypothermic effect for both buspirone and 1-PP is characterized in a quantitative manner using the set-point model. The results show that buspirone and 1-PP behave as partial 5-HT1A agonists in vivo relative to 5-HT1A agonists such as R-8-OH-DPAT and flesinoxan. Pharmacodynamic analysis using both a sigmoidal Emax, a non- and a competitive interaction model revealed that in the rat, the peak 1-PP concentrations reached are too low to contribute significantly to the buspirone effect after intravenous administration.
| |
Acknowledgments |
|---|
We thank Erica Tukker for technical assistance in the animal experiments. The generous donation of buspirone and 1-PP by Bristol-Myers Squibb is highly appreciated.
| |
Footnotes |
|---|
Accepted for publication July 3, 2002.
Received for publication April 4, 2002.
1 Present address: Pharsight Corporation, Argentum, 2 Queen Caroline St., Hammersmith, W6 9DT London, UK.
2 Present address: Pliva Ïstra ivaèki Institut, Division of Pharmacology, Prilaz Baruna Filipoviaæ 25, 10000 Zagreb, Croatia.
3 Pfizer Global Research & Development, Discovery Biology, Ramsgate Road, Sandwich, Kent CT13 9NJ, UK.
DOI: 10.1124/jpet.102.036798
Address correspondence to: Meindert Danhof, Leiden/Amsterdam Center for Drug Research, Division of Pharmacology, Gorlaeus Laboratory, P.O. Box 9502, 2300 RA Leiden, The Netherlands. E-mail: m.danhof{at}lacdr.leidenuniv.nl
| |
Abbreviations |
|---|
5-HT1A, 5-hydroxytryptamine1A;
1-PP, 1-(2-pyrimidinyl)-piperazine;
PK-PD, pharmacokinetic-pharmacodynamic;
R-8-OH-DPAT, R-(+)-8-hydroxy-2-(di-n-propylamino)tetralin;
S-(
)-8-hydroxy-2-(di-n-propylamino)tetralin, WAY-100,635,
N-[2-[4-(2-methoxyphenyl)-1-piperazinyl]ethyl]-n-2-pyridinyl-cyclohexanecarboxamide;
PVP, polyvinylpyrrolidone;
HPLC, high performance liquid
chromatography;
2-AR,
2-adrenergic
receptor.
| |
References |
|---|
|
|
|---|
2-adrenoceptor antagonist activity of ipsapirone and gepirone is mediated by their common metabolite 1-(2-pyrimidinyl)-piperazine (PmP).
Eur J Pharmacol
151:
365-371[CrossRef][Medline].
2C-adrenoceptor expression in mice: influence on locomotor, hypothermic, and neurochemical effects of dexmedetomidine, a subtype-nonselective
2-adrenoceptor agonist.
Mol Pharmacol
51:
36-46This article has been cited by other articles:
![]() |
H. Wong, R. C. Dockens, L. Pajor, S. Yeola, J. E. Grace Jr., A. D. Stark, R. A. Taub, F. D. Yocca, R. C. Zaczek, and Y.-W. Li 6-Hydroxybuspirone Is a Major Active Metabolite of Buspirone: Assessment of Pharmacokinetics and 5-Hydroxytryptamine1A Receptor Occupancy in Rats Drug Metab. Dispos., August 1, 2007; 35(8): 1387 - 1392. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Yassen, E. Olofsen, A. Dahan, and M. Danhof Pharmacokinetic-Pharmacodynamic Modeling of the Antinociceptive Effect of Buprenorphine and Fentanyl in Rats: Role of Receptor Equilibration Kinetics J. Pharmacol. Exp. Ther., June 1, 2005; 313(3): 1136 - 1149. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Raghavan, D. Zhang, M. Zhu, J. Zeng, and L. Christopher CYP2D6 CATALYZES 5-HYDROXYLATION OF 1-(2-PYRIMIDINYL)-PIPERAZINE, AN ACTIVE METABOLITE OF SEVERAL PSYCHOACTIVE DRUGS, IN HUMAN LIVER MICROSOMES Drug Metab. Dispos., February 1, 2005; 33(2): 203 - 208. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. P. Zuideveld, P. H. Van der Graaf, D. Newgreen, R. Thurlow, N. Petty, P. Jordan, L. A. Peletier, and M. Danhof Mechanism-Based Pharmacokinetic-Pharmacodynamic Modeling of 5-HT1A Receptor Agonists: Estimation of in Vivo Affinity and Intrinsic Efficacy on Body Temperature in Rats J. Pharmacol. Exp. Ther., March 1, 2004; 308(3): 1012 - 1020. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Nakajima |