(Hypertension. 1995;25:564-569.)
© 1995 American Heart Association, Inc.
Articles |
From the Research Laboratories, Schering-Plough, Comazzo (E.B., E.O., C.C., A.M.), Milan, and the Department of Public Health and Biostatistics, University of Pisa (S.M.), Italy.
Correspondence to Dr Erminio Bonizzoni, Research Laboratories, Schering-Plough SpA, I-20060 Comazzo, Milan, Italy.
| Abstract |
|---|
|
|
|---|
Key Words: telemetry hemodynamics models, statistical hypertension rats, inbred SHR adenosine
| Introduction |
|---|
|
|
|---|
The telemetric monitoring system has the advantage of providing almost continuous BP and HR profiles (eg, every 5 minutes) from conscious, unstressed, and freely moving animals. So far, little effort has been made to develop parametric functions analogous to pharmacokinetic functions3 and suitable for modeling and defining the response profile to drugs by means of a few constants. Current methods for the analysis of telemetric data often waste information, provide a biased picture of the shape of BP and HR time course, and cause confusion in drug comparisons. Moreover, because peak response time (tmax) varies from one animal to another, chronograms reporting the means of the values obtained at each time point for each animal4 5 6 7 tend to flatten the response curves, which may lead to underestimation of the intensity of the drug effects.
To exploit fully the large amount of information provided by the telemetry system, we elaborated a mathematical model that can be useful in cardiovascular pharmacology studies. In this article, we present and discuss the properties of a family of mechanistic nonlinear models suitable for describing the hemodynamic response profile by means of a limited number of location, intensity, and shape constants, some of which have straightforward pharmacological meaning.
To test the reliability of the model, we analyzed the hemodynamic effects induced by adenosine agonists currently under investigation in our laboratories.8 9 These drugs, with different selectivity for either A1 or A2a adenosine receptor subtypes, have the advantage of producing clearly differentiated hemodynamic effects. After acute administration in conscious rats, A1 selective agonists, such as 2-chloro-N6-cyclopentyladenosine (CCPA), elicit a decrease in BP accompanied by bradycardia,10 whereas A2a selective agonists, such as 2-hexynyl-5'-N-ethylcarboxamidoadenosine (2-hexynyl-NECA), induce hypotension with reflex tachycardia.8 11 All these effects are characterized by a single peak of activity after drug administration. These models also may be applied to drugs that produce two opposite peaks of activity over time, as in the case of the nonselective adenosine agonist 5'-N-ethylcarboxamidoadenosine (NECA), which induces bradycardia followed by tachycardia.
| Methods |
|---|
|
|
|---|
Surgery
Rats were anesthetized with an intraperitoneal administration of
pentobarbital (30 mg/kg). A midline laparotomy was performed, and a
segment (about 1 cm below the renal arteries) of the abdominal aorta
was isolated. The catheter tip was inserted into the aorta just above
the iliac bifurcation, the transmitter was fixed to the muscles, and
the abdominal wall was sutured. After recovery from anesthesia, the
rats were housed individually in cages placed on the receivers.
Gentamicin (5 mg/kg twice daily) was administered for 2 days after
surgery to prevent infection.
Hemodynamic Measurements
BP and HR were recorded and analyzed by a Dataquest 4 system
(Data Sciences Intl). The system consists of BP sensors (model
TA11PAC40), receivers (model RA1010), a consolidation matrix (BCM100),
a personal computer (Everex 486/33), and Dataquest 4 software.
Hemodynamic recordings were taken every 5 minutes, starting 2 hours before drug administration and for 24 hours thereafter. Each recording lasted 10 seconds, and the hemodynamic values of all cardiac cycles within this period (about 50 at baseline) were averaged.
Drugs
The 2-hexynyl-NECA was synthesized at the University of
Camerino, Italy (Prof G. Cristalli). CCPA and NECA were purchased from
Research Biochemical Inc. Drugs were administered intraperitoneally,
dissolved in 2% dimethyl sulfoxide. The 2-hexynyl-NECA and CCPA were
tested at two dose levels to determine the sensitivity of the model for
detection of hemodynamic effects. Only one dose was tested for NECA to
illustrate the application of the biphasic function.
Statistical Evaluation
Rationale for and Characteristics of the Model
In control SHR, circadian rhythms of hemodynamic profiles were
found to be more prominent for HR (minimum, 280 beats per minute
[bpm]; maximum, 350 bpm, when lights were turned on and off,
respectively) than for BP (where differences ranged within 10 mm Hg).
Furthermore, "intrarat" random fluctuations (ie, variability
unexplained by circadian rhythms) were wide (SD, 12 mm Hg and 30 bpm
for BP and HR, respectively), as shown in Fig 1. This
"physiological noise" adds up to the effects caused by the
administration of an antihypertensive drug and sometimes tends to mask
them. In our experience, when the effect is much greater than noise,
the typical pattern of the response curve displays a single asymmetric
peak, with a prompt onset followed by a slow decrease toward the
baseline.
|
The above pattern may be modeled by the following family of four-constant exponential functions (model 1):
![]() |
where E[y(t)] is the expected value of the effect y(t)
recorded at time t since drug administration, g( ) is any monotonic
function of t,
is the maximum intensity of the effect (peak), and
is the time to peak. The shape of the curve depends on the function
g( ) and constant
(
>0), whereas ß (ß>0) expresses the
width of the peak: for a given g( ) and
, the larger the ß, the
narrower the peak.
A suitable choice for g( ) is the Box-Cox transformation
![]() |
restricting attention to
<0 or
=0; in the latter case,
g(t+
)=log(t+
). For a given g( ), increasing values of
reduce
the asymmetry of the peak but may lead to biased estimates of the
response, which is supposed to be null at t=0. Conversely, if the shape
constant
is nearly zero, the expected effect E[y(t)] tends toward
zero when time t approaches zero. In any case, E[y(t)] tends toward
zero when t
. In this article, we considered only the case when
=0, ie, the log transformation. Lower values of
(eg,
=-1,
the reciprocal transformation) may apply to peaks with more pronounced
asymmetry.
Peak value ymax (
) and time to peak
tmax (
) estimated as constants of the function
and other constants that may be derived from the model, such as area
under the curve and times required to observe a 50%
(t1/2)
or a 75%
(t1/4)
reduction of maximum effect, can be directly compared with the
homologous constants typical of pharmaco-kinetics.
In this context, however, it should be noted that kinetics of reversible pharmacological effects and kinetics of plasma levels of drugs are linked by complex nonlinear relations,12 so that neither shape of response curves nor kinetic constants are expected to overlap.
The adequacy of the above models to describe the response y(t) to the administration of a drug depends critically on what we consider to be a response. Under the assumption that physiological fluctuations and drug-induced effects are simply additive, the response should correspond to the hemodynamic value yD(t) recorded at time t after drug administration minus the value y0(t) expected at the same time t after vehicle administration. In this case, yD(t) may be estimated by fitting an appropriate model to the 24-hour hemodynamic profile recorded after vehicle. When the above assumption is not tenable or when circadian rhythms can be considered negligible compared with drug effects (throughout the 24-hour recording or at least in the period of maximum effect), y0(t) may be regarded as constant over t and set equal to the mean values recorded during the day before treatment (or during only a few hours before treatment). These approximations are reasonably accurate if drug administration and most significant responses occur in the morning when average HR and BP values undergo small, spontaneous variations and the recorded peak may be ascribed mainly or entirely to drug effects. Furthermore, because of their rigid structure, the above models are rather insensitive to fluctuations far from the peak.
A more complex eight-constant model (model 2) derived from the juxtaposition of two components (each defined as in model 1) may apply to drugs that induce bradycardia followed by tachycardia, such as the nonselective adenosine receptor agonist NECA:
![]() |
![]() |
where the constants
1,
ß1,
1, and
1
and
2, ß2,
2, and
2 (of the first and second
peaks, respectively) have the same meaning as in model 1, provided that
the second effect occurs after the end of the first without
overlapping.
Fitting Procedure and Analysis of Residuals
Least-squares estimates of the constants of the above models
were obtained by means of PROC NLIN,13 which
resorts to the iterative Gauss-Newton algorithm as modified by
Marquardt.14 This procedure allows the definition of
bounds within which the estimate must be restrained; eg, the estimates
of peak width (ß) and time to peak (
) should be positive, and
hypotensive effect (
) should be negative. Furthermore, the procedure
allows the definition of a grid of initial estimates for the constants,
from which the best combination (the combination that leads to the
lowest residual sum of squares) is selected as the starting point for
the iterative procedure. Bounds and grid are particularly useful when
low doses are tested and hemodynamic effects are weak.
The mean response curve for each drug dose was expressed as the "mean constant curve,"15 ie, a curve with constants that are the mean of the individual kinetic constants estimated in each rat. This curve indicates the mean location and height of the peak of the effect and the mean pattern of drug response. Mean constant curves are commonly used in pharmacokinetics and auxology because different subjects may have similar kinetics as regards serum concentration of a drug and growth velocity of a somatic trait but different locations of maxima and minima of the kinetic function.
Clearly, any model gives a rather approximate representation of the "true" shape of the kinetics of the response; therefore, fitted values can show, at fixed times, a systematic tendency to underestimate or overestimate the measured values. To detect these tendencies, we calculated for each rat the differences (the so-called residuals) between observed and fitted values.16 The average of residuals for all rats at a given time estimates the extent of the systematic departure of the model from the true kinetics of the response at that time. The plot of mean residuals versus time is a valuable tool to compare the adequacy of different models and to describe the pattern of circadian rhythms after drug administration.
| Results |
|---|
|
|
|---|
Fig 2 shows the kinetics of the hemodynamic effects induced by the adenosine A1 agonist CCPA (0.03 and 0.3 mg/kg IP). CCPA induced hypotension and bradycardia, which are dose-related for both peak effect and half-life (t1/2). Regarding time to peak (tmax), the dose-related effect is observed only for HR (Table 1).
|
|
The function models systolic BP (SBP) profiles adequately until almost complete recovery to baseline (Fig 2a and 2b). Systematic departures, ranging within 15 mm Hg, appear after 2 hours at the lower dose and after 6 hours at the higher dose.
For HR, fit is satisfactory during the first 6 to 8 hours after CCPA administration (Fig 2c and 2d). When lights were turned off, HR increased and the function, which models drug effects and not circadian rhythms, was constantly lower than raw HR profiles. In particular, at 0.3 mg/kg, HR values were higher by 25 to 80 bpm (Fig 2d).
Fig 3 shows the kinetics of hypotensive and tachycardic effects induced by the adenosine A2a agonist 2-hexynyl-NECA (0.003 and 0.03 mg/kg IP). Peak effect and half-life are dose-related, whereas time to peak is unaffected by dose.
|
At the lower dose, SBP fitting is satisfactory during the 24 hours (Fig 3a); at the higher dose, the analysis of residuals reveals that in the neighborhood of t1/2, the recorded hypotensive effects decreased more quickly than the function, even though differences were within 10 mm Hg. Subsequently, wide pressure fluctuations were apparent (Fig 3b).
As regards the HR fitting, at the lower dose the function was close to the HR profile for 1.5 hours until recovery to baseline. The following systematic departures, caused by circadian rhythms, were not modeled by the function (Fig 3c). At the higher dose, the function fits the HR profile adequately up to 3.5 hours after drug administration but overestimates the recorded values by 10 to 30 bpm between 4 and 8 hours (Fig 3d). At both doses, a sharp and transitory peak followed by a fluctuating pattern of HR occurred when lights were switched off. These phenomena, which seem to be unrelated to drug administration, did not affect the model.
Like the other adenosine agonists, NECA (0.1 mg/kg IP) induced hypotension. The SBP decrease showed a delayed time to peak (tmax=2.2 hours) and was longer lasting than that observed with the other compounds (Table 2).
|
Fig 4 gives the kinetics of the effects induced on HR by NECA (0.1 mg/kg). The initial A1-mediated bradycardia and the following prolonged period of reflex tachycardia were modeled by a double-exponential function (model 2). Peaks occur at 1.26 and 13.86 hours after drug administration (Table 2), and the fit is quite satisfactory during the whole 24-hour period (Fig 4).
|
| Discussion |
|---|
|
|
|---|
So far, studies based on the telemetry technique usually provided a qualitative description of the changes in hemodynamic parameters induced by drugs, focusing on recorded maximum effect (peak) and time to peak (tmax). Fourier analysis was recently applied as the fitting model, but this approach can be used only to evaluate the alterations of circadian rhythms.17 In the absence of a model, the peak value is often estimated by the mean response values recorded at the mean tmax.18 19 Because tmax is different for each subject in the same experiment, this approach leads to underestimation of the maximum effect. On the other hand, if the peak value is estimated by the mean of each peak effect,4 5 the estimate tends to be positively biased. In fact, the response is likely to show its apparent maximum when random fluctuation and drug effect have the same algebraic sign. Furthermore, in the absence of a model, the time required to observe a 50% reduction of peak effect is difficult to identify, especially when random fluctuations are large; the use of raw BP or HR values20 or a least-squares fit of values in the neighborhood of the peak21 may not be the best choice. For example, the use of raw HR values observed after administration of 2-hexynyl-NECA (0.003 mg/kg IP) led to t1/2=1.30 hours (0.60 to 2.01); this value is higher and has a wider confidence interval than the estimate based on our model (Table 1). According to statistical inference, the wider the confidence interval, the less reliable the estimate.
In most studies, different hemodynamic profiles have been compared at each time point. Such an approach involves a huge number of statistical tests, leading to high risk of type I or II errors and contradictory results (ie, close sequences of significant and nonsignificant differences when observations contiguous on the time axis are considered22 ). Conversely, with our approach, estimates of peak effect and tmax, obtained by fitting a rather rigid mathematical function, are based on the whole response profile and not on one single observation and thus are expected to be both unbiased and precise. With this model, other parameters, related either to persistence (t1/2 or t1/4) or amount (the area under the response curve at fixed times) of the effects, can be easily estimated. In addition, with our model, which summarizes the response curve in a limited number of constants, the effects of drugs can be compared using appropriate univariate or multivariate tests covering hemodynamic profile, peak effects, time to peak, and time to recovery.
After elaboration of mathematical functions, models of hemodynamic profiles were made with three drugs that have different selectivities for the A1 or A2a adenosine receptor. These drugs are good pharmacological tools for their different effects on the hemodynamic parameters examined. The A1 agonist CCPA markedly reduces both BP and HR10 ; the A2a agonist 2-hexynyl-NECA decreases BP and induces reflex tachycardia.11 The nonselective agonist NECA produces hypotension and biphasic effects on HR, ie, tachycardia at lower doses and initial bradycardia followed by tachycardia at higher doses.8 23 24
In the present study, the administration of the above drugs to SHR confirmed the expected changes of the hemodynamic profile. Elaboration of the data on an IBM 486 SX personal computer allowed us to calculate not only BP and HR profiles but also the hemodynamic kinetic constants for each drug, thus providing a quantitative estimate of maximum effect, its occurrence, and its persistence (t1/2). Furthermore, a dose relation of the drug effects was detected. Biphasic effects on HR produced by NECA might be a consequence of the different stimulation of A1 and A2a receptors achieved at different plasma concentrations.
We developed this mathematical function to provide cardiovascular pharmacologists with a sound statistical tool for the analysis of continuous recordings of hemodynamic parameters obtained by telemetry in laboratory animals.
In this study, we applied mathematical functions only to BP and HR profiles, but it may be useful to evaluate other cardiovascular parameters (eg, arterial blood flow or resistance). We administered compounds by the intraperitoneal route; however, this analysis may be generalized to other administration routes. Finally, this model may be implemented with commercially available programs for nonlinear curve fitting, such as PROC NLIN (SAS), SPSS, ALLFIT, and PHARMFIT. The technical background is the same as that required for fitting pharmacokinetic models.
For investigation of antihypertensive agents and, more generally, hemodynamic effects of new drugs, an in vivo screening is useful for the assessment of their duration of action, efficacy, and potency. In this context, the impact of the telemetry system is arousing interest2 and will surely be used more extensively in the near future. The mathematical models we have described aim to exploit the large amount of information provided by the telemetry technique and are therefore expected to be useful in characterizing the hemodynamic profiles of drugs of interest in cardiovascular pharmacology.
Received July 22, 1994; first decision September 8, 1994; accepted December 5, 1994.
| References |
|---|
|
|
|---|
2. Bazil MK, Krulan C, Webb RL. Telemetric monitoring of cardiovascular parameters in conscious spontaneously hypertensive rats. J Cardiovasc Pharmacol. 1993;22:897-905. [Medline] [Order article via Infotrieve]
3.
Mathot RAA, Van Schaick EA, Langemeijer MWE, Soudijn W,
Breimer DD, Ijzerman AP, Danhof M. Pharmacokinetic-pharmacodynamic
relationship of the cardiovascular effects of adenosine A1
receptor agonist N6-cyclopentyl-adenosine in the
rat.J Pharmacol Exp Ther. 1993;268:616-624.
4. Barres C, Cerutti C, Morin B, Sassard J. Cardiovascular effects of two new calcium antagonists, PY 108-068 and PN 200-110, in conscious spontaneously hypertensive rats. Br J Pharmacol. 1988;93:176-184. [Medline] [Order article via Infotrieve]
5. Burges RA, Dodd MG, Gardiner DG. Pharmacologic profile of amlodipine. Am J Cardiol. 1989;64:10I-20I. [Medline] [Order article via Infotrieve]
6. Dodd MG, Gardiner DG, Carter AG, Sutton MR, Burges RA. The hemodynamic properties of amlodipine in anesthetized and conscious dogs: comparison with nitrendipine and influence of beta adrenergic blockade. Cardiovasc Drugs Ther. 1989;3:545-555. [Medline] [Order article via Infotrieve]
7.
Asano M, Uchida W, Shibasaki K, Terai M, Inagaki O, Takenaka
T, Matsumoto Y, Fujikura T. Pharmacological profile of YM-16151-1 and
its optical isomers: a novel calcium entry blocking and selective
beta-1 adrenoceptor blocking agent. J Pharmacol Exp Ther. 1990;254:204-211.
8.
Casati C, Monopoli A, Dionisotti S, Zocchi C, Bonizzoni E,
Ongini E. Repeated administration of selective adenosine A1 and A2
receptor agonists in the spontaneously hypertensive rat: tolerance
develops to A1-mediated hemodynamic effects. J Pharmacol Exp
Ther. 1994;268:1506-1511.
9. Conti A, Monopoli A, Gamba M, Borea PA, Ongini E. Effects of selective A1 and A2 adenosine receptor agonists on cardiovascular tissues. Naunyn Schmiedeberg's Arch Pharmacol. 1993;348:108-112. [Medline] [Order article via Infotrieve]
10. Monopoli A, Conti A, Dionisotti S, Casati C, Camaioni E, Cristalli G, Ongini E. The pharmacology of 2-chloro-N6-cyclopentyladenosine, CCPA, a highly selective A1 adenosine receptor agonist. Arzneimittelforschung. 1994;44:1305-1312. [Medline] [Order article via Infotrieve]
11. Monopoli A, Conti A, Zocchi C, Casati C, Volpini R, Cristalli G, Ongini E. The pharmacology of 2-hexynyl-5'-N-ethylcarboxamidoadenosine, 2HE-NECA, a new selective A2a adenosine receptor agonist. Arzneimittelforschung. 1994;44:1296-1304. [Medline] [Order article via Infotrieve]
12. Gibaldi M, Perrier D. Pharmacokinetics. 2nd ed. New York, NY: Marcel Dekker Inc; 1975:189-213.
13. SAS Institute Inc. SAS/STAT Guide for Personal Computers. 6th ed. Cary, NC: SAS Institute Inc; 1989:1135-1193.
14. Marquardt DW. An algorithm for least squares estimation of nonlinear parameters. J Soc Ind Appl Math. 1963;11:431-441.
15. Arunlankshana O, Schild HO. Some quantitative uses of drug antagonists. Br J Pharmacol. 1959;14:48-52. [Medline] [Order article via Infotrieve]
16. Draper NR, Smith H. Applied Regression Analysis. 2nd ed. New York, NY: John Wiley & Sons; 1981:141-183.
17.
Calhoun DA, Zhu S, Wyss JM, Oparil S. Diurnal blood pressure
variation and dietary salt in spontaneously hypertensive rats.
Hypertension. 1994;24:1-7.
18. Bazil MK, Lappe RW, Webb RL. Pharmacologic characterization of an endothelinA (ETA) receptor antagonist in conscious rats. J Cardiovasc Pharmacol. 1992;20:940-948. [Medline] [Order article via Infotrieve]
19. Bazil MK, Webb RL. Hemodynamic effects of amlodipine and benazeprilat in spontaneously hypertensive rats. J Cardiovasc Pharmacol. 1993;22:897-905.
20. Abou-Mohamed G, Caldwell RW, Ibrahim TM, Tuttle RR. Integrative cardiovascular actions of a novel catecholamine, GP-2-128. J Cardiovasc Pharmacol. 1994;23:485-491. [Medline] [Order article via Infotrieve]
21. Kanda A, Haruno A, Miyake H, Nagasaka M. Antihypertensive effects of MPC-1304, a novel calcium antagonist, in experimental hypertensive rats and dogs. J Cardiovasc Pharmacol. 1992;20:723-730. [Medline] [Order article via Infotrieve]
22.
Ludbrook J. Repeated measurements and multiple comparisons in
cardiovascular research. Cardiovasc Res. 1994;28:303-311.
23.
Webb RL, McNeal RB, Barclay BW, Yasay GD. Hemodynamic effects
of adenosine agonists in the conscious spontaneously hypertensive rat.
J Pharmacol Exp Ther. 1990;254:1090-1099.
24.
Webb RL, Barclay BW, Graybill SC. Cardiovascular effects of
adenosine A2 agonists in the conscious spontaneously hypertensive rat:
a comparative study of three structurally distinct ligands. J
Pharmacol Exp Ther. 1991;259:1203-1212.
This article has been cited by other articles:
![]() |
J.-N. Yang, J.-F. Chen, and B. B. Fredholm Physiological roles of A1 and A2A adenosine receptors in regulating heart rate, body temperature, and locomotion as revealed using knockout mice and caffeine Am J Physiol Heart Circ Physiol, April 1, 2009; 296(4): H1141 - H1149. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. K. Dhalla, M.-Y. Wong, W.-Q. Wang, I. Biaggioni, and L. Belardinelli Tachycardia Caused by A2A Adenosine Receptor Agonists Is Mediated by Direct Sympathoexcitation in Awake Rats J. Pharmacol. Exp. Ther., February 1, 2006; 316(2): 695 - 702. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Negishi, J.-W. Xu, K. Ikeda, M. Njelekela, Y. Nara, and Y. Yamori Black and Green Tea Polyphenols Attenuate Blood Pressure Increases in Stroke-Prone Spontaneously Hypertensive Rats J. Nutr., January 1, 2004; 134(1): 38 - 42. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Monopoli, C. Casati, G. Lozza, A. Forlani, and E. Ongini Cardiovascular Pharmacology of the A2A Adenosine Receptor Antagonist, SCH 58261, in the Rat J. Pharmacol. Exp. Ther., April 1, 1998; 285(1): 9 - 15. [Abstract] [Full Text] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Hypertension Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 1995 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |