Vol. 27, Issue 6, 701-709, June 1999
Cardiovascular Effect and Simultaneous Pharmacokinetic and
Pharmacodynamic Modeling of Pimobendan in Healthy Normal Subjects
Kai-Min
Chu,
Oliver Yoa-Pu
Hu, and
Shyh-Ming
Shieh
Pharmaceutical Research Institute, National Defense
Medical Center (O.Y.-P.H.) and Division of Cardiology, Department of
Medicine, Tri-Service General Hospital, National Defense Medical
Center (K.-M.C., S.-M.S.), Taipei, Taiwan, Republic of China
 |
Abstract |
Pimobendan is a new inotropic agent with vasodilator properties. We
have reported the pharmacokinetics of enantiomers of pimobendan in
healthy humans. The present report focuses on the pharmacodynamic effect of pimobendan and a simultaneous pharmacokinetic-pharmacodynamic modeling. Eight normal healthy volunteers were studied with oral administration of 7.5 mg and i.v. administration of 5 mg of racemic pimobendan. Concentrations of enantiomers of pimobendan were determined by high performance liquid chromatography. Cardiovascular effects of
pimobendan were evaluated by echocardiography. Oral pimobendan significantly reduced 29.0% and 16.5% of the left ventricle
end-systolic dimension (LVESD) and end-diastolic dimension,
respectively. The mean velocity of circumferential fiber shortening,
ejection fraction, and fractional shortening significantly increased
105.9%, 29.8%, and 46% from their baseline values, respectively. The
cardiovascular effects of i.v. pimobendan were similar but to a lesser
extent. Plots of effect versus plasma concentration (Cp)
showed counterclockwise hysteresis loops. A hypothetical effect
compartment was established and incorporated into a sigmoid
Emax model to describe the time courses of Cps
of pimobendan and effects on LVESD. The maximal changes
(Emax) in LVESD would be 2.60 ± 0.51 cm and 2.30 ± 0.13 cm as estimated from plasma data of (+)- and (
)-pimobendan,
respectively. The estimated effect-site concentrations corresponding
with 50% of the maximal effect (Ce50) were 6.54 ± 1.35 and 6.64 ± 1.35 ng/ml for (+)- and (
)-pimobendan,
respectively. A simultaneous pharmacokinetic-pharmacodynamic
model could be established to suppress the hysteresis loop and to
predict the pharmacological effect based on Cp.
 |
