Introduction

Top Abstract Introduction Formulation of the Model Behavior of the Model Discussion Appendix References

The majority of the agonists and antagonists used to study pharmacological receptors (and for therapeutic intervention) are competitive, i.e., they bind to the same site as the endogenous ligand on the receptor, the orthosteric site. Allosteric modulators are molecules that exert their effects via a site on the receptor protein that is distinct from the orthosteric site. With the notable exception of muscarinic receptors (e.g., Lazareno et al., 1998), allosteric modulators have largely been ignored as targets for drug discovery programs. The allosteric effects of benzodiazepine receptor ligands at -aminobutyric acid_A receptors were discovered after their in vivo effects were known (Ehlert et al., 1983). There are arguments, however, that suggest that allosteric modulators may, under certain circumstances, provide better drugs than their orthosteric counterparts. In particular, the effects of allosteric modulators are saturable and thus there is less likelihood of adverse effects from overdose. This property can be capitalized upon to provide increased duration of action, by increasing dosage, while still avoiding unwanted effects. Also, allosteric activators of receptors may have advantages over direct agonist molecules. The effects of this type of compound require the presence of the endogenous agonist and are thus likely to mimic more closely the normal physiological effects of that agonist. Allosteric sites are also likely to be less well conserved than the ligand binding site in a particular receptor family thus potentially allowing the design of ligands with greater subtype selectivity (Tuek and Proka, 1995).

The muscarinic acetylcholine receptor family has been widely studied as a model system for the interaction of allosteric modulators of G-protein-coupled receptors (Stockton et al., 1983; Ellis et al., 1991; Lee and El-Fakahany, 1991; Lazareno and Birdsall, 1995; Proka and Tuek, 1995). In this case, a ternary complex model (Fig. 1a) has been proposed to describe the effects of these allosteric modulators on the binding of ligands to the orthosteric site (Stockton et al., 1983; Tuek and Proka, 1995). A number of adaptations of this basic model have also been proposed to explain more complex behaviors (Waelbroeck, 1994; Lazareno and Birdsall, 1995; Proka and Tuek, 1995; Hoare and Strange, 1996). However, there are, to my knowledge, no reports of models of allosteric modulation in which modern theoretical models of receptor activation are explicitly considered, although Ehlert's work is important in terms of the more classical models (Ehlert, 1988). The two-state model of receptor activation (originally presented in Karlin, 1967; Thron, 1972; and Colquhoun, 1973; and recapitulated by Leff, 1995) (Fig. 1b) is one of the simpler models of receptor activation that qualitatively describes our current understanding of receptor behavior. In particular, unlike earlier models (Stephenson, 1956; Furchgott, 1966; Black and Leff, 1983), it allows receptors to have a low level of activity in the absence of agonist. Another advantage of this model, for the present discussion, is that it models only the behavior of the receptor and makes no assumptions about subsequent signal transduction steps, thus making it potentially applicable to any type of receptor. This is in contrast with models such as the ternary complex (De Lean et al., 1980) and cubic ternary complex models (Weiss et al., 1996a-c), which are specifically formulated to describe the behavior of G-protein-coupled receptors.

View larger version (9K): [in this window] [in a new window]   Fig. 1.   a, the ternary complex model of allosteric modulation. The two ligands, A and B, bind to different sites on the same receptor, R. This results in two binary complexes AR and RB and the ternary complex ARB. The affinity of ligand for the complementary binary complex may be different from its affinity for the free receptor. The equilibrium constants are explained in the text and defined in Table 1. b, the two-state model of receptor activation. The receptor exists in two forms, an inactive state, R, and an active state, R*. The ligand, A, can bind to either of these forms, to generate AR or AR*, and may discriminate between them. Again, the equilibrium constants are defined in Table 1.

