| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
BACKGROUND: Inhaled anesthetics have been postulated to act at multiple receptors, with modest action at each site summing to produce immobility to noxious stimulation. Recent experimental results affirm prior findings that inhaled anesthetics interact additively. Synergy implies multiple sites of action by definition. In this essay, we explore the converse: does additivity imply a single site of action? METHODS: The interaction of one versus two ligands competing for the same binding site at a receptor was explored using the law of mass action. Circuits were then constructed to investigate how the potency of drugs and the steepness of the concentration versus response relationship is amplified by the arrangement of suppressors into serial circuits, and enhancers into parallel circuits. Assemblies of suppressor and enhancer circuits into signal processing units were then explored to investigate the constraints signal processing units impose on additive interactions. Lastly, the relationship between synergy, additivity, and fractional receptor occupancy was explored to understand the constraints imposed by additivity. RESULTS: Drugs that compete for a single receptor, and that similarly affect the receptor, must be additive in their effects. Receptors that bind suppressors in serial circuits, or enhancers in parallel circuits, increase the apparent potency of the drugs and the steepness of the concentration versus response relationship. When assemblies of suppressor and enhancer circuits are arranged into signal processing units, the interactions may be additive or synergistic. The primary determinant is the relationship between the concentration of drug associated with the effect of interest and the concentration associated with 50% receptor occupancy, kd. Effects mediated by very low concentrations are more likely to be additive. Similarly, inhaled anesthetics that act at separate sites are unlikely to exhibit additive interactions if anesthetic drug effect occurs at concentrations at or above 50% receptor occupancy. However, if anesthetic drug effect occurs at very low levels of receptor occupancy, then additivity is expected even among anesthetics acting on different receptors. CONCLUSIONS: Additivity among drugs acting on different receptors is only likely if the concentrations responsible for the drug effect of interest are well below the concentration associated with 50% receptor occupancy.
Several investigations studying anesthetic pairs find that inhaled anesthetics interact with other inhaled anesthetics in an additive manner.1–4 Using the minimum alveolar concentration (MAC) that suppresses movement to noxious stimulation as the measure of anesthesia drug effect, a companion manuscript in this issue of Anesthesia & Analgesia describes the additivity of 16 inhaled anesthetic pairs with divergent effects on known ion channels, and with a rigorous definition of additivity.5 The authors suggest that the lack of synergy is consistent with a unitary mechanism of inhaled anesthetic action, rather than a mechanism requiring diverse actions synergistically summing to generate immobility in response to noxious stimulation. A second companion paper in this issue of Anesthesia & Analgesia explores the interactions among many drug classes of IV and inhaled drugs used in anesthetic practice.6 The general rule, albeit with exceptions, is that drugs with different sites of action interact synergistically. Eckenhoff and Johansson7 proposed that inhaled anesthetics produce immobility by summation of modest effects on multiple receptor systems. Given the synergy observed in human and animal trials with many drug combinations that act on different receptors,6 if Eckenhoff and Johansson are correct we would expect at least some inhaled anesthetic pairs that have divergent effects on specific receptors would act synergistically. However, this is unambiguously not the case.5 The interactions are additive. This essay explores the hypothesis that synergy implies separate sites of action, and the corollary hypothesis that additivity implies a unitary mechanism of action.
Definitions Drug interactions are classically defined using isobolographic analysis. Assessment of the interaction requires that we precisely define drug effect. It is experimentally and mathematically convenient to select half of the maximum possible effect because concentration versus effect curves are typically steepest at half-maximal effect, permitting the most precise measurement of the concentration associated with the stated effect, while requiring the fewest animals to assess the concentration associated with a 50% probability of an animal generating a binary (yes/no) response (e.g., movement to noxious stimulation). For inhaled anesthetics, the effect is usually MAC, the alveolar concentration at steady-state associated with a 50% probability of moving in response to noxious stimulation. As shown in Figure 1, the next step is to determine the concentration of each drug that achieves the effect of interest. Concentration may be an instantaneous plasma concentration, a steady-state plasma concentration, or an effect-site concentration. For inhaled anesthetics, the "concentration" is often the steady-state end-tidal concentration of the anesthetic. For IV anesthetic drugs, "concentration" is typically the calculated effect-site concentration, which is, by definition, the plasma concentration at steady-state. In some studies a bolus dose or infusion rate may substitute for concentration.
In Figure 1, the X and Y axes units are fractions of the concentration of drugs A and B, respectively, that produce the defined drug effect. For inhaled anesthetics, the "concentration" in Figure 1 corresponds to the MAC fraction. "Isoboles" are lines that define the concentrations (in the present case, the concentrations normalized to the concentration that produces the target effect) of two drugs that, given together, produce exactly the same effect. For an additive relationship, a proportional increase in the concentration of drug B exactly compensates for a proportional reduction in the concentration of drug A, resulting in a straight (additive) isobole. For example, the point marked X on Figure 1 implies that half of the required concentration of drug A (e.g., a 0.5 fraction of the concentration of drug A that produces the target effect) given with half of the required concentration of drug B (e.g., the same 0.5 fraction for drug B), produces the drug effect produced by the required concentration of drug A alone or the required concentration of drug B alone. This defines additivity. Experimentally, additivity is considered to result if test findings lie within a defined area surrounding the additivity isobole. Evidence for synergy results if the combined (normalized) concentrations of both drugs required to produce the stated drug effect are significantly less than predicted by additivity (Fig. 1). Evidence for "infra-additivity" results if the combined concentrations of both drugs required to produce the stated drug effect significantly exceed that predicted by additivity. Infra-additivity differs from classical pharmacological antagonism. For example, naloxone is a classical pharmacological antagonist of opioid drug effect. No concentration of naloxone causes analgesia, which precludes constructing an isobole that relates naloxone analgesia and opioid analgesia. What we will call "infra-additivity" in this analysis is often called "antagonism" in drug interaction studies. We have no quibble with this definition, but in this essay we use the term "infra-additivity." Many investigators define synergy or infra-additivity as occurring when a combination differs significantly from the line of additivity. We will focus our analysis on three points of synergy in Figure 2. Point X demonstrates statistically significant synergy, as implied by the confidence bars for the respective concentrations of drugs A and B. The synergy demonstrated by point X is probably clinically significant. Point Y also deviates significantly from the line of additivity, but this statistically significant synergy is so slight that clinically it cannot be distinguished from additivity. Point Z is the most problematic. Despite a distance from point Z to the line of additivity that suggests a clinically important synergy, statistical significance is not achieved. Many investigators would conclude that, since point Z does not differ statistically significantly from the line of additivity, point Z represents additivity. However, point Z may demonstrate true synergy, as discussed below.
For this essay, and the companion papers, we define "synergy" as resulting when the sum of normalized fractional concentrations of two drugs that combined produce a target effect (e.g., MAC) is significantly <0.9 (i.e., a 10% or more decrease from additivity that reaches statistical significance). The 10% definition is arbitrary, but is consistent with the companion manuscripts.5,6 Similarly, there is "additivity" when the fractional concentrations of the two drugs, in combination, sum to between 0.9 and 1.1. Finally, there is "infra-additivity" when the fractions of the concentration of each drug individually that produces the same effect, sum to significantly exceed 1.1 (i.e., a 10% increase from additivity that reaches statistical significance). If the fractional sum of the concentrations is <0.9, or exceeds 1.1, but this difference is not statistically significant, we define this as "additivity by default" meaning that the data do not allow exclusion of synergy (if <0.9) or infra-additivity (if more than 1.1). The most accurate assessment of synergy, additivity, and infra-additivity is achieved by giving approximately equipotent concentrations, indicated by the line of equipotent dosing in Figure 1. Detection of synergy becomes difficult or impossible if the concentration of either drug approaches the concentration that given alone produces the defined effect. By our definition, synergy cannot be detected if one gives 91% of the concentration of drug A that produces the desired effect because the concentration fraction to produce that effect is already 0.91 (i.e., within 10% of 1.0). The sum of the concentration fractions, even if the concentration of drug B is 0, still exceeds 0.9, eliminating synergy. Referring again to Figure 2, X represents "synergy" by our definition because the fractional concentrations sum to <0.9, and the 95% confidence bounds exclude the line of additivity. Although the 95% confidence bounds for Y exclude the line of additivity, Y represents "additivity" because the fractional sum lies between 0.9 and 1.1, Z represents "additivity by default." Z strongly suggests synergy. The relationship suggested by point Z may be synergistic, but the study has insufficient statistical power to assert this at P < 0.05. However, at the same time it is not correct to conclude that Z represents additivity. Indeed, because point Z is the location with the highest probability of defining the relationship, the relationship is probably synergistic. By definition, additivity is how a drug interacts with itself. The combination of a normalized 0.5 concentration of drug A with a normalized 0.5 concentration of drug A, obviously produces exactly the same effect as a 1.0 normalized concentration of drug A, because the subject received a 1.0 normalized concentration of drug A. By tautology, a drug's interaction with itself must be additive. Note that "synergy," "additivity," and "infra-additivity" always refer to adding drug concentrations (or doses). "Synergy," "additivity" and "infra-additivity" are inappropriate for discussing drug effects. Most concentration versus response relationships are nonlinear (e.g., doubling the concentration does not double the effect). The effect of adding two concentrations, either of one drug or multiple drugs, can deviate from the sum of the individual effects, depending on the position of each concentration in the concentration versus response relationship. Reverting briefly to the dose domain, 100 mg of vecuronium, added to 100 mg of rocuronium, produces a less than additive effect. Each dose produces 100% blockade, but the combination can only produce 100% blockade, not 200% blockade. Thus, the effects are not additive. Indeed, most drugs show a ceiling effect. Combining two drugs in the same class at doses that individually produce the ceiling effect would still give the same ceiling, which is clearly less than additivity. So "additivity" and "synergy" strictly refer to adding drug concentrations (or doses), never adding drug effects. Returning to Figure 1, the lines separating additivity from synergy and infra-additivity indicate 10% deviations from the additivity isobole. The region of additivity is just 13% of the total area of the graph. Indeed, the only meaningful deviation is that which occurs along the line of equipotent dosing, and that is just 10% of the total possible deviation from additivity. This is the first suggestion of the punch line of the manuscript: additivity is unusual, except for a drug acting on itself, which is additive by definition.
Mass Action and Synergy: Individual Receptors to Simple Circuits to Single Processing Units The point relative to our examination of anesthetic mechanisms is that competition of full agonists for the same site of action will yield a strictly additive interaction. We will now explore the converse: when can two drugs display additivity of normalized concentrations yet act at separate sites?
Receptors Grouped into Circuits
Circuits Grouped into Signal Processing Units
Several parameters affect the performance of these signal processing units, including the concentration of drug associated with the effect of interest, the probability of a signal propagating past an unbound receptor, the number of receptors at Site A and Site B, and the efficacy of a drug in producing the maximal possible response (either reducing the probability of signal propagation to 0 for a suppressor, or increasing the probability of signal propagation to one for an enhancer). Appendix C demonstrates that the primary factor that determines whether the circuit demonstrates additivity or synergy is the concentration of suppressor or enhancer associated with the drug effect of interest, normalized to kd, the dissociation constant (i.e., the concentration associated with 50% receptor occupancy). If the normalized concentration is <1%, then additivity is expected. If the normalized concentration is in the range of (or exceeds) kd, the concentration associated with 50% receptor occupancy, then synergy is expected. The probability of signal propagation past an unbound receptor, the number of receptors in series or parallel, the intrinsic efficacy, and the geometry of the circuit are of considerably less importance than the relationship between the concentration associated with the drug effect of interest relative to the concentration associated with 50% receptor occupancy.
Receptor Occupancy Versus Concentration–Effect Relationships
, where f is the fraction of receptors occupied (see Y axis on the left), [L] is the concentration of the ligand, and kd is the dissociation constant. The axis marked "concentration" shows the concentration of the ligand, normalized to the concentration producing 50% receptor occupancy, i.e., "concentration" = [L]/kd. As a result, the concentration at "X" is 1.0, which is the concentration associated with 50% occupancy. The steeper curve on the right is a hypothetical curve for inhaled anesthetic effect defined by the relationship
where Pno response is the probability of no movement (see Y axis on the right), MAC is 10 (i.e., 10 fold higher than the concentration producing 50% receptor binding), and
Now consider two concentrations of the same inhaled anesthetic, designated by the two vertical arrows. A concentration of five (
Figure 5 shows the same two curves, except that MAC is now 0.1, one-tenth of the concentration associated with 50% receptor occupancy. Again consider two concentrations of a single drug. The first concentration (0.05; first upward arrow) produces a 5% receptor occupancy and, as before, a negligible 6% probability of no response. The second concentration (0.1), "SYMBOL", is double the first and increases receptor occupancy to 9%, again with a huge increase in drug effect to 50% probability of no response.
As previously shown, neuronal circuits amplify drug potency and increase the steepness of concentration versus response relationships, compatible with the predictions of Figure 5. For this reason, MAC probably is far below kd (as illustrated in Fig. 5), rather than greatly exceeding kd (as in Fig. 4.) Additionally, a MAC occurring at concentrations far exceeding 50% receptor occupancy (Fig. 4) would require a circuit exquisitely tuned to abruptly amplify the signal with almost no change in fractional receptor occupancy. Further diminishing the likelihood of this scenario, the many organisms that show similar sensitivity to inhaled anesthetics would need to conserve such circuitry, to no apparent purpose.
Two Drug Interaction
The lower example in Figure 6 considers two targets of drug action, one for drug A, and one for drug B. As before, a
Figure 7 shows two possibilities for drug interactions when MAC is 10-fold less than the concentrations associated with 50% receptor occupancy. In the upper example, there is a single target of drug action, just as in Figure 5. A concentration of
The lower example considers two targets of drug action, one for drug A, and one for drug B. As before, It is easy to envision that two anesthetics that each produce 5% receptor occupancy at separate targets could have the same effect as twice the concentration of either drug alone acting on a single target to produce 9% receptor occupancy. Either way, the system will require considerable gain to translate the receptor occupancy (9% at a single target, versus 5% at 2 separate targets) to the huge increase in drug effect. However, the complexity of the circuitry is entirely in increasing the gain, which is an expected property of neural circuitry: additivity when MAC << kd does not require the complex constraints that exist if MAC vastly exceeds kd. This analysis suggests that if MAC exceeds kd, then additive interactions among inhaled anesthetics strongly imply a single site of anesthetic action. However, if MAC is well below kd, then additive interactions among inhaled anesthetics do not necessarily imply multiple sites of action.
Kaminoh et al. examined pathways of receptors arranged in combined serial and parallel circuits and performed an analysis similar to ours.9 They considered the effect of a single antagonist. Similar to our results, receptors in series increased the potency of antagonists. They also found that receptors in parallel decreased the potency of antagonists, an analysis we did not perform. Most critically, they found that it required a combination of parallel and serial circuits of antagonists to generate the steep concentration response curve associated with inhaled anesthetics. For example, a circuit consisting of 16 parallel pathways, each pathway containing 16 elements in series, would have a Hill coefficient of 4.5, similar to what is observed in MAC determinations for inhaled anesthetics. However, they and we predicted that circuits consisting only of elements in series, or elements in parallel, cannot generate a Hill coefficient more than 1.5. This is not surprising, because we used the same mathematical analysis for the probability of signal propagation, both over units of the circuit, and for the arrangement of individual circuits in serial and parallel assemblies. The Kaminoh et al. analysis, our analysis, and the analysis by Eckenhoff and Johansson7all suggest that the steepness of the concentration versus response relationship for anesthetics is partly derived from multiple receptors acting in circuits. Figures A2 and A3 show that the increasing steepness is associated with increasing potency. Thus, steep concentration versus response relationships for anesthetic drugs imply that a half-maximal effect (i.e., the clinical effect) occurs at concentrations much less than kd. We developed simple circuit examples to explore interactions among drugs acting at separate sites. Such circuits provide theoretical examples to understand basic principles. We have no illusions that they describe real neurological circuits. However, our choice of simple parallel or serial circuits was not arbitrary. Each neuron receives multiple dendritic connections. Presuming that more than one of these is activated, multiple dendritic connections are parallel assemblies. Groups of neurons running in parallel also may receive input simultaneously, again representing parallel transmission. However, presynaptic receptors arranged along a single axon are arranged in series. Additionally, nerves connect with each other in series to propagate a signal. Since movement of an animal or human in response to noxious stimulation may occur at a distance from the stimulation (e.g., pinch the tail and the upper extremity moves), the signal must pass through serial assemblies of neurons, just as we have postulated serial assemblies. Although circuits are most readily envisioned as neural structures, the analysis need not relate to nerves at all. Signals can propagate within a single cell through cascading enzymatic reactions or ion fluxes. Signals may even propagate within a single protein through multiple sites of receptor attachment. Our analysis demonstrates that additivity is expected when the drug effect of interest occurs at concentrations associated with low fractional receptor occupancy. This conclusion was reached from two different lines of reasoning: analysis of the probability of signal propagation as a function of the concentration of suppressors or enhancers (Section 2, and Appendices A–D), and an analysis of the relationship between fractional receptor occupancy and clinical concentration versus response relationships (Section 3). The next step in this line of reasoning is to verify these predictions against known values of kd for different anesthetic drugs, and see whether our models accurately predict synergistic or additive reactions. Unfortunately, the value of kd is known for very few anesthetics. Indeed, the few reported values are suspect because of artificial experimental conditions.7For this reason, verification of these findings against experimental data is not included herein. Instead, it is a topic for further research. As have others,7 we conclude that the clinical potency of inhaled anesthetics is likely amplified by neuronal circuitry from the binding potency at the level of individual receptors. We further conclude that additivity is an expected result from competition of drug for a common binding site, or when the drug effect occurs at concentrations that are much less than the concentration associated with 50% receptor occupancy, kd. If the drug effect is only found at concentrations close to or exceeding kd, then the interactions are necessarily synergistic. Thus, inhaled anesthetics either share a single site of action, or the effect is a result of binding to multiple sites of action at concentrations that leave the vast majority of sites unbound.
The authors express their appreciation to the reviewers of this article for their considerable contributions. Insights and suggestions from the reviewers identified critical deficiencies in the first submission and resulted in virtually a complete reanalysis and rewriting of the article. Further efforts by the reviewers helped anchor the theoretical arguments and vet the presentation to make this complex analysis as clear as possible.
When a single drug binds to a receptor, the binding follows the law of mass action:
where [D] is the concentration of the drug, [RU] is the concentration of unbound receptor, and [DR] is the concentration of bound receptor (actually, these are activities rather than concentrations, but for this essay we use the latter term). The rate of formation of [DR] is
. At steady-state, which is nearly instantaneous, the net rate of formation is 0, and thus [D][RU]kon=[DR]koff. If we define kd, the dissociation constant, as koff/kon, and f, fractional receptor occupancy, as
, then
This is an asymptotic relationship that reaches 50% receptor occupancy when [D] = kd. This can be expanded to two drugs competing for the same receptor. With simple competition, drug A and drug B independently follow the law of mass action:
, and
. The steady-state relationships are therefore [DA][RU]kon,A=[DAR]koff,A and [DB][RU]kon,B=[DBR]koff,B. If we define kd,A as koff,A/kon,A, kd,B as koff,B/kon,B, and f, fractional receptor occupancy, as
, then
If we select 50% occupancy as the drug effect of interest, then we can solve for the concentration of drug B associated with 50% receptor occupancy as a function of drug A:
. If we normalize the concentration of drug A to kd,A,
, and if we normalize the concentration of drug B to kd,B,
, then we can simplify the relationship between the concentration of drugs A and B associated with 50% receptor occupancy to
This relationship is shown in Figure A1, and obviously defines strict additivity. Thus, the law of mass action dictates that if 1) two full agonists independently compete for the same receptor, and 2) fraction receptor occupancy determines drug effect, then the interaction between the two drugs follows strict additivity.
The derivations of equations 1, 2, and 3 are given in Appendix D, which is available as a web supplement at www.anesthesia-analgesia.org.
Consider the action of a drug that alters the probability of propagation of a signal. The "signal" could be an action potential, an intracellular cascade, or the stepwise conformational transformation of a single protein. The critical notions are that the "signal" propagates information over time, and that propagation is binary. The signal either propagates, or it doesn't. We will define P as the probability of signal propagation. We will define "suppressors" as drugs that decrease the likelihood of signal propagation, and "enhancers" as drugs that increase the likelihood of signal propagation. We define Pun as the probability that a signal will pass an unbound receptor. If a suppressor is applied, the probability (P) will decrease from Pun towards 0. If an enhancer is applied, the probability (P) will increase from Pun towards 1. We define E as the efficacy of the drug action. E defines how far the drug is capable of altering the probability of signal propagation towards the limits of P = 1 (for enhancers) or P = 0 (for suppressors). E ranges from 0 (no efficacy) to 1 (full efficacy). For notational convenience, we define the concentration of a suppressor, Su, in units of the concentration that produces 50% receptor binding:
when the ligand, [L] suppresses signal propagation. We similarly define the concentration of an enhancer, En, in units of the concentration that produces 50% receptor occupancy:
when the ligand, [L], enhances signal propagation. Note that [L] and kd must have the same units. Although this would typically be moles per liter, the actual units are not relevant since they cancel out in the definitions of Su and En. Although Su and En will be referred to as concentrations, they could also be considered fractions or multiples Kd, depending on whether they are less or greater than one. As mentioned, the probability that a signal will pass an unbound receptor is Pun. Pun is most readily envisioned as one for a suppressor, and 0 for an enhancer. However, it is included in the equations to provide greater generality. We will assume that if a particular receptor is bound to a suppressor or enhancer, then the suppressor or enhancer is 100% likely to exert its effect. In other words, the probability of a signal traversing a receptor bound to a suppressor is Pun – E · Pun, and the probability of a signal traversing a receptor bound to an enhancer is Pun + E · (1 – Pun). Since the probability that a single receptor is bound at a certain drug concentration is exactly the same as the fraction of receptors bound in a population at that concentration, the law of mass action states that the probability of binding for a suppressor is
, and the probability of binding for an enhancer is
. These equations can be derived by rearranging Equations 4 and 5 and substituting into Equation 1 (see Appendix D in the online supplementary material). Given the above definitions, the probability that a particular signal passes a receptor in the presence of a suppressor is
, where Su is the concentration of the suppressor, as defined above, and E is the efficacy of the suppressor at reducing the probability of propagation to 0. The probability that a particular signal passes a receptor in the presence of an enhancer is
, where En is the concentration of the enhancer, as defined above, and E is the efficacy of the enhancer at increasing the probability of propagation to one. The upper graph in Figure A2 shows multiple receptors arranged in series, each containing a binding site for a suppressor. The probability that a signal passes multiple receptors occupied by a suppressor is the product of the individual probabilities. If N is the number of receptors in series, then the probability that the signal propagates past N receptors is
To explore this relationship we arbitrarily define Pun as 1 and E as 0.2. The choice of 0.2 permits demonstration of amplification of drug effect with physiologically reasonable numbers of receptors. We can graph this relationship as a function of Su, the concentration of suppressor, and N, the number of receptors in series. The bottom graph in Figure A2 shows this relationship. Although the effect of binding a single receptor is modest, with each bound receptor only reducing the probability of propagation by 20%, the effect of binding multiple receptors is amplified over the course of transmission. With four suppressors in series, the probability of propagation decreases by 50%. With 16 suppressors in series the probability of propagation decreases to nearly 0. Figure A2 also shows that the increase in effect is accompanied by an increase in the apparent potency of the suppressor. For N = 1, half-maximal effect occurs at 1 (by definition). For N = 16, half maximal effect occurs at a concentration of approximately 0.25. The increasing effect and potency is accompanied by a modest increase in the steepness of the concentration versus response relationship. The relationship for N = 16 can be modeled with a standard Hill equation,
, where Su50 is the suppressor concentration associated with 50% maximum drug effect, Max is the maximum drug effect, and
, similar to the analysis by Eger et al.8 However, Su50 is 0.25, a roughly four fold increase in potency, and We next turn to an analysis of enhancers. We cannot simply invert our model of suppressor in series to a model of enhancers in series. In an inverted model, Pun would be 0, and the model would be
, where 0.2 is the incremental increase in probability of propagation when the enhancer is bound into the receptor, and En is the concentration of the enhancer, normalized to the concentration associated with 50% receptor binding. In the inverted model, the probability of signal propagation is only 0.2 when N = 1. The probability of signal propagation further decreases for N > 1. Even if one assumes a very high prior probability of signal propagation past the unbound receptor, e.g., P =
, if N = 50, and En = 1 (50% receptor occupancy), the probability of signal propagation is <1%. One can overcome this by permitting the enhancer to nearly saturate the receptor (e.g., En >> 1), but as explained in the manuscript, if En >> 1 then the additivity necessarily implies a single site of action. This is why the circuit for enhancers, shown in Figure A3, is arranged as a parallel circuit, rather than in a series circuit. Enhancers acting in parallel can be analyzed with similar assumptions to suppressors acting in series, except the overall probability of propagation in a series circuit of N receptors is PN, while the overall probability in parallel circuit of N receptors is 1 – (1 – P)N. Given N receptors in parallel, the probability that the signal will propagate in a parallel circuit of N receptors is
. We will postulate a weaker effect for an enhancer, as the parallel circuit permits more enhancers to participate in the overall propagation of a single signal. If we define Pun as 0 (the signal cannot propagate without the enhancer) and E as 0.1, we can graph this relationship as a function of En, the concentration of enhancer normalized to the concentration associated with 50% receptor binding, and N, the number of receptors in parallel. This relationship is shown in the bottom graph in Figure A3. Although the effect of each bound enhancer is weak, the effect is amplified by the parallel circuit. With 16 enhancers in parallel, the maximum probability of propagation increases to approximately 80%. With 256 enhancers in parallel, the probability of propagation reaches nearly 100%. This is accompanied by an increase in the apparent potency of the enhancer, as well as an increasingly steep concentration versus response relationship. Increasing the number of enhancers 256 fold increases the apparent potency by 40 fold, as evidenced in the leftward shift of the dose versus response relationship in Figure A3. This is accompanied by an increase in the steepness of the concentration versus response curve. The relationship for N = 256 can be modeled with a standard Hill equation,
where En50 is the suppressor concentration associated with 50% maximum drug effect, Max is the maximum drug effect, and
. However, En50 is 0.025, an increase in potency of approximately 40 fold, and
Figure 3 shows two simple signal processing units, each with a site of action for drug A, and a site of action for drug B. These "sites of action" are either serial circuits comprised of many receptors for drugs that suppress signal propagation, "suppressors" as shown in Figure A2, or parallel circuits comprised of many receptors that enhance signal propagation, "enhancers" as shown in Figure A3. We will examine the behavior of four signal processing units, each comprised of two drugs: 1) Two suppressor drugs, whose sites of action are in a series circuit relative to each other, 2) Two suppressor drugs, whose sites of action are in a parallel circuit relative to each other, 3) Two enhancer drugs, whose sites of action are in a series circuit relative to each other, and 4) Two enhancer drugs, whose sites of action are in a parallel circuit relative to each other. An Excel spreadsheet, "circuits.xls," accompanies this manuscript in a supplement on the Anesthesia & Analgesia web site. Interested readers should download this spreadsheet and test the circuits to help understand their behavior. The mathematics were developed using Mathematica (Wolfram Research, Champaign, IL). Appendix D gives the derivation of the equations, which was performed using Mathematica. Appendix D and the Mathematica notebooks are available as Web supplements (www.anesthesia-analgesia.org) for interested readers. Series Signal Processing Unit: Two Suppressors in Series As explained in Appendix B, suppressors are likely to be arranged serially. The probability of a signal crossing N suppressors is
where Pun is the probability of passage past an unbound receptor, E is the efficacy of a single bound receptor suppressing signal propagation on a scale from 0 (no action) to 1 (complete suppression), and Su is the concentration of the suppressor as a fraction of the concentration associated with 50% receptor occupancy. Consider two different drugs, each of which suppresses signal propagation, arranged in series relative to each other. That would be the arrangement shown in the upper graph of Figure 3. The probability of passage is the probability that the signal would pass the assembly of receptors for suppressor A,
, times the probability the signal would pass the assembly of receptors for receptor B,
, namely
. Notice that the receptor for suppressor A and the receptor for suppressor B have unique probabilities for passage of a signal past an unbound receptor (Pun,A and Pun,B), as well unique efficacies (EA and EB). In our analysis we will consider how the probability of passage past an unbound receptor and the intrinsic efficacy of a single bound receptor affect whether drugs are additive or synergistic in their interaction. For simplicity we assume that Pun,A = Pun,B, and EA = EB. With the simplifying assumptions that Pun,A = Pun,B and EA = EB, suppressor A and suppressor B have identical behavior in this circuit, except that they need not exist in the same concentration. Given our simplifying assumptions, X units of suppressor A given alone has exactly the same effect as X units of suppressor B given alone. What happens if we give suppressors A and B concurrently? Let's call M the Multiplier on X, so that M x X units of suppressor A, given concurrently with M x X units of suppressor B produces exactly the same effect as X units of suppressor A alone or X units of suppressor B alone. If M = 0.5, then we have strict additivity, because 0.5 X units of A + 0.5 X units of B produces the same effect as 1.0 X of either drug alone. We have defined synergy as the sum of fractions being <90% of the normalized dose (i.e., X in this example) of either drug alone. If we can produce the target drug effect with M <0.45, then we have synergy, because 0.45 X units of drug A plus 0.45 X units of drug B gives 0.90 units of X, exactly at the border between synergy and additivity. The effect of a single suppressor in this circuit, given at concentration Su, is
given that Su of the other drug is 0, reducing equation 6 to PunN. The effect of two suppressors, each given at a concentration of M Su, is
We are looking for the value of M such that these two effects are equal. Thus, we want to find M so that
=
The solution for M is shown in Appendix D, which is available in the online supplementary material. Having solved for M, the fraction of two drugs that, given together, produces the same effect as either drug given alone, we can now examine the behavior of M as a function of:
If we could create graphs in five dimensional space we could capture the relationship in a single figure. Since this is not possible, Figure A4 explores the relationship as a series of three-dimensional graphs.
In every graph in Figure A4, the Z axis is M, the multiplier on the concentration, Su, of a single suppressor such that two different suppressors, each at a concentration M Su, produces the same effect as either suppressor alone at a concentration of Su. In every graph the X axis is the concentration of suppressor producing the desire effect by itself. The Y axis is Pun (top row), E (second row) or N (third and fourth rows). The most obvious relationship shown in A4 is that regardless of the value of Pun, N, and E, M approaches 0.5 (perfect additivity) for very small values of Su. Since Su = 1 is the concentration associated with 50% receptor occupancy, Figure A4 shows that regardless of the other parameters, two suppressor circuits, arranged in a signal processing unit in serial with respect to each other, will always show simple additivity if the drug effect of interest occurs at concentrations that are <1% of those associated with 50% receptor occupancy. This is not a surprising finding, as fractional receptor occupancy increases almost linearly with concentration at very low levels of receptor occupancy. On the other hand, if the drug effect of interest requires 50% receptor occupancy (Su = 1), then profound synergy is expected, regardless of the values of Pun, N, and E. Looking at more subtle findings for this circuit, the top row of graphs shows that relationship between M and Su is independent of Pun for this circuit. The second row of graphs shows that for values of Su greater than 0.1, E, the intrinsic efficacy, influences whether the relationship is additive (values of E near 1) or synergistic (values of E <0.5). However, there is almost no change in the relationship as E varies from 0.001 to 0.5. The last two rows of graphs explore the influence of N, the number of suppressor receptors in series at Site A, and at Site B (i.e., N = NA = NB). There is virtually no influence of N on whether the relationship demonstrates synergy or additivity. Thus when a signal processing unit consists of two suppressors arranged in series, the relationship between the two drugs can be synergistic or additive, depending on the concentration at which the drug exerts the effect of interest relative to the concentration associated with 50% receptor occupancy. Series Signal Processing Unit: Two Suppressors in Parallel Two suppressors in a parallel circuit have a more complex behavior. The effect of a single suppressor in this circuit, given at concentration Su, is
. The effect of two suppressors, each given at a concentration of M Su, is
. We are looking for the value of M such that these two effects are equal. Thus, we want to find M so that
= 1
. The solution for M is shown in Appendix D, which is available in the online supplementary material. Figure A5 shows the relationship between M, Su, Pun, N, and E. Overall, we see the same pattern as in Figure A4. The relationship is additive for values of Su that are <1% of the concentration associated with 50% receptor binding, regardless of the other parameters. At concentrations higher than 0.01, the relationship is more likely to be synergistic. As seen for suppressors in series, Pun has little effect. Interestingly, the effect of E is the opposite of what was seen in Figure A4, in that values of E nearly 1 are associated with profound synergy. There is an effect of N, but it is modest compared to the effect of Su. It is interesting that several of the graphs are incomplete. At large values of N these signal processing units produce very small numbers for probability of signal propagation, typically in the range of 10– 20. Such small numbers challenge the accuracy of the computation, particularly when the final calculation, M, involves the ratio of two very small numbers. The graphs become irregular, and then unsolvable, due to numeric instability when challenged by large values of N in combination with very low values of Pun and E.
Series Signal Processing Unit: Two Enhancers in Serial This was the most challenging circuit to analyze. Looking at the top circuit in Figure 3, one might guess that if Site A and Site B are assemblies of enhancers arranged in parallel, and the probability of passage was at each site was
, then the probability of passing both sites, each with a concentration of En was
. This is not correct. The reason is that at Site B the "incoming" probability is not one, but is the probability coming out of Site A, which is
. Every one of the parallel enhancer circuits sees this probability, not one, in the calculation for the probability effects of site B. This can be readily appreciated with the Excel spreadsheet "circuits.xls" provided with this manuscript as a supplement. An interesting aspect of this is that enhancer A and enhancer B are not exchangeable. Enhancer B is more potent than enhancer A. Enhancers A and B also have different peak effects and different half-maximal effects, even when Pun, N, and E are identical. As explained in Appendix D, this circuit can be analyzed as follows. First, propose an effect that occurs at a certain concentration of enhancer A. Because enhancer A is less potent than enhancer B, any effect produced by enhancer A can also be produced by enhancer B. The effect from EnA units of enhancer A is
(See appendix D in the online supplementary material). Second, find the concentration of enhancer B, EnB, that produces the same effect. A closed form solution exists and is given in Appendix D. This establishes the concentrations of enhancers A and B that produce a given drug effect. Third, find M, so that M · EnA units of enhancer A, given concurrently with M · EnB units of enhancer B, produces exactly the same drug effect as EnA units of enhancer A, or EnB units of enhancer B. There is no closed form solution for M. It must be solved numerically. Appendix D shows how the numeric solution was found with Mathematica. Figure A6 shows the relationship between M and EnA, Pun, N, and E. As was seen for suppressors arranged in serial and parallel signal processing units, the most important variable for two enhancers arranged in a serial signal processing unit is EnA, the concentration associated with the effect of interest. The relationship is almost always additive for values of EnA that are <1% of the concentration associated with 50% receptor binding, although the graph in the upper left corner shows that this is not always the case. Figure A6 also demonstrates that when two enhancers are arranged in a serial circuit, M decreases (and hence synergy increases) increases as Pun decreases, and E increases. The effects of N are modest, and not easily characterized (e.g., the first figure in the third row shows decreasing M at both low and high values of N).
Series Signal Processing Unit: Two Enhancers in Parallel The effect of a single enhancer in this circuit, given at concentration En, is
. The effect of two enhancers, each given at a concentration of M En, is
. We are looking for the value of M such that these two effects are equal. Thus, we want to find M so that 1–(1–Pun)N
. The solution for M is shown in Appendix D in the online supplementary material. Figure A7 shows the relationship between M and En, Pun, N, and E. Even though the function for two enhancers arranged in parallel is not the same as the function for two suppressors arranged in serial, Figure A7 is indistinguishable from Figure A4. As was true for the other circuits, the relationship is strictly additive if the effect of interest is observed at values of En that are <1% of the value associated with 50% receptor occupancy. The values of Pun and N have no effect on whether the interaction is additive or synergistic. Very high values of E are associated with a modest decrease in synergy. Summary Perhaps anticlimactically, all four simple circuits generated qualitatively similar results. When two suppressors or two enhancers are arranged into simple signal processing units, they demonstrate synergy when the drug effect of interest approximates or exceeds 50% receptor occupancy. However, if the drug effect occurs at concentrations much <50% receptor occupancy, then the drug effect is more likely to be additive.
Appendix D involves several very large equations that are not easily typeset. It is available in the online supplementary material at www.anesthesia-analgesia.org.
This article has supplementary material on the Web site: www.anesthesia-analgesia.org. Accepted for publication March 27, 2008. Dr. Steven L. Shafer, Editor-in-Chief, was recused from all editorial decision, related to this manuscript. Dr. Flood is the wife of Dr. Shafer, Editor-in-Chief of Anesthesia & Analgesia. This manuscript was handled by James Bovill, former Section Editor of Anesthetic Pharmacology and Dr. Shafer was not involved in any way with the editorial process or decision.
This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|