Introduction |
Pimobendan
[4,5-dihydro-6-2-(4-methoxyphenyl)-1H-benzimidazole-5-yl-5-methyl-3(2H)-pyridazinone]
(pimo)1 is a new benzimidazole-pyridazinone
derivative that exhibits both positive inotropic and vasodilator
properties (Ruegg et al., 1984
; Von Meel, 1985
; Verdouw et al., 1986
)
by specific type III phosphodiesterase inhibition (Scholz and Meyer,
1986
; Schmitz et al., 1989
) and by directly increasing the calcium
sensitivity of the cardiac myofilaments (Ruegg, 1986
; Fujino et al.,
1988
). Oral and i.v. pimo has been shown to increase cardiac
output (CO), heart rate (HR), ejection fraction (EF), left ventricular
maximum rate of change in pressure over time (dP/dt), and
decrease systemic vascular resistance, left ventricle filling pressure,
and pulmonary capillary wedge pressure in a dose-dependent manner in
animals and patients with chronic congestive heart failure (Walter et al., 1988
; Renard et al., 1988
; Hagemeijer et al., 1989a
,b
; Baumann et
al., 1989
; Remme et al., 1989
). However, no hemodynamic data are
available from healthy volunteers.
We have reported the pharmacokinetics of enantiomers of pimo in
healthy humans (Chu et al., 1995a
). The present report focuses on the
pharmacodynamic (PD) effect of pimo. In addition, a simultaneous pharmacokinetic (PK)-PD modeling was established to describe and predict the hemodynamic effect of pimo after oral and i.v. administration.
Materials and Methods
Procedures.
The study protocol was approved by the Human Subjects Institutional
Review Board of the National Defense Medical Center and the Department
of Health, Executive Yuan, Republic of China. The criteria for
the selection of the eight volunteers, the laboratory tests, and the
general procedure have been described previously. Briefly, three 2.5-mg
(7.5 mg) capsules of pimo (Boehringer Ingelheim GmbH, Ingelheim,
Germany) were administered in the oral study. Two weeks after
the oral study, 5 mg of pimo was given i.v. Blood samples were analyzed
for enantiomers of pimo by a coupled achiral-chiral normal-phase
high-performance liquid chromatography (Chu et al., 1992
). Data were
fitted iteratively reweighted as 1/(predicted value) (Sheiner and Beal,
1985
; Sheiner, 1985
) to two- or three-compartment models by using a
nonlinear regression program PCNONLIN (Weiner, 1986
). A lag time was
included in each model in the oral study. In all subjects the
2-compartment model gave the better fit to the data as assessed by the
Akaike information criterion (Yamaoka et al., 1978
). Blood pressure
(BP) and HR were checked and the hemodynamic effects of pimo were
evaluated by echocardiography at 0.5, 1, 1.5, 2, 3, 4, 5, 6, 8, and
10 h after drug administration. A Hewlett-Packard Ultrasound
Imaging System-1500 (Hewlett-Packard Co., Palo Alto, CA) was used for
echocardiographic recording.
Echocardiography Recording and Analysis.
Cardiac function indices including cardiac index (CI), EF, mean
velocity of circumferential fiber shortening (MVcf), systolic time
intervals, and peak flow velocities were evaluated by a Hewlett-Packard SONOS ultrasound system with phased-array Doppler (model 1500). M-mode,
two-dimensional echocardiographic images, and pulsed-wave Doppler
velocity signals were recorded on a Panasonic standard VHS video
cassette recorder (model AG7300) for a later playback and analysis. The
subjects were examined in a slight left-lateral decubitus position
during quiet respiration. Left parasternal long-axis, left parasternal
short-axis, apical 4-chamber, and apical 5-chamber views were recorded.
Mitral valve-flow peak velocity (MVPV) and aortic flow peak velocity
(AOPV) were recorded using a pulsed-wave Doppler. Transmitral flow
velocity recordings were obtained from an apical 4-chamber or apical
2-chamber view using a 2.5 MHz imaging transducer. The sample-volume
cylinder, approximately 1 cm in axial length, was superimposed on the
two-dimensional image at the level of the mitral valve orifice
and oriented parallel to the apparent long-axis of blood flow. The
transducer was then angled to record maximal transmitral flow
velocities. Doppler tracings of aortic flow were obtained using the
apical 5-chamber view, with the sample-volume placed just beyond the
valve leaflets within the proximal aortic root. Slight adjustments in
transducer angulation or sample-volume position were at times required
to maximize the audio and graphic quality of the Doppler signal. The
aortic valve motion on parasternal long-axis view and velocities of
mitral and aortic flow were recorded over several cardiac cycles at a
monitor sweep speed of 100 mm/s. M-mode echocardiographic measurements followed the standard method as recommended by the American Society of Echocardiography (Sahn et al., 1978
). Left ventricle end-diastolic dimension (LVEDD) and left ventricle
end-systolic dimension (LVESD) were measured in parasternal long-axis
view at the level between the papillary muscle and mitral leaflet tips. Left atrium (LA) was measured at the aortic valve level. Left ventricle
volume (V) was determined by the dimension (D) using Teichholz's equation (Teichholz et al., 1976
) V = [7.0/(2.4 + D)] × D3. Stroke volume (SV) equals LVEDV
LVESV; EF equals
SV/LVEDV and CI equals SV × HR/body surface area. Pre-ejection
period (PEP) was the time duration measured from the onset of
ventricular depolarization (onset of QRS complex on simultaneous
electrocardiogram) to the opening of the aortic valve. Left ventricular
ejection time (LVET) was measured from the opening to the closure of
the aortic valve. All measurements were made with the aid of an
off-line computerized analysis station equipped with internal calipers
and a programmable graphic analyzer (Dextra 300, Dextra Medical Inc.,
Long Beach, CA). The recorded images were played back through a
videocassette system equipped with a frame-by-frame bidirectional
search module (Panasonic model AG7300).
Reproducibility and Validation.
The variabilities of echocardiographic indices were evaluated in eight
other normal subjects before conducting this study. Four levels of
coefficient of variation (CV) were evaluated. Level-1 CV was determined
by repeatedly measuring the cardiac function indices in one cardiac
cycle 5 times to evaluate the error resulting from using the pointing
device of a computer digitizer. Level-2 CV was evaluated by measuring
the indices 5 times in five consecutive cardiac cycles to detect the
error caused by sonographer, interpreter, and examined subject in
addition to level-1 error. Level-3 CV was evaluated by examining each
subject at four time points (9:00 AM, 11:00 AM, 2:00 PM, and 4:00 PM)
in a day to identify the time-period variation in addition to the error
of the previous two levels. Level-4 CV detected the interobserver
variation by measurement of all indices by four experienced
echocardiographers. The percentage of coefficient of variation (CV%)
of levels-1 and -2 were less than 5% except in PEP (8.0% in level-2).
Level-3 CV% was less than 10% except in MVPV (11.1%) and CO
(12.0%). The interobserver CV% was less than 10% except in measuring
LV posterior wall dimension (10.6%).
In this study, all indices were the average of at least three
measurements in different cardiac cycles. An observer who was blinded as to PK data performed all measurements.
Simultaneous PK-PD Modeling.
Plasma and red blood cell concentrations of (+)- and (
)-pimo and
changes in LVESD were used to characterize the simultaneous PK-PD
model. Plots of (+)-pimo concentration versus percent changes in LVESD
in time sequence showed counterclockwise hysteresis loops. Therefore,
two different levels of effect may be seen at a single plasma
concentration (Cp) and two different
Cps of pimo may produce the same effect, which
seems against PD principles. However, this might happen in conditions
in which the site of action of pimo is kinetically distinguishable from
the plasma compartment. To collapse the hysteresis loop and determine
the relationship between Cp and the effect, one
may: 1) sample effect-site concentration (Ce), 2)
perform multiple steady-state experiments, or 3) model the effect site
as a kinetic "compartment" linked to the PK compartments. We
postulated a hypothetical effect compartment linked to the plasma
compartment by a first order process (link PK model; Fig. 1). A two-compartment model of pimo
disposition was used when modeling the PK data. It was assumed that
drug enters and leaves the effect compartment by a first order process
with no appreciable reflux of drug back into the PK system. In this
link PK model, the concentration and time course of drug at the effect
compartment are determined by the elimination rate constant of the
effect-compartment (ke0). The greater the
ke0, the less is the time lag between the central
and effect compartments. To determine an appropriate
ke0, a parametric approach with sigmoid maximal
effect (Emax) PD model was included in the
simultaneous PK-PD model.