In this report, therefore, the principles of the two-state model of receptor activation have been combined with those of the ternary complex model of allosteric modulation to generate a general qualitative model of the interaction between allosteric and orthosteric ligands at pharmacological receptors both in terms of binding and of functional activation of the receptor. Models of this form have previously been published to aid discussion of the effects of allosteric modulators in specific receptor contexts (Bruns and Fergus, 1990; Ehlert et al., 1983), however, the mathematical analyses were not presented and the possibility of cooperativity at the level of receptor activation was not considered.

	Formulation of the Model

Top Abstract Introduction Formulation of the Model Behavior of the Model Discussion Appendix References

In the ternary complex model (Fig. 1a), the receptor (R) can bind reversibly to an (arbitrarily assigned) orthosteric ligand (A) with affinity constant K, or to an allosteric ligand (B) with affinity constant M, resulting in two binary complexes AR and RB (the order denoting the binding site). The unbound ligand may subsequently bind to the complementary binary complex (e.g., A to RB), its affinity being modified by the allosteric constant . This parameter is therefore analogous to the cooperativity factor designated  by Ehlert (1988) and Tuek and Proka (1995). In this formulation of the model, a value of  greater than unity indicates positive cooperativity and a value of  less than unity indicates negative cooperativity between the two ligands. When  = 1 there is no interaction between the two compounds.

In the two-state model (Fig. 1b), the inactive receptor (R) enters the active state by undergoing a conformational change to form the active receptor (R*). Functional responses are assumed to be proportional to the amount of R*-containing species within the system (Leff, 1995), hence, strictly, it is the behavior of the pharmacological stimulus rather than the response that is determined in this model. In the absence of ligand, the proportion of receptors in the R* state is governed by the equilibrium constant L. The ligand (A) binds to the inactive state of the receptor with an affinity constant, K, and to the active state of the receptor with a modified affinity, K. Binding of the ligand shifts the isomerization equilibrium constant to L. In this model the allosteric constant  can be regarded as the intrinsic efficacy of the ligand. If A has a value of  greater than unity it favors the formation of R* and is an agonist. A value of  equal to unity indicates a neutral antagonist (i.e., no receptor state selectivity) and a value of  less than unity indicates a molecule that stabilizes the R state, i.e., an inverse agonist.

As shown in Fig. 1, both of these models can be represented by cyclic equilibria and appear to contain the common step of a ligand binding to the endogenous agonist binding site (A binding to R). However, the equilibrium constants governing these events differ in character between the two models. In the ternary complex model, K represents a macroscopic (and therefore experimentally measurable) equilibrium constant for the binding of the orthosteric ligand to the receptor. In contrast, in the two-state model, K represents the microscopic binding affinity of the ligand for the inactive state of the receptor, a quantity that cannot be determined experimentally [the measured affinity constant of A in the two-state model is given by K(1 + L)/(1 + L)]. The approach inherent in the ternary complex model can, however, be applied to the two-state model on the basis of an effect of the allosteric ligand on the affinity of the orthosteric ligand for the inactive receptor.

Inserting an allosteric ligand into the two-state model and completing the equilibrium to allow for all of the potential interactions leads to the cubic model shown in Fig. 2 (a ternary complex model and the two-state model form, respectively, the front and left-hand faces of this new model). In the combined model, the allosteric two-state model, the allosteric modulator could be an agonist or an inverse agonist in its own right. This extra property of the allosteric modulator is denoted by the intrinsic efficacy term  and forms the upper face of the cube. The allosteric interaction of the two ligands at the level of binding (to the inactive receptor) is again described by . However, it is also possible for the binding of one of the ligands to modify the ability of the other ligand to activate the receptor (i.e., change its intrinsic efficacy). This potential activation cooperativity is modeled by the parameter . The parameters of the model and their definitions are summarized in Table 1. It should be noted that, although it is based on the two-state model, this model allows the two ligands to bind to the binary complexes with affinities that are different from their affinities for R and R*. Strictly, therefore, it should be described as a multistate rather than a two-state model. This also brings the resulting model to a level of analysis that is intermediate between that of the two-state and ternary complex models. Its equilibrium constants are neither macroscopic because they are not measurable nor microscopic because conformational changes in the receptor are implicit in nonunit values of the allosteric constants  and  (e.g., when   1 the conformation of the receptor in AR must be different from that of R otherwise B would not distinguish between them). The cubic ternary complex model (Weiss et al., 1996a) is also formulated at this level of analysis.

