![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
NEUROPHARMACOLOGY
College of Science and Technology and School of Medicine, Temple University, Philadelphia, Pennsylvania
Received February 20, 2004; accepted June 1, 2004.
| Abstract |
|---|
|
|
|---|
When the isobole for a specified effect has been determined, the dose pair of the combination that experimentally produces this effect may plot as a point that is below or above the additive isobole line and thus indicate superadditivity or subadditivity, respectively. The graph containing the additive line and the plotted experimental point constitutes a graph known as an isobologram.
Since its employment by Loewe (1953
) the isobologram has been used in many studies of drug combinations (e.g., Gerthoffer, 1996
; Wei and Roerig, 1998
; Fairbanks and Wilcox, 1999
; Hurley et al., 1999
; Bolan et al., 2002
; Field et al., 2002
; Cichewicz and McCarthy, 2003
) as well as in experimental designs that administer a single drug at two different sites (Yeung and Rudy, 1980
). The popularity of the isobologram is undoubtedly related to the simplicity of constructing the line of additivity (linear under the assumption of constant relative potency) and the consequent visual assessment it affords when the experimentally derived dose combination for the specified effect is plotted. That experimental dose combination, usually obtained from a regression analysis of the combination, is then compared statistically with the corresponding combination dose on the additive line to determine a significant difference. Those statistical tests and additional discussions of isobolograms may be found in the monograph by Tallarida (2000
) and in previous reviews (Tallarida, 1992
, 2001
).
The main subject of this communication is the analysis of a drug combination displaying variable relative potency as is the case when two drugs (or ligands) yield different maximum effects, e.g., in experiments that employ a full agonist and a partial agonist. In this case the relative potency is not constant and we show here that the isoboles are not straight lines, a fact that has not been adequately addressed in the literature. Specifically, we show here how the fitted dose-effect curves of the individual agents are used in the calculation of the additive isoboles and, thus, the expected effect for any dose combination. We begin by examining the usual situation of constant R and subsequently derive and apply the equations to data with variable R.
| Theory |
|---|
|
|
|---|
![]() |
![]() |
![]() |
These are hyperbolas that appear S-shaped when the dose scale is calibrated logarithmically as in Fig. 1. The combination dose (a,b) of drugs A and B is equivalent to the dose b + b'of drug B, where b' is obtained by equating 100b'/(b' + B50) to Eca/(a + Ac).
|
Thus,
![]() | (1) |
![]() | (2) |
|
| Results |
|---|
|
|
|---|
Results of a Hypothermic Combination. If one of the two compounds lacks efficacy then the concept of additivity means that its concomitant use with the active compound results in no change in the dose-effect relation of the active compound. However, higher doses of the "inactive" compound may produce measurable effects and, thus, the use of the higher dose range presents a case of variable relative potency. In that regard, we illustrate with the use of data obtained in a previous study (Rawls et al., 2002
). As will be shown, the data from that study also require a more general curve-fitting procedure that, in turn, necessitated a generalization of eq. 2. It is, therefore, a very useful illustration. That study examined the effect of the NMDA antagonist dextromethorphan (DXM) on the hypothermic response to WIN 55212-2, a selective cannabinoid agonist, in rats. The cannabinoid produced dose-dependent hypothermia in the dose range of 1 to 10 mg/kg, i.m., whereas DXM evoked dose-dependent hypothermia only for doses equal to or greater than 30 mg/kg (i.m). [However, the lower doses of DXM (10 mg/kg) potentiated the hypothermic response to WIN 55212-2 (1, 2.5, or 5 mg/kg).] We use that data here to illustrate the calculations for the additive effect of a combination of WIN 55212-2, the full agonist, and DXM in its higher dose range, where it produces modest dose-dependent effects.
The dose-effect data for the partial agonist, DXM, and the full agonist, WIN 55212-2, were not well described by simple hyperbolas and, thus, required equations of the form E = EcAq/(Aq + Acq) and E = EBBp/(Bp + B50p), respectively, where p and q are curve-fitting parameters (Hill coefficients). These graphs are shown in Fig. 3. The term EB is the maximum effect (previously expressed as 100%), whereas all other terms have the same meaning as previously defined in eq. 2. In this case the isobolar equation corresponding to eq. 2 is
![]() | (3) |
|
Parameters of these fits for the hypothermic data were as follows: Ec = 1.58; q = 1.92; Ac = 65.8; EB = 4.17; p = 1.73; B50 = 3.99. Isoboles for several different effect levels, shown in Fig. 4, were calculated from eq. 3. Because the theoretical maximum effect of DXM is 1.58, the additive isoboles for the two largest effects shown do not have intercepts on the horizontal axis. Although this example was chosen mainly to illustrate the calculations, the additive isoboles derived provide a basis for future combination experiments with WIN 55212-2 that employ DXM in this higher dose range.
|
Variable Relative Potency for Two Full Agonists. The combination of a full and partial agonist presented above represents a clear example of a varying potency ratio, one that becomes extreme as the effect level approaches Ec. But a variable potency ratio also occurs with two full agonists. To illustrate the isobole equations (and their graphs) when two full agonists have a varying potency ratio, we have used the parameters from the previous hypothermic example but adjusted the efficacies to be the same. This was done to illustrate a common situation in which two full agonists have nonparallel regressions of effect on log dose. Toward that end we used values of drug B (the more potent drug) that give several different effect levels. Figure 5 shows the isoboles for these effects. It is seen that these are slightly curved and are typical of this situation of nonparallel regressions. Further mathematical detail for this case is provided in the Appendix. In contrast to the combination of drugs of different efficacies, this combination leads to isoboles that reach both axes and, because q and p differ only slightly, the curvature is minimal. When q/p > 1, as in this illustration, the isoboles (plotted as b against a) are concave downward. For q/p < 1, the isoboles are concave upward.
|
| Discussion |
|---|
|
|
|---|
The term denoted by
, defined in the Appendix, has a relation to the relative potency that becomes more evident when the dose component a of drug A is small. It is seen (for the simple hyperbolic case and low doses) that
is a nondimensional parameter that indicates the relative potency of drugs A and B: one unit of drug B is equivalent to
units of drug A. More generally, one unit of drug B is equivalent to an amount of drug A given by
1/q. Most published reports that use isobolographic analysis have involved drugs that are equal, or approximately equal, in efficacy. In certain other studies, one of the drugs has zero efficacy, a case that is readily analyzed by comparing the dose-effect curve of the active agent alone and when it is in combination as previously described. A significant difference then indicates a departure from simple additivity. That situation applied to low doses of DXM given concomitantly with WIN 55212-2. Another early example is contained in a study by Porreca et al. (1990
) in which it was shown that peripheral [Leu5]enkephalin had virtually no antinociceptive efficacy, but when administered with morphine, the combination enhanced the effects of morphine.
The combination of a full and a partial agonist, the main topic of this investigation, has not received adequate attention. The analysis presented here shows that this situation leads to curvilinear isoboles. That curvature could be misinterpreted to suggest synergism (or subadditivity) when the combination is actually additive.
Additivity is based on the concept of dose equivalence, i.e., how much of one of the drugs is equivalent to a dose of the other. This equivalence is the basis of the relation derived for the additive isobole so that the combination doses, a and b, can be expressed as a dose of either one of them. In our derivation, we chose to demonstrate the conversion of dose a of drug A into an equivalent of drug B. Thus, dose b plus this equivalent becomes the dose of drug B that yields the expected effect when this summed dose is used in the dose-effect equation of drug B. All tested dose pairs provide additive effect values, thereby producing the additive dose-effect relation. This relation is then viewed against the actual dose-effect relation that is experimentally determined. The additive and actual dose-effect relations, expressed in the usual way as regressions of effect on log (dose), thereby provide the additive and the experimental doses for the effect level selected (usually 50% of Emax). These are statistically compared as described by Tallarida (2000
). A significant difference indicates a departure from simple additivity, either synergism or subadditivity depending on the relative values of the additive and actual dose values. Combinations of the kind discussed here, i.e., those with variable relative potency, are likely to occur in detailed studies of mechanism, especially in mechanistic studies involving endogenous compounds. For example, endogenous opioids (enkephalins, endorphins, and dynorphin) have been the subject of many studies of their ability to activate G proteins as measured by the binding of [35S]GTP
S to membranes. Some of these compounds show high efficacy when compared with common agonists such as morphine and meperidine (Alt et al., 1998
). Marked differences in both efficacy and potency have also been seen among µ opioids using forskolin-stimulated cAMP production as shown by Gharagozlou et al. (2003
). That study demonstrated the rank of order of efficacy of certain compounds as follows:
endorphin = fentanyl > etorphine = morphine > cyclazocine. In contrast, their potencies were different in rank order: etorphine > cyclazocine > fentanyl > morphine [tmt]
endorphin. As more studies of intimate mechanisms are conducted it is likely that additional combinations involving endogenous and exogenous agents will be examined for interactions, and because of these marked differences in efficacy and potency, the analysis presented here is especially applicable. A further application is to studies involving a drug and its active metabolite(s). An example is afforded by buprenorphine, which produces an active metabolite, norbuprenorphine, of lower efficacy (Ohtani et al., 1995
; Cowan, 2003
). The difference in efficacies requires the methodology presented here in order to detect possible interactions between the compounds, an area currently being pursued in our laboratories. Aside from these mechanistic-type studies there will continue to be interest in joint drug application, and when these show varying relative potency, it is necessary to calculate the additive isoboles as described here (eqs. 2 and 3). Because the additive isoboles define the expected response surface of a combination experiment, this approach is quite useful in studying combination drug action. The use (misuse) of linear isoboles in cases of variable relative potency could lead to incorrect conclusions on the nature of the possible interaction of the agents.
| APPENDIX |
|---|
|
|
|---|
![]() | (A1) |
![]() |
Bc, and the additive isobole starting at (0,B) is an unbounded hyperbolic arc that decreases toward the horizontal asymptote b = B Bc.
When B
Bc, meaning effect levels below Ec, the additive isoboles are arcs that join intercepts (0,B) and (A,0) where A and B are unitary doses that are equally effective.
Intercept A may also be expressed by
![]() |
is evident; it is seen that for very small doses, hence very small values of A and B, this expression gives A =
B, thereby illustrating that
is a relative potency term for small doses; i.e., one unit of drug B is equivalent to
units of drug A.
Curve-Fitting Parameters. For the hypothermic combination, which uses curve-fitting parameters p and q, the form given by eq. A1 becomes
![]() | (A2) |
= (EBAcq)/(EcB50p). In this case, we also get an interpretation of
by considering very small doses: Aq =
Bp. For effects >Ec the additive isobole starts at (0,B) and is an unbounded arc that decreases toward the asymptote given by b = B Bc. For effects <Ec, the additive isobole originates at (0,B) and terminates at (A,0) where intercept A is given by
![]() |
Two Full Agonists with Different Regression Line Slopes. A compound whose dose-effect data are modeled with a curve-fitting parameter p according to E = 100Bp/(Bp + B50p), commonly has the data displayed as (log B,E), and these data are fitted by a linear regression of E on log B. (Common logarithms are most often used.) This transformation of dose to the logarithmic scale leads to the "logistic" form:
![]() | (A3) |
This common procedure, when limited to effects in the midrange, has a slope that is proportional to p. Specifically, at the 50% effect level it may be shown that the slope = 25 ln(10) p (approximately 57.6 p). Thus, two full agonists with different curve-fitting parameters, p and q, will have regression line slopes that differ. This situation is illustrated in the example described in the text and illustrated in Fig. 5. In this case each maximum is normalized to 100 and thus Ac = A50. The isobole equation corresponding to eq. A2 becomes
![]() |
![]() |
![]() | (A4) |
| Acknowledgements |
|---|
| Footnotes |
|---|
Article, publication date, and citation information can be found at http://jpet.aspetjournals.org.
ABBREVIATIONS: NMDA, N-methyl-D-aspartate; DXM, dextromethorphan; WIN 55212-2, (R)-(+)-[2,3-dihydro-5-methyl-3-(4-morpholinylmethyl)pyrrolo[1,2,3-de]-1,4-benzoxazin-6-yl]-1-naphthalenylmethanone; [35S]GTP
S, guanosine 5'-O-(3-thio)triphosphate.
Address correspondence to: Dr. Ronald J. Tallarida, Dept. of Pharmacology, Temple University School of Medicine, 3420 N. Broad St, Philadelphia, PA 19140. E-mail: ronald.tallarida{at}temple.edu
| References |
|---|
|
|
|---|
Alt A, Mansour A, Akil H, Medzihradsky F, Traynor JR, and Woods JH (1998) Stimulation of guanosine-5'-O-(3-[35S]thio) triphosphate binding by endogenous opioids acting at a cloned mu receptor. J Pharmacol Exp Ther 286: 282-288.
Bolan EA, Tallarida RJ, and Pasternak G (2002) Synergy between mu receptor subtypes: evidence for functional interaction among mu receptor subtypes. J Pharmacol Exp Ther 303: 557-562.
Cichewicz DL and McCarthy EA (2003) Antinociceptive synergy between delta(9)tetrahydrocannabinol and opioids after oral administration. J Pharmacol Exp Ther 304: 1010-1015.
Cowan A (2003) Buprenorphine: new pharmacological aspects. Int J Clin Pract Suppl 133: 3-8.
Fairbanks CA and Wilcox GL (1999) Spinal antinociceptive synergism between morphine and clonidine persists in mice made acutely or chronically tolerant to morphine. J Pharmacol Exp Ther 288: 1107-1116.
Field MK, Gonzalez MI, Tallarida RJ, and Singh L (2002) Gabapentin and the NK1 antagonist CI-1021 act synergistically in two rat models of neuropathic pain. J Pharmacol Exp Ther 303: 730-735.
Gerthoffer WT (1996) Agonist synergism in airway smooth muscle contraction. J Pharmacol Exp Ther 278: 800-807.
Gharagozlou P, Demirci H, Clark D, and Lameh J (2003) Activity of opioid ligands in cells expressing cloned mu opioid receptors. BMC Pharmacol 3: 1-15.[CrossRef][Medline]
Hurley RW, Grabow TS, Tallarida RJ, and Hammond DL (1999) Interaction between medullary and spinal delta-1 and delta-2 opioid receptors in the production of antinociception in the rat. J Pharmacol Exp Ther 289: 993-999.
Loewe S (1953) The problem of synergism and antagonism of combined drugs. Arzneim-Forsch 3: 385-390.
Ohtani M, Kotaki H, Sawada Y, and Iga T (1995) Comparative analysis of buprenorphine- and norbuprenorphine-induced analgesic effects based on pharmacokinetic-pharmacodynamic modeling. J Pharmacol Exp Ther 272: 505-510.
Porreca F, Jiang Q, and Tallarida RJ (1990) Modulation of morphine antinociception by peripherally-administered [Leu-5] enkephalin: demonstration of a synergistic interaction. Eur J Pharmacol 179: 463-468.[CrossRef][Medline]
Rawls SM, Cowan A, Tallarida RJ, Geller EB, and Adler MW (2002) N-Methyl-D-aspartate antagonists and WIN 55212-2 [4,5-dihydro-2-methyl-4(4-morpholinylmethyl)-1-(1-naphthalenyl-carbonyl)-6H-pyrrolo[3,2,1-i,j]quinolin-6-one], a cannabinoid agonist, interact to produce synergistic hypothermia. J Pharmacol Exp Ther 303: 395-402.
Tallarida RJ (1992) Statistical analysis of drug combinations for synergism. PAIN 49: 93-97.[CrossRef][Medline]
Tallarida RJ (2000) Drug Synergism and Dose-Effect Data Analysis, Chapman-Hall/CRC Press, Boca Raton, FL.
Tallarida RJ (2001) Drug synergism: its detection and applications. Perspectives in pharmacology. J Pharmacol Exp Ther 298: 865-872.
Wei ZY and Roerig SC (1998) Spinal morphine/clonidine antinociceptive synergism is regulated by protein kinase C, but not protein kinase A activity. J Pharmacol Exp Ther 287: 937-943.
Yeung JC and Rudy TA (1980) Multiplicative interaction between narcotic agonisms expressed at spinal and supraspinal sites of antinociceptive action as revealed by concurrent intrathecal and intracerebroventricular injections of morphine. J Pharmacol Exp Ther 215: 633-642.
This article has been cited by other articles:
![]() |
A. S. Braverman, R. J. Tallarida, and M. R. Ruggieri Sr. The Use of Occupation Isoboles for Analysis of a Response Mediated by Two Receptors: M2 and M3 Muscarinic Receptor Subtype-Induced Mouse Stomach Contractions J. Pharmacol. Exp. Ther., June 1, 2008; 325(3): 954 - 960. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. D. Fischer, E. I. Zimmerman, M. J. Picker, and L. A. Dykstra Morphine in Combination with Metabotropic Glutamate Receptor Antagonists on Schedule-Controlled Responding and Thermal Nociception J. Pharmacol. Exp. Ther., February 1, 2008; 324(2): 732 - 739. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. J. Tallarida An Overview of Drug Combination Analysis with Isobolograms J. Pharmacol. Exp. Ther., October 1, 2006; 319(1): 1 - 7. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. I. Lorenzo and P. Sanchez-Marin Comments on "Isobolographic Analysis for Combinations of a Full and Partial Agonist: Curved Isoboles" J. Pharmacol. Exp. Ther., January 1, 2006; 316(1): 476 - 478. [Full Text] [PDF] |
||||
![]() |
R. J. Tallarida Response to Comments on "Isobolographic Analysis for Combinations of a Full and Partial Agonist: Curved Isoboles" J. Pharmacol. Exp. Ther., January 1, 2006; 316(1): 479 - 479. [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||