View larger version (9K):
[in this window]
[in a new window]
|
Fig. 1.
This diagram represents the link PK model
defined as the link of the hypothetical compartment to the PK model.
By assuming that the transfer of the mass of drug from the central
compartment to the effect compartment is relatively small, the drug
moving back to the central compartment is negligible and will not add
another compartment to the original two-compartment PK model. GI:
gastrointestinal compartment; C: central compartment; P: peripheral
compartment; E: effect compartment.
|
|
To establish a simultaneous PK-PD model, the first step is to obtain PK
parameters. In all subjects the 2-compartment PK model gave the best
fit for the data as assessed by the Akaike information criterion
(Yamaoka et al., 1978
).
The PK model: (see Appendix eq. 1)
where Cp represents plasma
concentration, t is time, ka is the
first-order absorption rate constant,
1
and
2 are the disposition rate
constants, and
,
, and
represent hybrid constants.
The second step is to apply the PK parameter in the link PK model
(Sheiner et al., 1979
). In the link PK model, the following equation
was used to describe the time course of the concentration in the
effect-compartment (Ce):
The link PK model: (see Appendix eq. 3)
The third step is to integrate the link PK model into a PD
model, which becomes a simultaneous PK-PD model (Schuttler et al.,
1987
).
The simultaneous PK-PD model:
In this equation, E is effect, E0 is the
calculated baseline effect, Emax is the
maximum change in predicted response that can be produced by the drug.
Ce50 is the concentration at the effect
site causing 50% of Emax and n is the Hill
coefficient (or slope factor), which determines the sigmoidicity of the
concentration-effect curve.
Using a nonlinear curve fitting program PCNONLIN, the parameters:
E0, Emax, n,
Ce50, and ke0 were
obtained to best fit the model to the data of effect and time.
For i.v. administration the third exponential term in the PK and the
link PK models were omitted.
Statistical Analysis.
The changes in the PD data according to time were evaluated by repeated
measured analysis of variances followed by Dunnett's multiple
comparisons. A paired t test was used to compare the parameters of the simultaneous PK-PD model in oral and i.v. studies. Differences in the p values of < .05 (two-tailed) were
considered significant. Results were reported as mean ± S.E..
 |