View larger version (23K): [in this window] [in a new window]   Fig. 2.   The allosteric two state model. The free receptor exists in inactive and active forms, R and R*, as in the two-state model. The two ligands, A and B, bind to different sites on the receptor and may discriminate between the two states and between the free receptor and the binary complexes. In the model, A is assumed to bind to the orthosteric site and B to the allosteric site; however, this assignment is entirely arbitrary due to the symmetry of the model. The equilibrium constants are defined in Table 1.

	Behavior of the Model

Top Abstract Introduction Formulation of the Model Behavior of the Model Discussion Appendix References

The derivation of expressions for the binding of the orthosteric ligand to the receptor in the presence of an allosteric ligand and for the functional activation of the receptor are shown under Appendix, as are expressions for the midpoints of the saturation binding, competition binding, and activation curves and for the maxima and minima of these curves. These expressions were used to investigate the effects of the various parameters of the model on binding and activation curves. The effects of varying parameters , , or  in turn were investigated in simulations of binding and functional experiments. These parameters (and ) may take any nonnegative value; a value of unity implies the parameter has no effect, i.e., it removes any selectivity between the states that parameter governs. The magnitudes of values above and below unity should be considered in log terms, thus the converse of  = 100 is  = 0.01. In the simulations, the two allosteric constants that were not being varied were set to unity to remove their effects from the system.

Simulations of the effects of varying the two affinity constants (K and M) are not shown because they are entirely as would be expected from a change in the affinity of the radioligand or the allosteric modulator. Similarly, varying  or L have effects analogous to those of the original two-state model (Leff, 1995) and are not discussed. Obviously, because  is the intrinsic efficacy of the orthosteric ligand, it was necessary to vary this parameter when simulations required an orthosteric agonist or antagonist ligand. The value of L in the simulations shown below was chosen to produce very little activity in the absence of agonist (unless otherwise stated), because this most closely mimics the behavior of most experimental systems.

Simulations of Competition Binding Curves. In these simulations the concentration of the orthosteric ligand, which adopts the role of radioligand, was chosen to provide a control fractional occupancy of ~0.5, i.e., at the apparent K_D. This concentration is most frequently used in radioligand binding experiments and allows effects of the allosteric ligand in either direction to be seen clearly.

Effect of  (Intrinsic Efficacy of B). The effects of this parameter on binding are dependent on the pharmacological activity of the radioligand that is used (Fig. 3). With an agonist radioligand ( > 1), binding is increased with increasing concentrations of an allosteric agonist ( > 1) and is decreased if B is an inverse agonist ( < 1) (Fig. 3a). When B is neutral ( = 1), it has no effect on the binding of A. These effects can be explained in terms of the effect of B on the conformation of the receptor. An allosteric agonist will increase the level of R* in the system and make it energetically more favorable for an orthosteric agonist to convert more receptors into this state. R* is the state of the receptor with a high affinity for agonists and an increase in the concentration of this form of the receptor will therefore increase the amount of binding of an orthosteric agonist due to an increase in its apparent affinity (Fig. 3b).