Results |
Pharmacodynamics.
After oral administration of 7.5 mg of pimo, the dimension of the heart
decreased (Fig. 2 and Table
1). The LVESD and LVEDD significantly
reduced 29.0% and 16.5% of their baseline values at 2 to
5 h, respectively. This effect lasted for at least 10 h. In a
similar trend, the LA dimension, LVET, and PEP were significantly reduced 20 to 30%. The parameters of LV contractile force including MVcf, EF, and fractional shorting were significantly increased to 105.9%, 29.8%, and 46%, respectively. BP showed no significant change, whereas HR increased 38%. The AOPV, which was determined by
the SV, aortic valve area, afterload, and contractility, was also
increased significantly to 29.3% in 2 to 4 h. The stroke volume index (SVI) reduced insignificantly in 2 to 4 h and
returned toward baseline after 5 h. The CI, calculated by SVI × HR, increased gradually but not significantly until 5 h. This
increase in CI resulted from an increase in HR. The hemodynamic change
after i.v. administration of 5 mg of pimo showed a similar change in the parameter but to a lesser extent (Fig.
3). Because the bioavailability of pimo
was about 50% (Chu et al., 1995a
), a greater effect was seen with a
smaller available dose of pimo (3.75 versus 5 mg), which suggested that
the demethylated metabolite formed in first pass metabolism was
contributing to the overall hemodynamic effect.

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 2.
This plot represents the cardiovascular
effects of pimo in eight normal healthy volunteers after oral
administration of 7.5 mg of racemic pimo. SBP and DBP: systolic and
diastolic blood pressure.
Decrease in cardiac chamber dimensions and increase in
ejection-phase indices, HR, and CI developed gradually toward their
maxima in 2 to 5 h and lasted for 8 to 10 h. Data were
expressed in mean ± S.E.. *represents p < .05 by using repeated measured ANOVAs followed by Dunnett's multiple
comparisons.
|
|

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 3.
This plot demonstrates the cardiovascular
effects of pimo in eight normal healthy volunteers after i.v.
administration of 5 mg of racemic pimo.
Legends are the same as in Fig. 1. The changes are similar to
those in the oral study but to a lesser extent. Diastolic BP decreased
significantly in 2 to 6 h. CI did not change significantly.
|
|
The adverse effect was unremarkable. One subject experienced dizziness
while voiding at 1.75 h after the oral dose, whereas others had
lightheadedness or palpitations. There were no arrhythmias found on the
electrocardiogram tracings in echocardiographic recording.
Simultaneous PK-PD Model.
There was a delay between the peaking time of plasma pimo concentration
and the Emax (changes in LVESD; Fig.
4). When the effects were plotted against
Cps in time sequence, large
counterclockwise hysteresis loops were observed for oral and i.v.
administration, respectively (Fig. 5A and
B). After applying the simultaneous PK-PD model, the hypothetical
Ces were obtained along with the estimated
effect (Eest; Fig. 4). It is worth noting that
observed effect (Eobs) is unable to fit or
estimate unless the effect site concentration of drug is known. This
figure demonstrates a good fitting between Eobs
and Eest. When effects were plotted against the
effect-compartment concentrations, the hysteresis loop collapsed (Fig.
5, C and D). This demonstrated that concentration-effect relationship
or pharmacodynamics of pimo could be found by using the simultaneous
PK-PD model.

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 4.
Observed Cp of
(+)-pimo and Eobs and Eest versus time curves
after oral administration of 7.5 mg of racemic pimo in eight normal
healthy volunteers.
A simultaneous PK and PD model was established using observed
Cps and Eobs. There is a good
correlation between Eobs and Eest. The dotted
line represents hypothetical effect compartment concentration of
(+)-pimo. Data are presented as mean ± S.E.
|
|

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 5.
Eobs and Eest versus
Cps of (+)-pimo in plasma plotted in time
sequence in eight normal healthy volunteers after oral (A) and i.v. (B)
administrations of pimo showed counterclockwise hysteresis loops.
The curves of Eest showed good fitting with the
observed data. It was unable to fit the data without knowing the
Ce, which was obtained by applying the
simultaneous PK-PD model. When effect was plotted against the
hypothetical effect compartment concentration (C, D), the hysteresis
loop collapsed and the concentration-effect relationship could be
found.
|
|
Table 2 shows the parameters of the
simultaneous PK-PD model by fitting the effect with concentration of
enantiomers of pimo in plasma and in red blood cells. In the oral
study, the observed and estimated baseline values of LVESD were similar
(3.45 ± 0.11 cm versus 3.36 ± 0.13 cm as for (+)-pimo). The
LVESD may reduce 1.3 cm (Emax/2) when
Cp reaches the
Ce50 or 6.54 ng/ml [for (+)-pimo]. The
time to reach Emax or maximal concentration in
effect-compartment was 2.28 h, which was longer than
Tmax of Cp of pimo.
The "n" was about 2.4. The ke0 was 0.33 h
1. There was no statistical difference between
enantiomers. In the i.v. study, ke0 was
significantly smaller than for those in the oral study, yet other
parameters were similar to those in the oral study.
View this table:
[in this window]
[in a new window]
|
TABLE 2
Parameters of simultaneous PK and PD modeling of pimo (enantiomer
concentration versus LVESD) in eight normal healthy volunteers
|
|
 |
Discussion |
Pharmacodynamics.
We have reported the hemodynamic effects of pimo in patients with
congestive heart failure (Chu et al., 1995b
). In patients, the HR, SV,
and CO were increased and the BP and the filling pressure of the
cardiac chambers were decreased. In our study of normal volunteers, the decrease in chamber size and increase in HR and CO were also significant, but BP and SV were decreased insignificantly. These different effects on SV in patients and in normal subjects might
be due to different preload and inotropic response.
Simultaneous PK-PD Model.
One of the major goals in clinical pharmacology is to have a drug given
in a proper way, proper dose, and suitable duration to achieve and
maintain an ideal therapeutic effect in patients of variable disease
entities. To produce its characteristic effects, a drug must be present
in appropriate concentrations at its sites of action. In most cases,
the concentration of drug in the systemic circulation will be parallel
to that at the sites of action. There is then a direct link between PK
and PD models. For some drugs, there is no clear or simple relationship
between pharmacological effect and concentration in plasma. Many
factors, such as formation of active metabolites, enzyme induction,
enzyme inhibition, receptor desensitization, indirect or irreversible
pharmacological effect, delayed equilibrium in concentration between
plasma and sites of action, and enantiomers with different elimination
rates and potency (Holford and Sheiner, 1981
; Mehvar, 1992
) may account for the lack of correlation. Although the growth of pharmacodynamics is
less than that in pharmacokinetics because difficulties exist in the
quantitative determination of meaningful responses, it is possible to
use mathematical models to describe the effect-concentration relationship. The simultaneous PK-PD model is first proposed by Segre
(1968)
and later by Galeazzi et al. (1976)
, Hull et al. (1978)
, and
subsequently elaborated by Sheiner et al. (1979)
. These modeling
techniques have been applied to various drugs including neuromuscular
blocking agents (Evans et al., 1984
), anesthetics (Stanski et
al., 1984
), bronchodilators (Oosterhuis et al., 1984
), antihypertensives (Reid and Meredith, 1990
), benzodiazepines (Buhrer et
al., 1990
), antibiotics (Garrison et al., 1990
), steroids (Kong et al.,
1989
), and diuretics (Hammarlund et al., 1985
). We herein applied this
integrated model to a cardiotonic agent.
To quantitatively evaluate cardiac performance is a complicated
process. There are a number of methods to evaluate cardiac performance,
such as cardiac catheterization, radionuclide angiocardiography, echocardiography, etc., using a number of parameters based on pressures, flows, volumes, and dimensions. Echocardiography was used
because it is reliable, noninvasive, and relatively inexpensive for the
monitoring of the cardiovascular effect of a drug without neurohumoral
disturbances. However, there are a number of limitations. First, it is
difficult to clarify whether the increased cardiac performance was due
to either an inotropic effect or vasodilatation effect alone or both
without direct measurement of intracardiac pressure and volume changes.
Second, an echocardiography examination takes 5 to 10 min, or even
longer in some cases with a small "cardiac window", to obtain a
series of images one at a time. To monitor the rapidly changing effect,
as in the initial phase of some i.v. studies, only limited numbers of
views can be obtained. Third, calculations of the 3-dimensional volume
or CO from a linear dimension primary measurement will cube the error
and subsequently affect the model fitting. We therefore used the
changes in LVESD as a meaningful PD response to establish a combined
PK-PD model.
The slope factor "n" was about 2.4 in the oral study and 4 to 6 in
the i.v. study. Theoretically, an n value greater than one indicates
positive cooperativity in receptor binding. This means that the
occupancy of one binding site by a ligand enhances the likelihood that
other couple sites on the molecule will preferentially bind the same
ligand. However, interpreting this n value from a entire body point of
view, we would consider this is an overall systemic effect integrated
from the effects on different organ systems. The overall changes in
LVESD might result from an increased inotropic effect, a vasodilating
effect, and an increase in HR directly or indirectly. The more the
organs respond to the drug, the greater the n value could be.
Nevertheless, this cannot explain the higher n value in the i.v. study.
Examination of standard errors showed the greatest uncertainty in the
estimates of n. This might be due to the estimated
Ce, which was smaller than Ce50 and fell in the lower left half of the
sigmoidal effect-log concentration curve.
The ke0 obtained in the i.v. study is
generally smaller than that calculated in an oral study
(p < .05 in (
)-pimo). There are several
possibilities. First, it may be a statistical type-one error. Second,
this discrepancy was caused by overestimated disposition rate constants
and
in the i.v. study because fewer data points were measurable
in the elimination phase (Chu et al., 1995a
). In the link PK model, the
higher
and
value will lead to a lower
ke0 and hence
ke0 will be underestimated. Third, the
amount and peaking time of the active metabolite may affect the
overall effect and result in a relatively delayed peaking of changes in LVESD in the i.v. study. In isolated canine ventricular muscle, Endoh
et al. (1991)
found that the demethylated metabolite of pimo is
one-third the efficacy and 500 times the potency as compared with pimo.
However, to apply these PD parameters in humans one should be cautious
that a metabolite is usually hydrophilic and thus has a different
volume of distribution and tissue affinity. The ratio of a parent drug
to its metabolite in plasma may be different from that in the effect
site. Thus, the in vivo efficacy and potency may be different from
those of an in vitro study. Therefore, to determine the percentage of
contribution of pimo and its metabolite to the overall effect based on
their Cp is difficult. In our study, it is
worth noticing that the time to reach Emax is
over 3 h after i.v. administration, which seems longer than that
in oral administration. In a normal situation, the effect is usually
seen earlier in i.v. administration than oral administration with the
exception that the drug acts more like a prodrug requiring
biotransformation in the liver to become active. This suggests that the
metabolite of pimo contributes to the overall effect. However, it is
yet to be clarified by direct administration of the demethylated
metabolite in humans.
The effect compartment concentration has been compared with a
steady-state concentration. Schwartz et al. (1989)
studied verapamil in
22 normal subjects in a single dose and steady-state infusions to
compare the estimated and real effect site concentration. They found
that use of the PD model to estimate effect site concentration provided
a closer estimate of the true steady-state concentration than the
estimation from the postinfusion Cp.
However, the model was limited in describing higher concentration
versus effect relationships.
In summary, this study demonstrated the cardiovascular effect of
pimo in normal subjects after oral and i.v. dosings. A simultaneous PK-PD model was developed to suppress the hysteresis loop and to
predict the pharmacological effect based on
Cp. It can provide valuable information on
dose-effect relationships and on a choice of optimal dosing intervals
in drug development. It is possible that in patients with variable
disease entities and impaired hepatic and renal functions the changes
in the time to reach peak effect and the therapeutic duration might be
predicted based on different PK parameters, yet this remains unproved.
 |