View larger version (21K): [in this window] [in a new window]   Fig. 3.   Effects of parameter  (the intrinsic efficacy of B) on binding curves. a, effect of varying  on competition binding curves against an agonist radioligand ( = 1000). Other parameters were [A] = 0.1/K, K = M = 1, L = 0.01,  = 1000,  =  = 1, and  = 100 (), 30 (), 10 (), 3 (), 1 (*), 0.3 (), 0.1 (), 0.03 (), or 0.01 (). b, effect of various concentrations of an allosteric agonist ( = 100) on saturation binding curves (presented on a log scale) to an agonist radioligand ( = 1000). Note the increase in apparent affinity of the radioligand. Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.01 (), 0.03 (]), 0.1 (), 0.3 (), 1 (), and 3/M (). c, effect of  on competition binding curves against an orthosteric inverse agonist radioligand ( = 0.001). Other parameters are as described in a.

Effect of  (Binding Cooperativity). The effects of this parameter are independent of the pharmacological activity of the radioligand used (Fig. 4a). Thus, positive cooperativity ( > 1) results in an increase in the level of binding of the radioligand and negative cooperativity ( < 1) results in a decrease. These effects are again due to changes in the apparent affinity of the receptor for the orthosteric ligand caused by the allosteric ligand (Fig. 4b). This is of course exactly the behavior of the analogous parameter in the ternary complex model of allosterism (Tuek and Proka, 1995). It should be noted that negative allosteric modulators do not inhibit binding to the nonspecific level unless the negative cooperativity is very strong, that is if   1 (the maximal reduction of binding will of course depend on the concentration of the radioligand).

View larger version (16K): [in this window] [in a new window]   Fig. 4.   Effects of parameter  (the binding cooperativity between A and B) on binding curves. a, effect of varying  on competition binding curves against a neutral antagonist radioligand ( = 1) (for illustration, the effects of this parameter are independent of ). Other parameters were [A] = 1/K, K = M = 1, L = 0.01,  =  = 1, and  = 10 (), 3 (), 1 (), 0.3 (), 0.1 (), 0.03 (), or 0.01 (). b, effect of various concentrations of a positive allosteric modulator ( = 30) on saturation binding curves (presented on a log scale) to an antagonist radioligand ( = 1). Note the increase in affinity of the radioligand. Parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.03 (), 0.1 (), 0.3 (), 1 (), 3 (), and 10/M ().

Effect of  (Activation Cooperativity). The qualitative effects of this parameter are again independent of the pharmacological activity of the radioligand; however, the magnitude of the effect is affected. Thus, positive activation cooperativity ( > 1) results in an increase in the level of binding of a given ligand, whereas negative activation cooperativity ( < 1) can decrease it (Fig. 5). The magnitude of the effect of this parameter varies with , however. Thus, highly efficacious agonist ligands show very marked changes in their level of binding (Fig. 5a), whereas the effects are less pronounced as  decreases and only extremely large values of  change the binding characteristics of antagonists and inverse agonists (Fig. 5, b and c). As with , inhibitory effects seem more sensitive to  than potentiative effects.

View larger version (19K): [in this window] [in a new window]   Fig. 5.   Effect of parameter  (the activation cooperativity between A and B) on binding curves. a, effect of varying  on competition binding curves against an agonist radioligand ( = 10 000; [A] = 0.01/K). Other parameter were K = M = 1, L = 0.01,  =  = 1, and  = 100 (), 30 (), 10 (), 3 (), 1 (*), 0.3 (), 0.1 (), 0.03 (), or 0.01 (). b, effect of varying  on competition binding curves against a neutral antagonist radioligand ( = 1; [A] = 1/K). Other parameters were as described in a. c, effect of varying  on competition binding curves against an inverse agonist radioligand ( = 0.1; [A] = 1/K). Other parameters were as described in a.

Simulations of Activation Curves. In these simulations, the effects of various concentrations of the allosteric modulator were determined on concentration-response curves to an orthosteric agonist. This mimics the experimental protocol that would be used when intending to perform a Schild analysis. The effects of the different parameters can be distinguished by their effects on the asymptotes and midpoints of the concentration-response curves (CRCs).