Footnotes |
Received September 11, 1998; accepted February 10, 1999.
This work was partially supported by the Institute of
Biomedical Sciences, Academia Sinica, Republic of China.
Send reprint requests to: Dr. Oliver Yoa-Pu Hu,
Professor and Director, Pharmaceutical Research Institute, National
Defense Medical Center, P.O. Box 90048-508, Taipei, Taiwan, Republic
of China. E-mail: hyp{at}ndmc1.ndmtsgh.edu.tw
 |
Abbreviations |
Abbreviations used are:
pimo, pimobendan;
AOPV, aortic flow peak velocity;
BP, blood pressure;
Ce, effect-site concentration;
Ce50, effect-site concentration
corresponding with 50% of the maximal effect;
CI, cardiac index;
CO, cardiac output;
Cp, plasma concentration;
CV, coefficient of variation;
dP/dt, rate of change in
pressure over time;
EF, ejection fraction;
Emax, maximal
effect;
E0, calculated baseline effect;
E0(obs), observed baseline effect;
Eest, estimated effect;
Eobs, observed effect;
HR, heart rate;
ke0, elimination rate constant of the
hypothetical effect compartment;
LA, left atrium;
LVEDD, left ventricle
end-diastolic dimension;
LVESD, left ventricle end-systolic dimension;
LVET, left ventricular ejection time;
MVcf, mean velocity of
circumferential fiber shortening;
MVPV, mitral valve-flow peak
velocity;
PD, pharmacodynamic;
PEP, pre-ejection period;
PK, pharmacokinetic;
SV, stroke volume;
SVI, stroke volume index;
Tmax, lag-time between Tmax(e) and
Tmax;
Tmax, time to reach peak concentration;
Tmax(e), time to reach peak effect.
 |
Appendix |
Laplace transform and anti-Laplace transform of
2-compartment model with first order input and the link PK model (Fig.
1)
List of symbols:
F: bioavailability
D: dose
ins,c: Laplace transform for the input
function of central compartment
ds,c: Laplace transform for the disposition
function of central compartment
as,c: Laplace transform for the amount of drug in
the central compartment
ins,e: Laplace transform for the input function
of effect site compartment
ds,e: Laplace transform for the disposition
function of effect site compartment
as,e: Laplace transform for the amount of drug in
the effect site compartment
Xc: amount of drug in the central compartment
Xe: amount of drug in the effect site compartment
Cp: plasma concentration
Ce: effect site concentration
Vc: volume of distribution of central
compartment
Ve: volume of distribution of effect site
compartment
s: Laplace operator
ka, k21,
1,
2: rate constants
Laplace transform for 2-compartment model with first order input:
Anti-Laplace transform:
eq. 1:
Laplace transform for hypothetical effect compartment (so called
link PK model, assuming the effect site receiving negligible actual mass of drug and not affecting 2-compartment model):
Anti-Laplace transform:
eq. 2:
eq. 2 can be rearranged into:
When t = 0, (L + M + N + O) = 0, then eq. 2 becomes:
Combining
eq. 1 and 2, the effect site concentration is:
eq. 3:
 |