Effect of . The parameter  is an intrinsic efficacy term in the allosteric two-state model and it behaves in the same way as  in terms of the activation state of the receptor. Thus, when  > 1, the allosteric ligand is an agonist and causes activation of the receptor in its own right. When  < 1, the allosteric ligand is an inverse agonist, it therefore inhibits any constitutive activity in the system. These effects can be seen as changes in the lower asymptotes of the CRCs to A in Fig. 6.

View larger version (19K): [in this window] [in a new window]   Fig. 6.   Effect of parameter  on activation curves to A. a, effects of various concentrations of an allosteric agonist ( = 100) on CRCs to A ( = 1000). Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] =0 (), 0.01 (), 0.03 (), 0.1 (), 0.3 (*), 1 (), 3 (), 10 (), or 30/M (). b, effects of various concentrations of an allosteric inverse agonist ( = 0.01) on CRCs to A ( = 1000). Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.1 (), 0.3 (), 1 (), 3 (*), 10 (), 30 (), 100 (), or 300/M ().

Effect of . The only effect of  is on the position of the CRCs; it has no effects on the upper or lower asymptotes of the curves (Fig. 7). Thus, values of  > 1 result in left-shifting of the curves (potentiation), whereas values of  < 1 cause right-shifting. The lack of effect of this parameter on the asymptotes is a logical consequence of its function in the model. It represents the binding cooperativity of the two ligands. However, it doesn't affect the ability of the ligands to discriminate between R and R*, thus it does not affect the maximal extent of functional activation. Again, this is borne out mathematically under Appendix where  is present in the denominator of the midpoint function for A (eq. 4). Thus, high values of this parameter decrease [A]₅₀ (increase potency), whereas values less than unity increase it. However, the limiting value of the upper asymptote of the CRC to A is L/(1 + L), which does not contain  and is therefore unaffected by it.

View larger version (17K): [in this window] [in a new window]   Fig. 7.   Effect of parameter  on activation curves to A. a, effects of various concentrations of a positive allosteric modulator ( = 100) on CRCs to A ( = 1000). Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.01 (), 0.03 (), 0.1 (), 0.3 (*), 1 (), 3 (), 10 (), or 30/M (). b, effect of various concentrations of a negative allosteric modulator ( = 0.01) on CRCs to A ( = 1000). Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.1 (), 0.3 (), 1 (), 3 (*), 10 (), 30 (), 100 (), or 300/M ().

Effect of . This parameter allows the allosteric ligand to modify the intrinsic efficacy of the orthosteric ligand without having any intrinsic efficacy of its own. In other words the modulator is not an agonist (or inverse agonist) in its own right and therefore has no effect on the lower asymptote of the CRC to A. When  > 1 the effective intrinsic efficacy of the orthosteric ligand is increased and the effects of A are potentiated and augmented (Fig. 8a). When  < 1, the effective intrinsic efficacy of the orthosteric ligand is decreased and the curves are right-shifted and the maximal effect is depressed (Fig. 8b). In common with binding, these effects are related to the effect of  on the ability of a ligand to convert RB to R*B. Thus, when  > 1 an agonist will convert RB to R*B more efficiently, R*B is an active form of the receptor thus a greater level of activation will result. The converse is true when  < 1. 

View larger version (20K): [in this window] [in a new window]   Fig. 8.   Effect of parameter  on activation curves to A. a, effects of various concentrations of an allosteric activator ( = 100) on CRCs to A ( = 1000). Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.003 (), 0.01 (), 0.03 (), 0.1 (*), 0.3 (), 1 (), 3 (), or 10/M (). b, effect of various concentrations of an allosteric inhibitor ( = 0.01) on CRCs to A ( = 1000). Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.1 (), 0.3 (), 1 (), 3 (*), 10 (), 30 (), 100 (), or 300/M ().

View larger version (18K): [in this window] [in a new window]   Fig. 9.   Effect of various concentrations of an allosteric activator ( = 100) on activation curves to a neutral antagonist ( = 1). Other parameters were K = M = 1, L = 0.01,  =  = 1, and [B] = 0 (), 0.003 (), 0.01 (), 0.03 (), 0.1 (×), 0.3 (), 1 (), 3 (), or 10/M ().

	Discussion

Top Abstract Introduction Formulation of the Model Behavior of the Model Discussion Appendix References

The model described in this report is a simple and logical extension of the two-state model when considering allosteric ligands. As in the two state model (and indeed other more recent models), this model highlights the inter-relatedness of the pharmacological concepts of affinity and efficacy (Leff, 1995; Weiss et al., 1996a-c). It predicts that the measurable affinity of a pharmacological agent depends not only on its binding affinity for the receptor but also on its ability to induce conformational changes within the receptor and the resistance of the receptor to those conformational changes. The model is, of course, also a modification of the ternary complex model of allosteric interactions. In this case, the model provides a framework in which each of the receptor species in the ternary complex model can cause downstream functional effects with different efficacies.

The behavior of the model is intuitively reasonable and consistent with published experimental observations (Bruns and Fergus, 1990; Lazareno and Birdsall, 1995; Kollias-Barker et al., 1997; Litschig et al., 1999). It predicts that an allosteric ligand that only affects orthosteric ligand binding (  1,  =  = 1) will only effect the midpoint of an orthosteric agonist with no effect on the maximal response. Only allosteric agents that can affect receptor activation, either through their own intrinsic efficacy (  1) or through cooperative modification of the orthosteric ligand's intrinsic efficacy (  1), can affect the maximal response induced by an orthosteric ligand. This potential cooperativity at the level of receptor activation (parameter ) also predicts the phenomenon of coagonism (Corsi et al., 1996). As described under Applicability of the Model the distinct behavior of each parameter in the model can be used to suggest potential mechanisms for the action of allosteric compounds.

An important practical prediction of this model is that it is only possible to generate a system-independent measure of the cooperativity between allosteric and orthosteric ligands when certain conditions are met (under Appendix). Specifically, when using functional assays to characterize an allosteric ligand, the allosteric ligand must not have any agonist (or inverse agonist) activity in its own right or its effects become system-dependent. Agonism is of course easy to test for in a functional assay. Inverse agonism is more difficult to determine because not all experimental systems are constitutively active. Fortunately, the inaccuracies introduced by inverse agonist activity are likely to be very small. In eq. 12, under Appendix, when  < 1, the terms in L (L is also less than 1) approximate to unity and the expression reduces to , which is still an independent characteristic of the ligand pair.

When using a ligand binding assay to characterize allosteric interactions, the requirements for system-independent quantification are even more stringent. The allosteric ligand must not only be devoid of agonist activity ( = 1) but must also have no effects on the efficacy of the orthosteric ligand ( = 1). Functional assays would, therefore, still be required even when characterizing an allosteric agent in binding assays (to allow appropriate interpretation of the data generated) because ligand binding assays cannot give meaningful information about agonist efficacy. If characterization were to be performed in ligand binding assays, the endogenous agonist would obviously be the preferred radioligand if available (and of course suitable).

Effects of an allosteric ligand on an orthosteric agonist's efficacy would most likely be manifest as a change in the maximal response that it elicits in a functional assay. However, care must be exercised when testing effects on maximal responses because the limiting factor to response size in many experimental systems is not the activation of the receptor but the activation of the signaling cascade (a manifestation of receptor reserve). This would make it impossible to detect increases in the maximal response to a full agonist and may truncate increases in the maximal response to partial agonists. It would also result in unreliable data on depressions of the maximal response. It may be possible to diagnose these problems by comparing the effects of a modulator on the midpoint of the agonist CRC and on its maximal effect. Discrepancies between the midpoints of the effects of the allosteric agent on these two parameters would indicate that the data could not be interpreted reliably. The most suitable assay system may therefore be one with very low receptor expression where even highly efficacious agonists are unable to fully activate the signal transduction cascade (see the concluding paragraph under Appendix). When an allosteric ligand does not meet the above-mentioned criteria (i.e., when it is also an agonist), the least complicated parameter that could be used to quantify its effect seems to be the ratio of the asymptotes of the maximal effect ratios and the dose ratios (under Appendix) taken from a functional assay.