References |
-
Baumann G,
Ningel K and
Permanetter B
(1989)
Cardiovascular profile of UDCG 115 BS-pimobendane and reversibility of catecholamine subsensitivity in severe congestive heart failure secondary to idiopathic dilated cardiomyopathy.
J Cardiovasc Pharmacol
13:
730-738[Medline].
-
Buhrer M,
Maitre PO,
Crevoisier C and
Stanski DR
(1990)
Electroencephalographic effects of benzodiazepines. II. Pharmacodynamic modeling of the electroencephalographic effects of midazolam and diazepam.
Clin Pharmacol Ther
48:
555-567[Medline].
-
Chu KM,
Shieh SM,
Wu SH and
Hu OYP
(1992)
Enantiomeric separation of a cardiotonic agent pimobendan and its major active metabolite, UD-CG 212 BS, by coupled achiral chiral normal-phase high-performance liquid chromatography.
J Chromatogr Sci
30:
171-176.
-
Chu KM,
Shieh SM,
Wu Sh and
Hu OYP
(1995a)
Plasma and red blood cell pharmacokinetics of pimobendan enantiomers in healthy Chinese.
Eur J Clin Pharmacol
47:
537-542[Medline].
-
Chu KM,
Shieh SM,
Wu SH and
Hu OYP
(1995b)
Pharmacokinetics and pharmacodynamics of enantiomers of pimobendan in patients with dilated cardiomyopathy and congestive heart failure after single and repeated oral dosing.
Clin Pharmacol Ther
57:
610-621[Medline].
-
Endoh M,
Shibasaki T,
Satoh H,
Norota I and
Ishihata A
(1991)
Different mechanisms involved in the positive inotropic effects of benzimidazole derivative UD-CG 115 BS (pimobendan) and its demethylated metabolite UD-CG 212 Cl in canine ventricular myocardium.
J Cardiovasc Pharmacol
17:
365-375[Medline].
-
Evans MA,
Shanks CA,
Brown KF and
Triggs EJ
(1984)
Pharmacokinetic and pharmacodynamic modelling with pancuronium.
Eur J Clin Pharmacol
26:
243-250[Medline].
-
Fujino K,
Sperelakis and
Solaro RJ
(1988)
Sensitization of dog and guinea pig heart myofilaments to Ca2+ activation and the inotropic effect of pimobendan: comparison with milrinone.
Circ Res
63:
911-922[Abstract/Free Full Text].
-
Galeazzi RL,
Benet LZ and
Sheiner LB
(1976)
Relationship between pharmacokinetics and pharmacodynamics of procainamide.
Clin Pharmacol Ther
20:
278-289[Medline].
-
Garrison MW,
Vance-Bryan K,
Larson TA,
Toscano JP and
Rotschafer JC
(1990)
Assessment of effects of protein binding on daptomycin and vancomycin killing of Staphylococcus aureus by using an in vitro pharmacodynamic model.
Antimicrob Agents Chemother
34:
1925-1931[Abstract/Free Full Text].
-
Hagemeijer F,
Brand HJ and
Roth W
(1989a)
Cardiovascular effects and plasma level profile of pimobendan (UD-CG 115 BS) and its metabolite UD-CG 212 in patients with congestive heart failure after single and repeated oral dosing.
J Cardiovasc Pharmacol
14:
302-310[Medline].
-
Hagemeijer F,
Brand HJ and
van Mechelen R
(1989b)
Hemodynamic effects of pimobendan given orally in congestive heart failure secondary to ischemic or idiopathic dilated cardiomyopathy.
Am J Cardiol
63:
571-576[Medline].
-
Hammarlund MM,
Odlind B and
Paalzow LK
(1985)
Acute tolerance to furosemide diuresis in humans. Pharmacokinetic-pharmacodynamic modeling.
J Pharmacol Exp Ther
233:
447-453[Abstract/Free Full Text].
-
Holford NH and
Sheiner LB
(1981)
Understanding the dose-effect relationship: clinical application of pharmacokinetic-pharmacodynamic models.
Clin Pharmacokinet
6:
429-453[Medline].
-
Hull CJ,
VanBeem BH,
McLeod K,
Sibbald A and
Watson MJ
(1978)
A pharmacokinetic model for pancuronium.
Br J Anaesth
50:
1113-1123[Abstract/Free Full Text].
-
Kong AN,
Ludwig EA,
Slaughter RL,
DiStefano PM,
DeMasi J,
Middleton E, Jr. and
Jusko WJ
(1989)
Pharmacokinetics and pharmacodynamic modeling of direct suppression effects of methylprednisolone on serum cortisol and blood histamine in human subjects.
Clin Pharmacol Ther
46:
616-628[Medline].
-
Mehvar R
(1992)
Stereochemical considerations in pharmacodynamic modeling of chiral drugs [letter].
J Pharm Sci
81:
199-200[Medline].
-
Oosterhuis B,
Berge RJ,
Sauerwein HP,
Endert E,
Schellekens PT and
Boxtel CJ
(1984)
Pharmacokinetic-pharmacodynamic modeling of prednisolone-induced lymphocytopenia in man.
J Pharmacol Exp Ther
229:
539[Abstract/Free Full Text].
-
Reid JL and
Meredith PA
(1990)
Concentration-effect analysis of antihypertensive drug responses.
Hypertension
16:
12-18.
-
Remme WJ,
Wiesfeld AC,
Look MP and
Kruyssen HA
(1989)
Hemodynamic effects of intravenous pimobendan in patients with left ventricular dysfunction.
J Cardiovasc Pharmacol
14 Suppl 2:
S41-S44.
-
Renard M,
Walter M,
Liebens I,
Dresse A and
Bernard R
(1988)
Pimobendane (UD-CG 115 BS) in chronic congestive heart failure. Short-term and one-month effects of a new inotropic vasodilating agent.
Chest
93:
1159-1164[Abstract/Free Full Text].
-
Ruegg JC
(1986)
Effects of new inotropic agents on Ca2+ sensitivity of contractile proteins.
Circulation
73:
III78-III84.
-
Ruegg JC,
Pfitzer G,
Eubler D and
Zeugner C
(1984)
Effect on contractility of skinned fibres from mammalian heart and smooth muscle by a new benzimidazole derivative, 4,5-dihydro- 6-[2-(4-methoxyphenyl)-1H-benzimidazol-5-yl]-5-methy l-3(2H)- pyridazinone.
Arzneimittel forschung
34:
1736-1738[Medline].
-
Sahn DJ,
DeMaria A,
Kisslo J and
Weyman A
(1978)
Recommendations regarding quantitation in M-mode echocardiography: Results of a survey of echocardiographic measurements.
Circulation
58:
1072-1083[Abstract/Free Full Text].
-
Scholz H and
Meyer W
(1986)
Phosphodiesterase-inhibiting properties of newer inotropic agents.
Circulation
73:
III99-III108.
-
Schmitz W,
von der Leyen H,
Meyer W,
Neumann J and
Scholz H
(1989)
Phosphodiesterase inhibition and positive inotropic effects.
J Cardiovasc Pharmacol
14:
s11-s14.
-
Schuttler J,
Stanski DR,
White PF,
Trevor AJ,
Horai Y,
Verotta D and
Sheiner LB
(1987)
Pharmacodynamic modeling of the EEG effects of ketamine and its enantiomers in man.
J Pharmacokinet Biopharm
15:
241-253[Medline].
-
Schwartz JB,
Verotta D and
Sheiner LB
(1989)
Pharmacodynamic modeling of verapamil effects under steady-state and nonsteady-state conditions.
J Pharmacol Exp Ther
251:
1032-1038[Abstract/Free Full Text].
-
Segre G (1968) Kinetics of interaction between drugs and biological
systems. Farmaco [Sci] 23:907-918.
-
Sheiner LB
(1985)
Analysis of pharmacokinetic data using parametric models. II. Point estimates of an individual's parameters.
J Pharmacokinet Biopharm
13:
515-540[Medline].
-
Sheiner LB and
Beal SL
(1985)
Pharmacokinetic parameter estimates from several least squares procedures: superiority of extended least squares.
J Pharmacokinet Biopharm
13:
185-201[Medline].
-
Sheiner LB,
Stanski DR,
Vozeh S,
Miller RD and
Ham J
(1979)
Simultaneous modeling of pharmacokinetics and pharmacodynamics: application to d-tubocurarine.
Clin Pharmacol Ther
25:
358-371[Medline].
-
Stanski DR,
Hudson RJ,
Homer TD,
Saidman LJ and
Meathe E
(1984)
Pharmacodynamic modeling of thiopental anesthesia.
J Pharmacokinet Biopharm
12:
223-240[Medline].
-
Teichholz LE,
Kreulen T,
Herman MV and
Gorlin R
(1976)
Problems in echocardiographic volume determinations: echocardiographic-angiographic correlation in the presence or absence of asynergy.
Am J Cardiol
37:
7-11[Medline].
-
Verdouw PD,
Hartog JM,
Duncker DJ,
Roth W and
Saxena PR
(1986)
Cardiovascular profile of pimobendan, a benzimidazole-pyridazinone derivative with vasodilating and inotropic properties.
Eur J Pharmacol
126:
21-30[Medline].
-
Von Meel JC
(1985)
Cardiovascular effects of the positive inotropic agents pimobendan and sulmazole in vivo.
Arzneimittelforschung
35:
284-288[Medline].
-
Walter M,
Liebens I,
Goethals H,
Renard M,
Dresse A and
Bernard R
(1988)
Pimobendan (UD-CG 115 BS) in the treatment of severe congestive heart failure. An acute haemodynamic cross-over and double-blind study with two different doses.
Br J Clin Pharmacol
25:
323-329[Medline].
-
Weiner DL
(1986)
NONLIN84/PCNONLIN: software for the statistical analysis of nonlinear models.
Methods Find Exp Clin Pharmacol
8:
625-628[Medline].
-
Yamaoka K,
Nakagawa T and
Uno T
(1978)
Application of Akaike's information criterion (AIC) in the evaluation of linear pharmacokinetic equations.
J Pharmacokinet Biopharm
6:
165-175[Medline].
0090-9556/99/2706-0701-0709$02.00/0
DRUG METABOLISM AND DISPOSITION
Copyright © 1999 by The American Society for Pharmacology and Experimental Therapeutics