Comparison with the Cubic Ternary Complex Model. The model described in this report has very close parallels with the cubic ternary complex model proposed by Weiss et al. (1996a-c). Indeed, if B were a G-protein the scheme in Fig. 2 would be identical with that used to describe the cubic ternary complex model. There is, however, one important difference between these two models in terms of the behavior of the allosteric constants and the activation curves. This difference is due to the way the two models quantify functional activity. In the allosteric two-state model, as in the two-state model, functional activity is quantified as the stimulus, i.e., the proportion of the receptor complexes that contain R*. In the cubic ternary complex model, the functional response is assumed to be proportional to the concentration of activated receptors that are bound to G-protein (Weiss et al., 1996a). This would be equivalent to [R*B] + [AR*B] in Fig. 2. This difference radically changes the effect of parameter . In the present model, the only effect of  is on the affinity of the orthosteric ligand; however, in the cubic ternary complex model,  modifies the ability of the receptor to interact with the G-protein and is therefore able to modify the maximal response induced by the ligand. Thus, in the cubic ternary complex model,  is an efficacy defining parameter along with , , and , whereas in the allosteric two-state model it is not. Another difference is that, in the cubic ternary complex model, it is possible to define a simple condition for agonism ( > (1 + L + M[B]/(1 + L + M[B])) whereas it is not meaningful to derive an equivalent expression in the allosteric two-state model.

Applicability of the Model. The acid test of a mathematical model is that it can be applied, at least qualitatively, to real, experimental data. A number of reports have been published on binding interactions between allosteric and orthosteric ligands (Ellis et al., 1991; Lee and El-Fakahany, 1991; Lazareno and Birdsall, 1995; Proka and Tuek, 1995; Hoare and Strange, 1996; Leppik et al., 1998). Because it is based on the most basic equilibrium model of these binding interactions, the allosteric two-state model is also compatible with these data, assuming that the allosteric and orthosteric ligands are site specific, of course. In terms of function, the prediction of coagonism by the allosteric two-state model has already been mentioned.

View larger version (12K): [in this window] [in a new window]   Fig. 10.   a, simulations of the effect of CPCCOEt on glutamate-induced inositol phosphate production. Parameters were defined such that the allosteric ligand decreased the maximal response to the agonist to 25% of the control level without altering the midpoint of the CRCs. The parameters were K = M = 1, L = 0.001,  = 10 000,  = 1,  = 8.5, and  = 0.03. The concentrations of B were 0 (), 0.3 (), 1 (), 3 (), 10 (), 30 (), 100 (), or 300/M (). b, simulated allosteric ligand competition binding curve using the parameters defined in a. The concentration of A was 0.1/K, however, there was no displacement of A over this range of concentrations of B at any concentration of A simulated (1 × 105  1000/K).

Concluding Remarks. This report describes a mathematical model of the interaction between orthosteric and allosteric ligands at pharmacological receptors that explicitly includes an effect of the allosteric modulator on the activation of the receptor. To my knowledge this is the first detailed report of such a system, particularly in its application to G-protein-coupled receptors (with the caveat of the lack of G-protein already mentioned). The behavior of this model is qualitatively consistent with the behavior of allosteric modulators in radioligand binding assays. Importantly, this model can also simulate and therefore suggest a potential mechanism for the effects of an allosteric modulator of metabotropic glutamate receptors that inhibits functional activity but does not affect the binding of the agonist. It also suggests a number of caveats that should be considered when analyzing the interaction between allosteric and orthosteric ligands at pharmacological receptors.