Drug-Acceptor Interactions, Drug-Acceptor Interactions

##article.fontSize##:  ##icon.font_small.alt##  ##icon.font_medium.alt##  ##icon.font_large.alt##
Drug-Acceptor Interactions

Chapter 7: Cubic Reaction Schemes. ATSM and HOTSM

Some general considerations about cubic models are described in sub-chapter 7.1, followed by the derivation and application of some cubic models, especially the allosteric two-state model (ATSM) and the homotropic two-state model (HOTSM).

7.1 Introduction to Cubic Models

7.1.1. Two Simple Cubic Models

We have now advanced to two two-state models that are essential for our conception of mechanisms in modulatory regulation at the receptor level. In both models, two ligands can interact via a receptive unit and the receptor may exist unliganded in two isomeric forms, a reactive or an active conformation. The two models are the allosteric two-state model (ATSM) and the homotropic two-state model (HOTSM) (Fig. 7.1). Two different (heterotropic) ligands act in the ATSM, while two identical (homotropic) ligands are active in the HOTSM. One could argue that a better name for the allosteric two-state model would be the ‘heterotropic two-state model’ (HETSM), but currently I will adhere to the name settled in the literature – the allosteric two-state model (ATSM).

Citation: Bindslev N 2008. DOI: 10.3402/bindslev.2008.10

© 2008 N Bindslev. This book and all matter and items published therein are distributed under the terms of the Creative Commons Attribution-Noncommercial 3.0 Unported License (http://creativecommons.org/licenses/by-nc/3.0/), permiting all non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Figure 7.1. (A) The allosteric two-state reaction scheme (ATSM). (B) The homotropic two-state reaction scheme (HOTSM). Arrows indicate the flow of information or path of activation. Receptor symbols and parameters are listed and explained in Table 7.1 and Box 7.1.

The ATSM was first analyzed by Hall (2000) varying either the agonist or the modulator concentration separately. Although the concentration of the heterotropic ligands in ATSM may be varied simultaneously in a fixed ratio (see later), this possibility was not explored by Hall (2000).

The HOTSM was developed by Bindslev (2004). In this model, the agonist and the modulator concentration by necessity vary simultaneously, since agonists and modulators will be identical molecules.

As we shall see, when a mixture of agonist and modulator molecules with a fixed concentration ratio are employed for dose-responses in the functional form of ATSM, this model behaves as a HOTSM.

Both models are explored in more detail in the two papers mentioned above, and should be consulted for an in-depth understanding of the model characteristics and capabilities as tools for synagic analysis. Here I shall present some highlights and additional possibilities for the two models.

There are three principal conditions in both the ATSM and the HOTSM: (1) binding of one ligand, (2) binding of two ligands simultaneously, and (3) a conformational isomerization for activation. Reaction schemes to describe the possible combinations of these three conditions result in cubic networks with seven independent system constants, as shown in Fig. 1A + B.

The ATSM is a cubic combination of the intervention model (Chapter 2) and the cyclic two-state model (cTSM) (Chapter 5). The HOTSM is a cubic combination of the auto-intervention model (Chapter 3) and again the cTSM described in Chapter 5.

The ATSM and HOTSM are both parts of and thus included in the four-pane two-state model (FP-TSM) in Fig. 7.2. The designation of system constants in model-pair ATSM and HOTSM in Fig. 7.1 is transferred from the FP-TSM in Fig. 7.2 in order to maintain a consistency of terms for the ATSM and HOTSM with those of the mother-model, FP-TSM. Terms for the receptor complexes and the parameters for ATSM and HOTSM are listed in Table 7.1 and Box 7.1. Apostrophes and subscripts at certain parameter symbols secure the difference between the two models (compare panels A and B in Fig. 7.1 and Fig. 7.2).

Figure 7.2. ATSM and HOTSM as parts of the four-pane two-state model (FP-TSM). (A) Without system constants. (B) With system constants.

Table 7.1. Ligands, receptor conformations and parameters in ATSM and HOTSM
*The system constants are all forward constants.

Compare Table 7.1 with Box 7.1.

Box 7.1. Receptor conformations and parameters in ATSM and HOTSM

ATSM

Acronyms are as listed here for the eight receptor conformations together with symbols for the seven independent parameters of ATSM in Fig. 7.1A:

Receptor conformations in ATSM

R unbound reactive receptor
RS agonist-bound reactive receptor
R* unbound active receptor
R*S agonist-bound active receptor
MR modulator-bound reactive receptor
MR* modulator-bound isomerized receptor with activity
MRS agonist- and modulator-bound reactive receptor
MR*S agonist- and modulator-bound receptor with activity

Parameters in ATSM

L isomerization constant–efficacy constant
A ss equilibrium association constant at the primary site for agonist S
A mm equilibrium association constant at the secondary site for modulator M
a intrinsic efficacy constant with an agonist already bound
b intrinsic efficacy constant with a modulator already bound
c co-operativity coefficient for binding of a second ligand when a first ligand is already bound
d co-operativity coefficient for binding of a second ligand when a first ligand is already bound to an active receptor or an intrinsic efficacy constant when two ligands are bound

HOTSM

Acronyms are as listed here for the eight receptor conformations together with symbols for the seven independent parameters of HOTSM in Fig. 7.1B:

Receptor conformations in HOTSM

R unbound reactive receptor
RS substrate-bound reactive receptor
R* unbound active receptor
R*S substrate-bound active receptor
SR modulator-bound reactive receptor without activity
SR* modulator-bound isomerized receptor with or without activity
SRS substrate- and modulator-bound reactive receptor
SR*S substrate- and modulator-bound receptor with or without activity

Parameters in HOTSM

L isomerization constant–efficacy constant
A ss equilibrium association constant at the primary site for substrate S
A ms equilibrium association constant at the secondary site for substrate S
a intrinsic efficacy constant with a substrate already bound at ‘O’
b’ intrinsic efficacy constant with a substrate already bound at ‘M’
c co-operativity coefficient for binding of a second ligand when a first ligand is already bound
d’ co-operativity coefficient for binding of a second ligand when a first ligand is already bound to either site of an active receptor or an intrinsic efficacy constant when two ligands are bound

‘O’ stands for the orthosteric or primary binding site, and ‘M’ designates the modulatory or secondary binding site.All parameters are forward constants for both ATSM and HOTSM.

The ATSM and the HOTSM are close siblings. In fact, it may be said that the HOTSM is a special model variant of the ATSM, since from a functional aspect, as already indicated, for ATSM with mixes of the two different ligands, when their ratio of concentration is kept constant during application, the ATSM becomes a copy of the HOTSM. This aspect of the ATSM/HOTSM pair was analyzed by Bindslev (2004), and discussed further in Section 7.8.1. However, formulation of occupancy in the two models results in completely different expressions (see Section 7.8.2).

7.1.2. A Promising Model for the Analysis of Modulator Action

Hall's ATSM is a clear choice for analyzing modulator effects and will probably be exploited in the coming years. ATSM is a modification of the ternary complex model of ‘allosteric’ interaction by Ehlert (1988). In the words of Hall (2000), the ATSM ‘provides a framework in which each of the receptor species in the ternary complex model can cause downstream functional effects with different efficacies’. In a fixed concentration ratio mode, both ATSM and HOTSM have three plateaus of effect, and the models predict that an allosteric ligand which affect the affinity but not the efficacy of a primary ligand will only influence a second plateau of the model's dose-response curves, and neither its third plateau at high ligand concentrations nor its initial first plateau before ligand application (see later in Fig. 7.27). It is only modulator ligands with an efficacy constants L different from unity that can affect the third level of function induced by a primary ligand.

Thus far, only sporadic implementation of the ATSM has been reported (see, e.g., May et al. 2004; Jensen & Spalding 2004; Franco et al. 2006; Langmead & Christopoulos 2006; Ehlert & Griffin 2008; and partially by Hoare et al. 2008). However, with the accelerated development of allosteric drugs and their clinical application (Table 7.2) (Gao & Jacobson 2006, sub-chapter 7.12), the relevance of ATSM as an analytical tool justifies a more extended examination of this model. Therefore, in Chapter 7, I will scrutinize the ATSM and its sibling the HOTSM for some general principles and some details.

Table 7.2. Examples of receptors for which there is development of allosteric drugs for better therapy (see also Raddatz et al. 2007)
The table does not include allosteric compounds for enzymes or growth factor receptors.

*See also Baggio & Drucker (2007), especially pp. 2146–2149.

*Small peptide ago-allosteric modulator for a GPCR family B member.

7.1.3. Other Cubic Models and Other Models

Other reaction schemes have been presented which also result in cubic models when based on the three principal conditions mentioned above. For instance, for G protein coupled receptors (GPCRs), a cubic ternary-complex model (CTCM), was scrutinized by Weiss and co-workers (Weiss et al. 1996a b c). On its development see Fig. 7.3. In the CTCM it is merely the activated receptor complexes when coupled to G protein (GR*) that are considered functional in the transduction pathway (arrows in Fig. 7.4A). Nevertheless, even in the absence of agonists but with functionality linked to G proteins coupled to the receptor, the CTCM includes the possibilities for simulation of both spontaneous activity and inverse agonism, as inherent in ATSM and HOTSM.

Figure 7.3. Historic development of the cubic ternary-complex model (CTCM). The ternary complex model (TCM) (de Lean et al. 1980), and the extended ternary-complex model (ETCM) (Samama et al. 1993) appear on the road towards the CTCM. Taken from Weiss et al. (1996), p. 160, Fig. 7.12) with permission.

Figure 7.4. The cubic ternary-complex model (CTCM). (A) Arrows indicate the flow of information or path of activation. The symbol R stands for G protein coupled receptor (GPCR), G for G Protein and parameter g replaces parameter b of the ATSM. (B) The lower panel shows the extended ternary complex model (ETCM), which is part of the CTCM in (A). The allosteric two state model is included to demonstrate that ECTM can be derived from ATSM as well.

To date, the CTCM has been implemented by few research groups (e.g., Bruheim et al. 2003; Monczor et al. 2003; Fitzsimmons et al. 2004; O'Brien et al. 2004; Pineyro et al. 2005).

Adding another ligand to the CTCM scheme, that is, operating with three different ligands, but baring simultaneity in binding, i.e., excluding tri-ligand binding, was also considered by Weiss and co-workers (1996a) and visualized in their Fig. 5. Other variations without explicit conformational isomerization have been suggested on the theme of cubic complex modeling, but with three ligands interacting by simultaneous binding (Marvizon & Baudry 1993; Kukkonen et al. 2001).

A model more complete than the CTCM for GPCRs was suggested by Christopoulos et al. (1998), and again by Christopoulos & Kenakin (2002), as before, requesting G protein coupling for activity and including spontaneous activity, inverse agonism, but now with tri-ligand binding, the same as simultaneous binding of three different ligands including (1) an agonist, (2) the G protein, and (3) a modulator molecule. It is a so-called cubic quaternary-complex model (CQCM) (Fig. 7.5A), with four principal conditions and still, as mentioned, implicating both inverse agonism and spontaneous activity in the presence of bound G protein. The CQCM may alternatively be represented as a cube-in-a-cube (Fig 7.5B). The cube-in-a-cube was discussed in a different setting as a ‘hyper-cube’ model with four ligands, that is, replacing the possibility of conformational isomerization by one extra ligand binding, leading to a pentanary-complex model (van Rijn & Willems-van Bree 2004).

Figure 7.5. The cubic quaternary-complex model (CQCM). (A) From Christopoulos et al. (1998) with permission. (B) In a double cubic version, cube-in-a-cube (courtesy of CHETAN® 2003). Both versions of the model assume the possibility of three ligands bound simultaneously, excluding co-lateral binding, and a conformational switch between a reactive and an active state. The model has 16 receptor conformations and 15 independent parameters.

The CTCM and its offspring, the CQCM, are further characterized in sub-chapter 7.11.

7.1.4. Contents of Chapter 7

As usual, both the ATSM and the HOTSM are analyzed for their characteristics related to binding studies (sub-chapters 7.2 and 7.5), as well as functional experimentation (sub-chapters 7.3 and 7.6). Further aspects on HOTSM are in sub-chapter 7.7. Overlap between ATSM and HOTSM is considered and visualized in sub-chapter 7.8.

A combination and extension of ATSM and HOTSM into the FP-TSM (Fig. 7.2), and this model's possibilities for simulating synagics are considered in sub-chapter 7.9. A subject related to the FP-TSM, the so-called OFCOR principle, is looked into in sub-chapter 7.10. Differences and similarities between ATSM/HOTSM and other cubic models are debated in sub-chapter 7.11. Finally, sub-chapter 7.12 exposes how a 10-year old market for allosteric drugs in therapy is longing for the implementation of a two-state model à la Hall's ATSM. For an overview, together with other models, some cubic models are listed in Table 7.3.

Table 7.3. Some examples of one-, two- (dimer), or multi-sited models in one-, two-, or three states involving one-, two-, or more ligands-and some in cubic representations (See also, for example Casadó et al. 2007)
*Cross-talk between binding sites.

7.2. Aspects for Binding with the ATSM

7.2.1. Basics of the Model

The ATSM by Hall (2000) is shown in Fig. 7.1A. Receptor conformations and parameters for the ATSM are given in Table 7.1 and Box 7.1. A distribution equation for the binding aspect of ATSM with ao as actual occupancy is:

Compare Eq. 7.1 with the corrected Hall Eq. 3 (Hall 2000, Erratum in: Mol Pharmacol 59: 161, 2001). In Hall's analysis it is the concentration of either an agonist or a modulator that is varied as the independent variable.

7.2.2. Variation in Single Parameters

Variations in single model parameters for the binding aspect of ATSM are shown in Fig. 7.6 for five different values of each parameter while increasing the agonist concentration.

Figure 7.6. Examples of parameter changes in ATSM for binding with the agonist [S] as independent variable. (A) In panels A–F, parameters b, c, d, L, A ss, and A mm, as indicated, vary in five steps between 10–2 (——) and 102 (–··–), except for L between 10–6 and 102. The other parameters are kept constant at unity except for a=100 and L=1/100. Varying parameter A m in panel F is the same as changing the modulator concentration [M]. As can be observed, changing [M] with the above fixed parameter values has no effect on the concentration–occupancy curve. Inserts in panel F, b>1 and c>1. A m or [M] is varied as before, but now with b=100 and a=1000 in the left insert. Compare with Fig. 3b in Hall (2000). In the right insert c=30 with a=1. Compare with Hall's Fig. 4b (2000).

7.2.3. Conclusion on Aspects of Binding as Simulated with the ATSM

The ATSM can model the effects of a modulator molecule on reduced as well as enhanced binding of an orthosteric ligand by changes in apparent affinity. The ATSM is well suited for a simulation of a modulator molecule that can either increase or decrease (a) binding of a radio-ligand, and it can also imitate (b) positive or negative effect on binding affinity of a neutral ant-agonist, and (c) positive or negative effects on binding affinity of an inverse agonist. All these effects can be observed in laboratory experiments with agonist, ant-agonist, and inverse agonist both as radio-ligands and non-labeled compounds. For binding, compare Fig. 7.7A with Hall's Fig. 3a, Fig. 7.7C with Hall's Fig. 4a, and Fig. 7.8A with Hall's Fig. 5a.

Figure 7.7. Variation of parameters b and c in ATSM for binding with the modulator [M] as independent variable. (A) System constant b is varied in five steps between 10–2 (——) and 102 (–··–) and the other parameters were a=1000, c=1, d=1, L=0.01, A ss=1, and A mm=1 with S fixed at =0.1. Compare with Hall's Fig. 3a (2000). (B) A 3-D presentations of ATSM with the same parameter values as in panel A for b varying in 3 steps: 10–2, 100, and 102. The 3-D presentations are cut off at S ([Agonist at ‘O’]) above 0.1. (C) System constant c is varied in five steps between 10–2 (——) and 102 (–··–) and the other parameters were, a=1, b=1, d=1, L=0.01, A ss=1, and A mm=1 with [S] fixed at =1. As indicated by Hall (2000), a = 1 may be taken as a neutral ant-agonist and c is independent of a. Compare with Hall's (2000) Fig. 4a. (D) A 3-D presentation of ATSM with the same parameter values as in panel C for c varying in 3 steps: 10–2, 100, and 102. The 3-D presentation is cut off at S ([agonist at ‘O’]) above 1.

Figure 7.8. Variation of parameter d in ATSM for binding with the modulator as independent variable. (A) System constant d is varied in five steps between 10–2 (——) and 102 (–··–) and the other parameters were a=10,000, b=1, c=1, L=0.01, A ss=1, and A mm=1 with [S] fixed at =0.01. Compare with Hall's (2000) Fig. 5a. (B) A 3-D presentation of ATSM with the same parameter values as in panel A for d varying in 3 steps: 10–2, 100, and 102. The 3-D presentation is cut of at S ([Agonist at ‘O’]) above 0.01. (C) The same as in panel A but with [S] = 1 instead of 0.01 and a=0.1 instead of 10,000. Compare with Hall's (2000) Fig. 5d. (D) A 3-D presentation of ATSM with the same parameter values as in panel C for d varying in 3 steps: 10–2, 100, and 102. The 3-D presentation is cut off at S ([Agonist at ‘O’]) above 1.

7.3. Aspects for Function with the ATSM

7.3.1 Distribution Equation

The conservation principle for function in the ATSM (Fig. 7.1A) may be formulated in a distribution equation given here as:

Symbols used are listed in Box 7.1 and Table 7.1. Compare Eq. 7.2 with Hall's equivalent expression (Hall 2000, Eq. 10). As for binding, in Hall's analysis it is the concentration of either an agonist or a modulator that is varied as independent variable, not concomitantly.

7.3.2 Variation in Single Parameters

Examples of variation in single model parameters for the functional ATSM (listed in Box 7.1) are shown in Fig. 7.9 for increasing concentrations of the agonist. See also the legend to Fig. 7.9 for further details.

Figure 7.9. Examples of parameter changes in ATSM for function with the agonist [S] as independent variable. (A) In panels A–F, parameters b, c, d, L, A ss, and A mm, as indicated, vary in five steps between 10–4 (——) and 104 (–··–), for L between 10–6 and 102. The other parameters are kept constant at unity except for a=100 and L=1/100. Varying parameter A m in panel F is the same as changing the modulator concentration [M]. As can be observed, changing [M] with the above fixed parameter values has no effect on the dose-response curve. Changing [M] with either b, c, d, or L different from unity is demonstrated in Figs. 7.10 11 12 13 .

7.3.3. Conclusion on Functional Aspects as Simulated with the ATSM

Observing the effects of parameter variation in the ATSM demonstrates how ATSM can simultaneously explain both positive and negative heterotropic modulation of function; by variation in b, c and d. For function, compare Fig. 7.10A + C with Hall's Fig. 6a,b, Fig. 7.11A + C with Hall's Fig. 7a,b, and Fig. 7.12A + C with Hall's Fig. 8a,c. In Fig. 7.13, there are plots of the functional ATSM for parameter L below and above 1/10 with modulator concentration [M] varying in five steps. Notice the spontaneous activity in plots of Fig. 7.13C + D for L at 0.2.

Figure 7.10. Variation of modulator concentration [M] in ATSM for function with the primary ligand S as independent variable. In A and B: b>1, in C and D: b<1. System constants were as follows: a=30, b=100 or 0.1, c=1, d=100, L=0.01, A ss=1, and A mm=1. [M] was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C by a factor 102 and in the 3-D plots arrows indicate [M] at 10–4 and 104. The arrows in the 3-D concentration plane illustrate the [M] used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D. Compare panels A and C with Hall's (2000) Fig. 6a,b.

Figure 7.11. Variation of modulator concentration [M] in ATSM for function with the primary ligand S as independent variable. In A and B: c>1, in C and D: c<1. System constants were as follows: a=300, b=10, c=100 or 0.01, d=1, L=0.01, A ss=1, and A mm=1. [M] was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C by a factor 102 and in the 3-D plots arrows indicate [M] at 10–4 and 104. The arrows in the 3-D concentration plane illustrate the [M] used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D. Compare panels A and C with Hall's (2000) Fig. 7a,b.

Figure 7.12. Variation of modulator concentration [M] in ATSM for function with the primary ligand S as independent variable. In A and B: d>1, in C and D: d<1. System constants were as follows: a=300, b=3, c=1, d=100 or 0.01, L=0.033, A ss=1, and A mm=1. [M] was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C by a factor 102 and in the 3-D plots arrows indicate [M] at 10–4 and 104. The arrows in the 3-D concentration plane illustrate the [M] used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D. Compare panels A and C with Hall's (2000) Fig. 8a,b.

Figure 7.13. Variation of modulator concentration [M] in ATSM for function with the primary ligand S as independent variable. In A and B system constants were as follows, L=0.01, a=1000, b=1, c=8.5, d=0.03, A ss=1, and A mm=1. In C and D, L=1/3 and the other system constants were as follows: a=50, b=0.5, c=1, d=0.01, A ss=1, and A mm=1. [M] was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C by a factor 102 and in the 3-D plots arrows indicate [M] at 10–4 and 104. The arrows in the 3-D concentration plane illustrate the [M] used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D. Compare panel A with Hall's (2000) Fig. 10a.

7.4 Further Aspects for the ATSM

Hall (2000) has analyzed the behavior of the ATSM for binding and function when either the agonist concentration or the modulator concentration varies as the independent variable. Some aspects on his analyses are presented and debated below.

7.4.1. Co-agonism

Co-agonism is seen for ligand-gated receptor channels (LGC) where a simultaneous presence of two different ligands is necessary to open the channel (Johnson & Ascher 1987; Dingledine et al. 1990; Marvizon & Baudry 1993; Corsi et al. 1996, Hall 2000). Other systems have been suggested with co-agonism, for instance, GABA-A channel receptors (Chang & Weiss 1999; Rusch et al. 2004; Rusch & Forman 2005).1

7.4.2. Discrepancy of Modulator Effects on Binding and Activity

In the ATSM with the right selection of parameters, as shown by Hall (2000), the occupancy with primary ligand S in a double-sited receptive unit can be almost independent of the concentration of an allosteric ligand, ranging over several orders of magnitude. Thus, the activity of the receptive unit may vary with changing concentrations of the allosteric ligand M, [M] (Fig. 7.13A and Hall's Fig 10a) with parameters a=10,000, b=1, c=8.5, d=0.03 and L=0.001, while the binding of S is almost constant for certain parameter values over an extended span of concentrations for [M] (Fig. 7.14A and Hall's Fig. 10b). This is an illustration of the situation for allosteric compound CPCCOEt at the human mGluR1b receptor, where there is no detectable change in glutamate binding for changing concentrations of CPCCOEt, while the activity changes with the change in [CPCCOEt] (Litschig et al. 1999; Hall 2000).

Figure 7.14. The fractional binding in ATSM as a function of modulator concentration [M]. (A) Parameter values are a=10,000, b=1, c=8.5, d=0.03, L=0.001, A ss=1, and A mm=1 with S varying between 10–5 and 103. Thus, for certain values of parameters, the occupancy appears independent of modulator concentration over a large range of agonist concentrations. Compare with Hall's Fig. 10b. Conversely, as demonstrated in Fig. 7.13A, the activity changes with variation in [M] as a function of [S]. See Hall's Fig. 10a. This discrepancy between occupancy and function of the ATSM explains observed effects with modulator CPCCOEt (Litschig et al. 1999). (B) However, observe that with [S] for instance fixed at =0.1, the concentration-occupancy relation in panel A is not absolutely straight as demonstrated by expanding the ordinate just above 50% in panel B. (C) To keep the concentration-occupancy constant for variations in [M] is a delicate balance of parameter selection as in panel A. It is a balance which easily breaks down with minor changes, here for instance in parameter d, with [S] still at 0.1 and the other parameters as in panel A. The ATSM is a ‘cliff-hanger’ model for effects of modulators as CPCCOEt.

However, for this to occur with the ATSM, demands are high on selected parameters for the ATSM. Parameter L must be ≤1/1000, a must be ≥10,000, and concomitantly only narrow value ranges for d and c will qualify. The span for the right selection of parameter values is minute, as already indicated in Fig. 7.14B. Additionally, observe for instance the effects of minute variations of parameter d in Fig. 7.14C. The ATSM parameter values have to be within extremely narrow ranges to follow the action of CPCCOEt, therefore the above explanation for the behavior of CPCCOEt does not appear the most likely one, although it is a possibility.

7.5. Aspects for Binding with the HOTSM

7.5.1 A Distribution Equation for Occupancy in HOTSM

The HOTSM is shown in Fig. 7.1B, and its parameters and receptor conformations are listed in Table 7.1 and Box 7.1. A distribution equation for receptor occupancy in HOTSM based on the conservation principle must include factorization of receptor conformations where two identical ligands are bound at the same time.

We may list two slightly different occupancy models for the homotropic two-state reaction scheme. One in which we double the value of conformations in the reaction scheme with due respect to the fact that in conformations with two ligands bound, SRS and SR*S, they count twice. This first model merely includes the possible conformations. Therefore, the binding version of HOTSM is formulated with ao as actual occupancy:

The second formulation of binding in HOTSM further takes into account that all conformations have a double possibility of binding ligands, therefore binding is related to a double set of conformations by simply doubling all conformations in the denominator of the distribution equation equal to all possible binding situations (see, e.g., Cornish-Bowden 2004, Eq. 11.30 or Cornish-Bowden 1995, Eq. 9.27).2 This distribution equation is thus given by:

Note that in Eqs. 7.3 and 7.4, some system constants now appear with an apostrophe in order to differentiate them from similar parameters in the ATSM. Further, here for HOTSM, the association constant A ms replaces A mm of the ATSM.

A general form of Eq. 7.4 is given by Cornish-Bowden (2004), Eq. 11.32, 1995, Eq. 9.29) and a variant form is derived by Kurganow (1982), Eq. 3.41). Compare this with the equations for the random reaction scheme for binding in Chapter 6–Eq. 6.4 for two ligands and Eq. 6.8 for n ligands.

7.5.2. Effects of Variation in Single Parameters of Binding-HOTSM

A preliminary analysis of varying single parameters in the occupancy version of the HOTSM has been carried out (Bindslev 2004). Meanwhile in binding, this analysis did not take into account that some of the receptor conformations have to count twice when two agonist molecules are bound simultaneously, and as a further possibility that all bound conformations should be related relative to all possible conformation multiplied by a factor 2, as in Eq. 7.4 (see the arguments for Eq. 7.4 above).

An analysis of these aspects of the HOTSM in binding is presented here. A demonstration of the effects of varying single parameters in occupancy-HOTSM by the two principles in Eqs. 7.3 and 7.4 are shown for parameters b’, c’, d’, L, A ss, and A ms in Figs. 7.15 16 17 18 19 20 and detailed comments may be found in the figure legends.

Figure 7.15. Variation of parameter b’ in HOTSM for binding with the agonist/modulator as independent variable. A and B are based on the conformation model for HOTSM, C and D on the probability model for binding (e.g., Cornish-Bowden 2004, Eq. 11.30). Parameter b’ was varied in five steps between 10–2 (——) and 106 (–··–) in panels A and C, while the other system constants were as follows: a=1, c’ = 1, d=1, L=0.1, A ss=1, and A ms=1. Compare with Bindslev's Fig. 10 (2004). In panel B and D the two surface plots are for b’=10–2 and 106. The arrows in the 3-D concentration plane illustrate the agonist/modulator concentration used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.16. Variation of parameter c’ in HOTSM for binding with the agonist/modulator as independent variable. A and B are based on the conformation model for HOTSM, C and D on the probability model for binding (e.g., Cornish-Bowden 2004, Eq. 11.30). Parameter c’ was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C, while the other system constants were as follows: a=1, b’ = 1, d’ = 1, L=0.1, A ss=1, and A ms=1. Compare with Bindslev's Fig. 10 (2004). In panel B and D, the two surface plots are for c’=10–4 and 104. The arrows in the 3-D concentration plane illustrate the agonist/modulator concentration used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.17. Variation of parameter d’ in HOTSM for binding with the agonist/modulator as independent variable. A and B are based on the conformation model for HOTSM, C and D on the probability model for binding (e.g., Cornish-Bowden 2004, Eq. 11.30). Parameter d’ was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C, while the other system constants were as follows: a=1, b’ = 1, c=1, L=0.1, A ss=1, and A ms=1. Compare with Bindslev's Fig. 10 (2004). In panel B and D, the two surface plots are for d’ = 10–4 and 104. The arrows in the 3-D concentration plane illustrate the agonist/modulator concentration used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.18. Variation of parameter L in HOTSM for binding with the agonist/modulator as independent variable. A and B are based on the conformation model for HOTSM, C and D on the probability model for binding (e.g., Cornish-Bowden 2004, Eq. 11.30). Parameter L was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C, while the other system constants were as follows: a=1, b’ = 1, c=1, d’ = 1, A ss=1, and A ms=1. Compare with Bindslev's Fig. 10 (2004). In panel B and D, two overlaid surface plots are for L=10–4 and 104. The arrows in the 3-D concentration plane illustrate the agonist/modulator concentration used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.19. Variation of parameter A ss in HOTSM for binding with the agonist/modulator as independent variable. A and B are based on the conformation model for HOTSM, C and D on the probability model for binding (e.g., Cornish-Bowden 2004, Eq. 11.30). Parameter A ss was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C, while the other system constants were as follows: a=1, b’ = 1, c=1, d’ = 1, L=0.1, and A ms=1. Compare with Bindslev's Fig. 10 (2004). In panel B and D the two surface plots are for A ss=10–4 and 104. The arrows in the 3-D concentration plane illustrate the agonist/modulator concentration used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.20. Variation of parameter A ms in HOTSM for binding with the agonist/modulator as independent variable. A and B are based on the conformation model for HOTSM, C and D on the probability model for binding (e.g., Cornish-Bowden 2004, Eq. 11.30). Parameter A ms was varied in five steps between 10–4 (——) and 104 (–··–) in panels A and C, while the other system constants were as follows: a=1, b’ = 1, c=1, d’ = 1, L=0.1, and A ss=1. Compare with Bindslev's Fig. 10 (2004). In panel B and D, the two surface plots are for A ms=10–4 and 104. The arrows in the 3-D concentration plane illustrate the agonist/modulator concentration used for the two extreme plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

With the extension in Eq. 7.4 compared to Eq. 7.3, the theory predicts a plateau level at 50% occupancy for certain parameter values, shown for b’ in Fig. 7.15, for c’ in Fig. 7.16, for A ss in Fig. 7.19, and for A ms in Fig. 7.20. However, in spite of the modification of Eq. 7.3 to Eq. 7.4, the introduced ‘flexibility’ of a plateau at 50% is not worth much in binding displacement experiments, which do not seem restrained to a plateau phase at exactly 50% occupancy; though, it may reasonably be argued to include this correction for the distribution equation of the binding-HOTSM.

7.5.3 Conclusion on Aspects of Binding as Simulated with the HOTSM

Positive co-operativity in occupancy can be simulated with HOTSM. A variant form of the HOTSM can describe terraced concentration binding curves, but merely terraced for binding around 50% occupancy. Also, compare these 50% plateaus with similar published plateaus for the square geometry in the KNF model by Koshland et al. (1966).

How is binding in HOTSM different from binding in ATSM? Binding in the HOTSM includes double occupancy with co-lateral binding, whereas in the ATSM co-lateral binding is excluded (compare plots in Figs. 7.15 7.16 7.17 7.18 7.19 7.20 with plots in Fig. 7.6).

The positive co-operativity of the HOTSM is insufficient for positive co-operativity in, for instance, hemoglobin that operates with four binding sites. Here, the HOTSM needs to be developed further.

7.6. Aspects for Function with the HOTSM

7.6.1 Simulation with HOTSM

The functional form of the homotropic two-state reaction model has several and surprising possible simulations available (Bindslev 2004). To illustrate, a research example, that is discussed in more detail in sub-chapter 7.7, is brought up here in the introduction (section 7.6.2), for which the HOTSM seems relevant in simulation of experimental self-inhibited and self-enhanced effects as reverse bell-shaped dose-responses.

Intrinsic self-enhancement of oxygen binding to hemoglobin following a deviation from the ordinarily hyperbolic load relation is not relevant for the functional form of HOTSM, while self-inhibition in the form of substrate inhibition, as described in Chapter 3, is relevant. The former is designated as positive co-operative, while the latter may be denoted as negative co-operative. However, since substrate inhibition can result in bell-shaped dose-responses (see Figs. 3.2 and 3.3), a designation such as ‘auto-ant-agonism’, ‘auto-inhibition’, or ‘negative auto-intervention’ for the decaying leg of bell-shaped dose-responses are more valid than ‘negative co-operativity’ because negative co-operativity is usually a reference to merely shallow dose-responses (Chapter 15) (Bindslev 2004). A general term when assuming genuine allostery which covers all the observed deviations from simple load relations, even reversed bell-shaped dose-responses, is auto-modulation. We may operate with positive and negative auto-modulation. HOTSM is well suited to simulate both positive and negative functional auto-modulation, including, as we shall see, reverse bell-shaped and reverse terraced dose-response relationships.

7.6.2 Examples of Reverse Bell-shaped Dose-responses

Due to significant spontaneous activity, typically in model systems of transfected cell, it is possible to observe reverse bell-shaped dose responses (Migeon & Nathanson 1994; Michal et al. 2001; Christopoulos et al. 2001; Accomazzo et al. 2002; Nasman et al. 2002; Hornigold et al. 2003; Holmqvist et al. 2005). These reverse bell-shaped dose-responses are not necessarily explained, as usual, by a G protein signal-bifurcation between for instance Gs and Gi complexes, since the bell-shaped dose-response may be present even after eliminating the function of Gi proteins by treatment with PTX of the receptor-transfected cells (Hornigold et al. 2003). Meanwhile, the observed reverse bell-shaped behavior is easily modeled by the HOTSM, suggesting an intrinsic receptor modulation as explanation. Thus, the HOTSM may also be used as an analytical tool for reverse bell-shaped dose-responses. This aspect of the HOTSM together with its inherent reversed terraced dose-responses is further discussed in sub-chapter 7.7. Models are discussed in Section 7.11.3 for the reverse bell-shaped dose-response of GPCRs where only a single type of G protein is present.

7.6.3. The Basics of HOTSM

The HOTSM is shown in Fig. 7.1B. It is a combination of the auto-intervention model described in Chapter 3 and the cTSM presented in Chapter 5. The receptor conformations and parameters for HOTSM are listed in Table 7.1 and Box 7.1, and its distribution equation for function can be formulated as:

see Bindslev (2004), Eq. 2). Again, note the parameters with apostrophe in order to differentiate them from similar parameters in the ATSM. The association constant A ms replaces the A mm of ATSM.

7.6.4. Variation in Single Parameters

In contrast to the ATSM, both concentrations of ligands at either binding site are varied simultaneously as they change concomitantly. This is compulsory for the HOTSM. I shall return to how this can be mimicked in the ATSM in sub-chapter 7.8.

For the functional form of the HOTSM, variations in parameters b’, c’, d’, L, A ss, and A ms are shown in Figs. 7.21 22 23 24 25 26 . Details of the parameter variations are commented on in the figure legends.

Figure 7.21. Variation of parameter b’ in HOTSM for function with increasing concentration of agonist/modulator as independent variable. For b’ > 1 in A and B, and for b’ < 1 in C and D. Parameter b’ was varied in five steps between 100 (——) and 104 (–··–) in panel A, and between 10–4 (——) and 100 (–··–) in panel C, while the other system constants were as follows: a=1000, c’ = 1, d=1, L=0.01, A ss=1, and A ms=1. Compare with Bindslev's Fig. 3 (2004). The three surface plots in panel B are for b’ = 100, 102, and 104 and in panel D for b’ = 10–4, 10–2, and 100. The arrows in the 3-D concentration plane in panels B and D illustrate the agonist/modulator concentration used for the plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.22. Variation of parameter c’ in HOTSM for function with increasing concentration of agonist/modulator as independent variable. For c’ > 1 in A and C, and for c’<1 in C and D. Parameter c’ was varied in five steps between 100 (——) and 104 (–··–) in panel A, and between 10–4 (——) and 100 (–··–) in panel C, while the other system constants were as follows: a=1000, b’ = 1, d=1, L=0.01, A ss=1, and A ms=1. Compare with Bindslev's Fig. 7 (2004). The three surface plots in panel B are for c’ = 100, 102, and 104, and in panel D for c’ = 10–4, 10–2, and 100. The arrows in the 3-D concentration plane in panels B and D illustrate the agonist/modulator concentration used for the plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.23. Variation of parameter d’ in HOTSM for function with increasing concentration of agonist/modulator as independent variable. For d’ > 1 in A and B, and for d’ < 1 in C and D. Parameter d was varied in five steps between 100 (——) and 104 (–··–) in panel A, and between 10–4 (——) and 100 (–··–) in panel C, while the other system constants were as follows: a=30, b’ = 3, c=0.1, L=1/9, A ss=1, and A ms=100. Compare with Bindslev's Fig. 4 (2004). The three surface plots in panel B are for d’ = 100, 102, and 104, and in panel D for d’ = 10–4, 10–2, and 100. The arrows in the 3-D concentration plane in panels B and D illustrate the agonist/modulator concentration used for the plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.24. Variation of parameter L in HOTSM for function with increasing concentration of agonist/modulator as independent variable. Parameter L varied in five steps from 0.2×10–2 (——) to 0.2×102 (–··–) by a factor 10 in panels A and C and in three steps from 0.2×10–2 to 0.2×102 in panels B and D by a factor 102. Other parameter values were as follows in A and B: a=1000, b’ = 1, c=1, d’ = 1, A ss=1, and A ms=1. In panels C and D for a=1000, b’ = 1, c=0.01, d’ = 0.01, A ss=1, and A ms=1. See for instance Fig. 2 in Bindslev (2004). The arrows in the 3-D concentration plane in panels B and D illustrate the agonist/modulator concentration used for the plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.25. Variation of parameter A ss in HOTSM for function with increasing concentration of agonist/modulator as independent variable. A ss was varied in five steps between 100 (——) and 104 (–··–) in panel A, and between 10–4 (——) and 100 (–··–) in panel C, while the other system constants were as follows: a=1000, b’ = 1, c’ = 1, d=1, L=0.01, and A ms=1. Compare with Bindslev's Fig. 5 (2004). The three surface plots in panel B are for A ss=100, 102, and 104 and in panel D for A ss=10–4, 10–2, and 100. The arrows in the 3-D concentration plane in panels B and D illustrate the agonist/modulator concentration used for the plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

Figure 7.26. Variation of parameter A ms in HOTSM for function with increasing concentration of agonist/modulator as independent variable. A ms was varied in five steps between 100 (——) and 104 (–··–) in panel A, and between 10–4 (——) and 100 (–··–) in panel C, while the other system constants were as follows: a=1000, b’ = 1, c’ = 1, d=1, L = 0.01, and A ss=1. Compare with Bindslev's Fig. 6 (2004). The three surface plots in panel B are for A ms=100, 102, and 104, and in panel D for A ms=10–4, 10–2, and 100. The arrows in the 3-D concentration plane in panels B and D illustrate the agonist/modulator concentration used for the plots in panels A and C, thus, the cut off in topography in 3-D by a plane raised at the arrows will produce the indicated plots in 2-D.

7.6.5 Conclusion on Functional Aspects as Simulated with the HOTSM

The functional HOTSM is a versatile tool for simulation of single-ligand dose-responses rendering load-deviant curves. Experimental curve forms that may be simulated with HOTSM are summarized in Fig. 7.27. In the past, such load-deviant curve forms of experimental synagics triggered use of the Hill and sums of the Hill equation as analytical tools (e.g., Bronnikov et al. 1999; Accomazzo et al. 2002; Hornigold et al. 2003). The HOTSM offers a mechanistic approach to the analysis, while the Hillian approach is nearly always a pure mathematical description lacking any mechanistic relevance (Chapter 10). Hence, future modeling ought to consider variations on the themes of the HOTSM, replacing versions of the Hill equation.Fig. 7.28

Figure 7.27. The response level of the functional-HOTSM has three plateaus. A first response-plateau in the absence of ligands is equal to 1/[1 + 1/(L)]. A third response-plateau at high agonist concentrations is given by 1/[1 + 1/(L·a·b’·d’)]. Finally, a second-response plateau at intermediate ligand concentrations is dependent on all seven system constants. The three plateaus may be above, at, or below each other, panels A–D, rendering dose-response curves that are positive or negative co-operative, bell-shaped or terraced, as well as reverse bell-shaped or reverse terraced. Parameter values in panels A–D were arbitrarily chosen in order to demonstrate several aspects of the model. A ss varied in three steps from 10–2 (----) to 102 (-·-·-·) by a factor 100 in panels A and B, while A ms varied in three steps from 0.3 (----) to 3000 (-·-·-·) by a factor 100 in panels C and D. The remaining parameters were fixed. Thus in panel A: L=1/4, A ms=20, a=100, b’ = 0.01, c’ = 0.01, and d’ = 17. In panel B: L=4, A ms=20, a=10, b’ = 0.01, c’ = 0.01, and d’ = 0.7. In panel C: L=1/4, A ss=1, a=100, b’ = 0.01, c’ = 0.01, and d’ = 17. In panel D: L=4, A ss=1, a=10, b’ = 0.01, c’ = 0.01, and d’ = 0.7. (After Bindslev (2004), Fig. 2).

Figure 7.28. An additional example of HOTSM in function, where b’ vary. Parameters were a=30, c’ = 0.1, d’ = 1, A ss=1, A ms=100, and L=1/9. Parameters b’ was varied as indicated in the panels. Unbroken line-curves represent the lowest b’ value.

7.7. Further Aspects for the HOTSM

With the right combination of parameter values, the functional HOTSM possesses both bell-shaped and terraced as well as reverse bell-shaped and reverse terraced relationships (Fig. 7.27). A prerequisite for observing inverse agonism, reverse terraced curves, and reverse bell-shaped dose-response relations is spontaneous activity which is detectable for L>1/100. Parameter L is the sole system constant to determine the initial level of activity in both HOTSM (Fig. 7.23 and Fig. 7.24) and ATSM (Fig. 7.13C + D).

7.7.1. Reverse Bell-shaped Synagics

Hornigold et al. (2003) described a system only expressing a single receptor subtype, m3, coupled to producing cAMP as spontaneous activity and affected in a reverse bell-shaped dose-response relation by methacholine (Fig. 7.29A). Even after treatment with PTX in order to eliminate possible coupling to inhibitory G proteins, this system still displays reverse bell-shaped synagics (Fig. 7.29B). This excludes bifurcation via Gs and Gi as an explanation for the observed reverse bell-shaped dose-response relationship.

Figure 7.29. A fit of HOTSM to experimental data. (A) Dose-response data of cAMP production in Chinese hamster ovary (CHO) cells expressing muscarinic m3 receptor subtypes and a HOTSM curve fitted by adjusting parameters. Fitted parameters were a=0.401, b’ = 0.534, c’ = 0.145, d’ = 11.1, A ss=0.142, and A ms=27.3 with L sat at 1/3. Data taken from Hornigold et al. (2003). (B) Control and petussis toxin (PTX), effects on cAMP production in CHO cells expressing only m3 receptor subtypes. If any effect, PTX augments the concave bell-shaped dose-response relation. Copy of Fig. 5a in Hornigold et al. (2003).

The concave bell-shaped dose-response in Fig. 7.29A can be simulated by the HOTSM (see curve in the figure) and may thus be interpreted as a receptor intrinsic modulation in line with the evidence of more than one binding site for ligands in functional muscarinic receptors (Huang & Ellis 2007).

Therefore, for a possible alternative explanation of the observed reverse bell-shaped dose-response curve, as mentioned above, the HOTSM is a simple, elegant, and sufficient analytical tool, as also pointed out by Bindslev (2004), Figs. 2 and 9).

7.7.2. The HOTSM Against the Hill Equation

The change between three plateau levels of the HOTSM in function (Fig. 7.27), and also the co-operative rise of the concentration-occupancy relation for the model in binding ( Figs. 7.15 7.16 7.17 7.18 7.19 7.20 ), are not as steep as dose-responses described by a single or a combination of two or more Hill equations. This is not surprising, as the whole analysis of the HOTSM is based on a system with only two binding sites. For a fair comparison, the HOTSM should be allowed to have as many sites as the Hill analyses in such situations (Bindslev 2004).

7.7.3. Pseudo-homotropic Two-state Model (PHOTSM)

*

A model similar to the extended two-state model for heterotropic interaction (ETCM) in Section 7.10.1, but now for homotropic behavior, co-operativity, has been put forward recently as two-state two-site model or a ‘two-state-dimer-model’ (Franco et al. 2005; Franco et al. 2007a b; Ehlert 2008). The former model, which is a fraction of the HOTSM, is capable of simulating reverse bell-shaped dose-response relations without an over-shoot (compare Fig. 3 in Franco et al. (2005) and Fig. 7.27). Obviously, the ‘two-state-dimer-model’ and the two-state two-site model are covered by the HOTSM. See Table 7.3 for more models.

7.8 Comparison between the ATSM and the HOTSM

From a functional viewpoint, the ATSM overlap with the HOTSM. From the aspect of occupancy, the two models cover separate issues and are very different.

7.8.1. Functional Studies

The HOTSM can be said to be a special version of the ATSM. The ATSM is congruent with the HOTSM when doses are applied from a mixture made of an orthosteric ligand and a modulatory ligand ( Figs. 7.21 7.22 7.23 7.24 7.25 7.26 ). As the concentration of ligands in the mixture is increased, a fixed proportion is kept between doses of the two ligands in the ATSM, just as for the HOTSM regime. Such an experiment with a primary ligand and a modulator in a fixed ratio will, in theory, follow the HOTSM when corrections are made for possible differences between affinity constants A mm and A ms, which in the model may simply be taken into account by multiplying either by a factor. Now the fixed ratio experiment has become a single-ligand study, and its analysis by the ATSM is identical to an analysis by the HOTSM. Aspects of the behavior of ATSM/HOTSM as system tools were presented by Bindslev (2004).

7.8.2. Binding Studies

The ATSM is a prime example for the simulation and analysis of dose-responses of occupancy in radio-ligand binding studies as affected by cold allosteric ligands, demonstrated by Hall (2000). In the HOTSM, application of a ligand that can bind to two sites simultaneously is different from the drug application described by Hall (compare Figs. 7.6 7.7 7.8 with Figs. 7.15 7.16 7.17 7.18 7.19 7.20 ).

To recapitulate, there is no feasible comparison between the ATSM and the HOTSM reaction schemes used in simulations for binding. The reaction scheme of ATSM can only bind ligands of interest for occupancy studies in a single site. Co-lateral binding is excluded. In the HOTSM, two ligands of interest can bind simultaneously in a co-lateral fashion by just adding one type of ligand. This clearly separates the two models from each other in binding studies.

7.9. Evaluation of FP-TSM as an Analytical Tool

7.9.1. Use of the FP-TSM Tool is Easy?

The FP-TSMs in Fig. 7.2 and 5.9 has 18 conformations and 17 independent system constants, and besides its genuine two-state isomerization, it involves two different ligands, an agonist or a substrate and a modulator, with mutual co-lateral binding. The possibility of co-lateral binding in FP-TSM differentiates it from the CQCM in Fig. 7.5.

Armed with some reasonable assumptions, the 17 constants can be reduced to around 12.

Still, to analyze the FP-TSM with 12 independent system constants sounds horrendous, but I believe it is only a matter of time before this model with unliganded spontaneous activity will be employed in the analysis of ligand–receptor interactions. How can I trust such a development? Well, consider all the efforts being put into drug design and drug development based and performed on the modeled dynamics of steric coordinates in ligand binding sites, and hold this against the immense degrees of freedom in possible alignment for receptive systems and their cognate ligands at an Angstrom scale. The possibilities are legio and will keep scientists and support floating for years – not necessarily with any foretold causal solutions.

7.9.2. Structure-activity-relationship (SAR)

The whole G protein-coupled receptor community eagerly awaits more details of the structural resolution of GPCRs, both for ligand binding sites and for the coupling between GPCRs and G proteins (Oldham et al. 2006; Kobilka & Scherkler 2008), as well as the interaction between GPCRs – G proteins and other proteins (Ferguson 2007), and looks enviously to the advancements of nearly 40 years in structural resolution of enzyme function. As the structural secrets are being revealed for transporters as well (Toyoshima 2007), porter-people are alert and expectations are high.

Meanwhile, even with recent promising breakthroughs in the structural resolution for receptors (Palechewski et al. 2000; Okada et al. 2002 2004; Palechewski 2006; Cherezov et al. 2007) and for several different transporters (Gouaux & MacKinnon 2005) as channels (Dutzler et al. 2002 2003; MacKinnon et al. 2003 2004; Unwin 2005; Jentsch et al. 2005a b; Miller 2006; Ramjeesingh et al. 2006; Jasti et al. 2007), pumps (Toyoshima et al. 2000; Obara et al. 2005; Jensen et al. 2006), exchangers (Dutzler et al. 2002 2003; Hunte et al. 2005; Jentsch et al. 2005a b; Accardi et al. 2006; Miller 2006; Nicoll et al. 2006; Nguitragool et al. 2006), ABC-transporters (Yu et al. 2003; Dawson & Locher 2006), and co-transporters (Yernool et al. 2004; Yamashita 2005; Indarte et al. 2007) leading the way to synthesis of new lead-compounds, dynamic possibilities in the binding pockets are still immense and overwhelming when it comes to the dynamic stereo-chemistry at a subatomic scale needed for the ultimate useful structure-activity-relationship (SAR).

7.9.3. The FP-TSM in SAR Analysis

In comparison with the efforts described in the preceding sections, to me it seems a modest task to develop the FP-TSM, which among other issues mechanistically may describe the reversal of auto-ant-agonism (reversal of negative auto-modulation or negative auto-intervention) by a classical competitive ant-agonist. An example of such reversal by a competitive ant-agonist has been observed for atropine at muscarinic receptors (Figure 30A) (Winding & Bindslev 1993), where atropine can reverse auto-inhibition induced by acetylcholine within one to two minutes, thus reactivating an auto-inhibited acetylcholine response of a cell membrane displayed muscarinic receptor pool (Fig. 7.30B).

In developing a scheme for this kind of ‘ant-agonism’, the terminology of designators in FP-TSM such as ‘competitive’ inhibitor, ‘non-competitive’ inhibitor, ‘interventor’, and ‘modulator’ may start to mingle. Atropine is a competitive ant-agonist in the classical sense (Brimblecombe 1974, Chapter 1; Carrijo et al. 1977; Hulme et al. 1978), but atropine will also appear as a non-competitive inhibitor or an interventor in intervention models, or a modulatory ligand in the ATSM and the FP-TSM. Although atropine is a ‘neutral’ ant-agonist, as it does not in itself evoke or prevent any activity in the absence of agonists or inverse agonists, it may have the characteristics of all four types of ligands mentioned above. However, atropine does not qualify as an ago-modulator, since it does not act on its own.

In the realm of G protein function, a start has been made with the implementation of a pre-FP-TSM (Downing et al. 2006). Furthermore, an alternative FP-TSM reaction scheme has been invoked for GABA-A channels (Rusch et al. 2004). The future looks bright for the FP-TSM as an analytical tool.

7.10. The Optimal-fixed-concentration-ratio (OFCOR)

7.10.1. Inhibition of Negative Auto-modification

In systems with negative auto-intervention or negative auto-modulation in the form of bell-shaped dose-responses at high agonist concentrations, this inhibition by the agonist itself may be counteracted by introducing a modulator or even a competitive ant-agonist, both with access to more than one agonist binding site.

As described above, a reactivation of such a negative auto-modulatory response was found for the ligand-pair acetylcholine/atropine in tracheal secretion (Fig. 7.30) (Winding & Bindslev 1993). There are several possible molecular explanations for this reactivation mechanism by atropine, one being competition at a secondary binding site with access for both inhibitory primary ligand (auto-inhibitory agonist acetylcholine) and a modulator (atropine). The FP-TSM is suited for such an explanation (see below).

Figure 7.30. Auto-inhibition by acetylcholine of tracheal secretion reactivated by atropine. (A) Low dose activation and high dose auto-inhibition of secretion (Isc) in hen tracheal epithelium by acetylcholine (ACh)?, followed by a reactivation-inhibition by atropine?. (B) Inhibition of secretion (Isc) by acetylcholine at high doses, total 512 M, and fast reactivation, t½ ≈60 s, by a single dose of atropine?, 100 M. ‘?’ = bumetanide at 10 M. Modified from Winding & Bindslev (1993, Figs. 1 and 3) with permission.

Figure 7.31. Examples of dose-response relation in the cubic ternary-complex model for G proteins (CTCM). Parameter d was varied in five steps. Values for parameters were as follows: a=30, b=3, c=0.1, Parameter d was varied in five steps between 10−5 (——) and 103 (–··–) in both panels A and B, while the other system constants were as follows: a=30, b=3, c’ = 0.1, L=1/9, A ss=1, and A ms=100. The concentration of the G protein as modulator, [G protein] at different expression levels, was fixed at 1 (arbitrary units) in A and at 100 (arbitrary units) in B as indicated by an arrow in the 3-D panels. The related contour plots are for d=10–1 in both panels with arrows indicating the G protein concentration. Compare the 3-D plots here with a 3-D plot in Fig. 7.23.

Figure 7.32. All kinds of models and mock-ups are possible, but not necessarily correct. (A) Parameter values for the ATSM a=30, b=3, c=vary in five steps by a factor 100 between 10–9 and 10–1, L=1/9, A ss=1, and A ms=100. (B) Can the allosteric two-state model describe Sidney's Opera house with the right architect?

Internalization of receptors as an explanation for the auto-inhibition in this system seems to be ruled out by the speedy reactivation within a minute or two by atropine (Fig. 7.30). Desensitization, converted to an attenuating dose-response, is another obvious possibility for the auto-inhibition by acetylcholine. Meanwhile, in this scenario, we have to explain how atropine behaves as a ‘competitive ant-agonist’ and how it resensitizes the auto-inhibition by acetylcholine in a time dependent manner.

As mentioned, the reactivation behavior can be addressed by a FP-OSM, which includes co-lateral binding at either site (Winding & Bindslev 1993). Ligand-relief of negative auto-modulation may of course also be simulated with the FP-TSM, which will be useful in case systems moreover comprise spontaneous activity.

7.10.2. Finding OFCORs

In case agonist concentrations in clinical treatment elicit auto-inhibition, thereby tempting a higher dose application with possible adverse effects, such auto-inhibition may be prevented by interventors, modulators, or ‘competitive ant-agonists’ with a ‘positive’ effect elicited from a secondary site. It will be possible to avoid adverse effects due to high concentrations of primary ligands that also elicit auto-inhibition by use of a single mixture of 'agonist' and 'auto-modulator' with an optimal concentration ratio of the two drugs. Of course, using such mixtures in dose-response experiments with optimum is the same as using OFCORs. The challenge is, of course, to find the right OFCORs for a given individual; his of hers drug-mixture.

Determining an optimal concentration ratio for experiments and treatments may be obtained by assaying the dose-response relations for and between the agonist and the modulator in several combinations of the two drugs. As an example, the OFCOR for a mixture was determined and verified in the study by Winding and Bindslev (1993) using a version of the FP-OSM. To date, a similar analysis using the FP-TSM is not in the literature.

7.11. More on Cubic Models and GPCR-related Models

7.11.1. The Extended Ternary Complex Model (ETCM)

In order to describe the dose-response relations obtained for GPCRs with constitutive activity and varied G protein levels, Samama and co-workers developed a model combining the original ternary-complex model equal to an intervention model (Chapter 2) for G proteins, with the cTSM (Chapter 5), in the simplest possible scheme (Samama et al. 1993; Lefkowitz et al. 1993). The ETCM for GPCRs is shown in Fig. 7.4B. Development of the ETCM was a stage on the road to the thermodynamically full model with modulators, including G proteins, and two-states combined in the CTCM (Fig. 7.3, Table 7.3) (Weiss et al. 1996a).

7.11.2. ATSM and CTCM and CQCM

When occupancy is considered, the ATSM (Hall 2000) for allosteric modulators is in principle identical to the CTCM for ‘modulation’ of GPCRs by G proteins (Fig. 7.1A and 7.4A) (Weiss et al. 1996a b). Examples of occupancy that cover both models are shown in Fig. 7.6.

The parallel occupancy between ATSM and CTCM is true also in simulation of experiments where the number of total receptors differs between radio-ligand agonist and radio-ligand ant-agonist data due to G protein-induced affinity changes (Kenakin 1997; Baker & Hill 2007; Ferguson 2007). Here, adjustments of parameter A mm come in to play.

On the other hand, in functional studies the two tools, ATSM and CTCM, behave somewhat differently (Weiss et al. 1996a c; Hall 2000). An illustration of this can be acknowledged by comparing the surface plot in 3-D panels of Figs. 23 and 31, where both concentrations of an agonist and the G protein as modulators are varied, as independent variables yielding the 3-D surface plots. Following the arrows in the 3-D concentration planes in Fig. 7.31A + B at an arbitrary concentration (expression level) of G protein exemplified with ‘1’ in Fig. 7.31A or at ‘100’ in Fig. 7.31B results in 2-D plots in A and B of Fig. 7.31, and should not be compared with the 2-D plots in Fig. 7.23 which are for the HOTSM. However, with simple modifications, the ATSM may as well be a model for GPCRs.

Experiments with varying expressions (concentrations) of G proteins are less common but have been performed, modeled, and discussed (Chidiac et al. 1996; Kenakin 1997; Tucek et al. 2001 2002; Kukkonen et al. 2001; Kukkonen 2004b).

With the CTCM it is possible to study GPCR dose-response relations with varying concentrations of the agonist as the independent variable and G proteins at different levels. The CTCM for G protein coupled to activated receptors, even in the absence of agonists GR*, obviously has possibilities of simulating spontaneous activity and inverse agonism, but in relation to non-G protein regulators, the CTCM is clearly unsatisfactory. Therefore, models such as the CTCM with an added third ligand are in demand and have been suggested as the CQCM (Christopoulos et al. 1998; Christopoulos & Kenakin 2002; May et al. 2004), but, to date, are not fully developed.

An amputated version of the CQCM with two G proteins, a GPCR, and an agonist yielding quaternary complexes3 but without an isomerization between two states, was implemented by Kukkonen et al. (2001) and Kukkonen (2004a), introducing a factor F controlling the level of receptors chelated by G proteins. Other simpler models, also without isomerization, including variations in G protein expression and a factor controlling maximal inhibition (MI) and maximal stimulation (MS) have been suggested (Tucek et al. 2002).

Needless to say, the CQCM illustrated in Fig. 7.5 should be developed sooner rather than later.

7.11.3. Models of ‘Ligand’ Numbers Matching ‘Receptor’ Numbers

In the latest development of models addressing the reverse bell-shaped dose-response relations for GPCRs, ‘dual effects’ (Jones et al. 1991; Kashihara et al. 1992; Migeon & Nathanson 1994; Vogel et al. 1995; Jakubik et al. 1996), two parallel approaches have been advocated for a limited G protein activation (Tucek et al. 2001 2002; Nasman et al. 2001; Kukkonen et al. 2001; Kukkonen 2004a). Tucek et al. (2002) and especially Kukkonen (2004a) give elegant descriptions of these models’ development and simulation capabilities; notwithstanding, arithmetical treatment in these two sources mutually overlook each others similarities in analysis.

Following Kukkonen, three principal models are invoked. The shuttling model, earlier described as a mobile receptor model (Jacobs & Cuatrecasas 1976) or collision coupling model (Tolkovsky & Levitzki 1978 1981), the complexing model (Chidiac 1998), and the precoupled model (Neubig 1994; Kukkonen 2004a). The latter comes in two versions; a precoupled-dependent and a precoupled-independent model. Somewhat different models on precoupling were also presented by Jacobs and Cuatrecassas (1976) and by Tolkovsky and Levitzki (1978).

In the recent models, the variation in a limited number of single type G protein that can be activated, is accomplished by a coefficient resembling the intrinsic activity factor of Ariens (1954) both in the Czech and the Swedish approach (Tucek et al. 2002; Kukkonen 2004a b). Thus, the limited number of active G proteins and the ratio between Gs and Gi proteins result in simulations that can match the observed reverse bell-shaped dose-responses. The simulated ‘spontaneous’ activity at 100% was obtained experimentally by forskolin stimulation.

7.11.4. Models with More than One-state for GPCRs

Model approaches to cover multi-states in GPCRs were recently presented, Table 7.3. For instance, one by Leff and co-workers (Leff et al. 1997; Scaramellini & Leff 1998 2002), another by Pin and co-workers (Parmentier et al. 2002) covering both two-states of a binding domain (see also Rovira et al. 2008) and two-state of an effector domain, and a third by Giraldo (2004). The models by Leff et al. and Giraldo both play on the theme of combining two cTSM for the sake of simulating GPCR functional behavior. These models operate with three or more states for the unliganded receptive unit, and may be elaborated further into multi-state models – a project left for the reader with some help from a recently published four-state model (Ehlert &Griffin 2008). A double-cubic model encompassing two states for both the receptor and its effector was already suggested in the 1970s (de Haen 1976).

7.11.5. Models for Binding with Dimerization as Isomerization

Dimerization of receptive units may lead to an allosteric interaction between binding sites in the dimer complex, thus resembling the switch between a reactive and an active state for a double-sited receptor. Recently, Durroux (2005) has suggested different and attractive models for such dimerization with ‘allosteric’ interaction, including cubic reaction schemes and four-pane two-state-like models. It will be interesting to follow development of Durroux’ models as they appear most relevant for a description and possible simulation of activity in dimerized GPCRs.

7.12. Use of Modulators in Therapy

7.12.1. Benefits of Modulators Compared with Orthosteric Ligands

An early example of modulators used in therapy with great success and few side-effects is the treatment of anxiety, depression, epilepsy, and sleeplessness with benzodiazepines, which started more than forty years ago. Nevertheless, the clinical aspect of the use of modulators or allosteric compounds is a relatively new issue, and has been advocated since the mid 1990s (Proska & Tucek 1994; Tucek & Proska 1995, Birdsall et al. 1995; Lazareno & Birdsall 1995). Some of the advantages of using modulator compounds instead of the natural endogenous agonist or artificial agonists with receptor subtype specificity can be listed as follows (Birdsall & Lazareno 2005):

(a) As the function of a modulator to a great extent is dependent on the presence and action of the physiological signal molecule, it is likely to mimic normal biological effects. The modulator is playing on the activity of an endogenous primary ligand.

(b) The enhancement of natural agonist activation or inhibition by modulators is less likely to reach adverse effects as for agonist drugs, even at high concentrations of the modulator compound, and therefore prolonged stimulation with high modulator concentrations may be induced without adverse effects as often seen with high concentrations of developed agonist drugs. There is ‘ceiling’ and ‘tune down’.

Besides reuptake and decay inhibitors, some of the expected therapeutic advantages with use of modulatory drugs may also be obtained with partial and partial inverse agonists (Fisas et al. 2006).

In clinical work with modulators, it is the functional aspect that is in focus. The picture of modulatory drug therapy has developed rapidly within the last few years for several receptor sub-families, and many developed drugs are in various phases of evaluation for introduction to the clinics (Table 7.2) (Gao & Jacobson 2006). Examples of biotech firms that develop such drugs are the NPS Pharmaceutical Company with a modulator for Ca sensors, and Xytis with a glycine site enhancer for NMDA channels under development (see homepages of NPS Pharmaceuticals Inc. and Xytis Inc.).

Of course, examples may appear where modulators are less beneficial than, for instance, weak partial agonists (Metha & Ticku 1999).

(c) Highly selective modulators for specific receptor subtypes are more likely to be found than subtype specific agonists, primary ligands, due to the evolutionary conservation of the primary binding sites for endogenous ligands (Birdsall et al. 1999; Spalding et al. 2002; Lanzafame et al. 2006; Langmead et al. 2006; Langmead & Christopoulos 2006). This is especially true for the muscarinic bindings sites. Actually, development of allosteric modulator drugs is an exploding field involving several effector molecules as indicated by recent publications listed in Table 7.2.

The receptor subtype selectivity by modulators, when judged to be based mainly on the conformational changes by induced fits or stabilization rather than on (high) affinity, is a kind of conformational change dubbed ‘absolute subtype selectivity’ (Lazareno et al. 2004).

During evolution, modulation may have been adapted for several secondary ligands with a more loose structure of their secondary binding sites. It has been difficult to find selectivity much above a factor 10 between primary ligands developed for muscarinic receptor subtypes (Alexander et al. 2006), whereas it seems likely that much higher selectivity may be found for modulators at the various muscarinic receptor subtypes (Birdsall & Lazareno 2005; Langmead et al. 2006). Other examples of allosteric modulators are for adenosine (Jacobson & Gao 2006; Childers et al. 2006), metabotropic glutamate (Kew 2004; Ritzen 2005; Shipe et al. 2005; Vauquelin & Liefde 2005; Vieira et al. 2005; Pin et al. 2005; Foster & Kemp 2006), and dopamine (Schetz 2005) receptor subtypes. Examples for ligand-gated ion channels (LGICs) are reviewed by Hogg et al. (2005) with benzodiazepines as the earliest prototype of modulators with allosteric responses at the GABA-A receptor (Olsen et al. 2004). For transporters, see, e.g., Maki & Dey (2006), Ghosh et al. (2006) and Borst et al. (2006), and other examples mentioned in sub-chapter 7.9 and Table 7.2.

Protocols and screens for new allosteric modulators for GPCRs are evolving (Langmead 2007; Raddatz et al. 2007; May et al. 2007; Avlani et al. 2007). Moreover, analytical tools are created to solve complex behavior in GPCR signaling (Fig. 7.31; Ehlert & Griffin 2008).

7.12.2. Summary of Modulator Benefits

A summary of the advantages of developed modulators compared with drug agonists is: (1) they mimic the function of endogenous primary ligands, (2) they are less likely to elicit adverse effects, and (3) they have possible greater receptor subtype selectivity than agonists and other molecules operating at a primary binding sites with its evolutionary conserved geometry.

Due to this rapidly expanding and potentially huge market for allosteric drugs, in the future, on theoretical grounds, it will be necessary to firmly implicate the ATSM and HOTSM for analysis of (allosteric) modulators in basic as well as applied approaches.

References

Accardi,  A, Walden,  M, Nguitragool,  W, Jayaram,  H, Williams,  C & Miller,  C.  Separate ion pathways in a Cl–/H+ exchanger.  J Gen Physiol  126:  563–570, 2005.  [Crossref]

Accardi,  A, Lobet,  S, Williams,  C, Miller,  C & Dutzler,  R.  Synergism between halide binding and proton transport in a CLC-type exchanger.  J Mol Biol  362:  691–699, 2006.  [Crossref]

Accomazzo,  MR, Cattaneo,  S, Nicosia,  S & Rovati,  GE.  Bell-shaped curves for prostaglandin-induced modulation of adenylate cyclase: two mutually opposing effects.  Eur J Pharmacol  454:  107–114, 2002.  [Crossref]

Alexander,  SPH, Mathie,  A & Peters,  JA.  Guide to receptors and channels.  Br J Pharmacol  147:  S8

2006.  [Crossref]

Ariëns,  EJ & de Groot,  WM.  Affinity and intrinsic-activity in the theory of competitive inhibition. Part III. Homologous decamehonium-derivatives and succinyl-choline-esters.  Arch Int Pharmacodyn  99:  193–205, 1954.

Avlani,  VA, Gregory,  KJ, Morton,  CJ, Parker,  MW, Sexton,  PM & Christopoulos,  A.  Critical role for the second extracellular loop in the binding of both orthosteric and allosteric G protein-coupled receptor ligands.  J Biol Chem  282:  25677–25686, 2007.  [Crossref]

Baggio,  LL & Drucker,  DJ.  Biology of incretins: GLP-1 and GIP.  Gastroenterology  132:  2131–2157, 2007.  [Crossref]

Baker,  JG & Hill,  SJ.  A comparison of the antagonist affinities for the Gi- and Gs-coupled states of the human adenosine A1-receptor.  J Pharmacol Exp Ther  320:  218–228, 2007.  [Crossref]

Baraldi,  PG, Iaconinoto,  MA, Moorman,  AR, Carrion,  MD, Cara,  CL, Preti,  D, Lopez,  OC, Fruttarolo,  F, Tabrizi,  MA & Romagnoli,  R.  Allosteric enhancers for A1 adenosine receptor.  Mini Rev Med Chem  7:  559–569, 2007.  [Crossref]

Benneyworth,  MA, Xiang,  Z, Smith,  RL, Garcia,  EE, Conn,  PJ & Sanders-Bush,  E.  A selective positive allosteric modulator of metabotropic glutamate receptor subtype 2 blocks a hallucinogenic drug model of psychosis.  Mol Pharmacol  72:  477–484, 2007.  [Crossref]

Bindslev,  N.  A homotropic two-state model and auto-antagonism.  BMC Pharmacol  4:  11

2004.  [Crossref]

Birdsall,  NJ & Lazareno,  S.  Allosterism at muscarinic receptors: ligands and mechanisms.  Mini Rev Med Chem  5:  523–543, 2005.  [Crossref]

Birdsall,  NJ, Cohen,  F, Lazareno,  S & Matsui,  H.  Allosteric regulation of G-protein-linked receptors.  Biochem Soc Trans  23:  108–111, 1995.

Birdsall,  NJ, Farries,  T, Gharagozloo,  P, Kobayashi,  S, Lazareno,  S & Sugimoto,  M.  Subtype-selective positive cooperative interactions between brucine analogs and acetylcholine at muscarinic receptors: functional studies.  Mol Pharmacol  55:  778–786, 1999.

Borst,  P, Zelcer,  N & van de Wetering,  K.  MRP2 and 3 in health and disease.  Cancer Lett  234:  51–61, 2006.  [Crossref]

Brauner-Osborne,  H, Wellendorph,  P & Jensen,  AA.  Structure, pharmacology and therapeutic prospects of family C G-protein coupled receptors.  Curr Drug Targets  8:  169–184, 2007.  [Crossref]

Brimblecombe,  RW    Drug Action on Cholinergic Systems .  London:  The Macmillan Press Ltd, 1974.

Bronnikov,  GE, Zhang,  SJ, Cannon,  B & Nedergaard,  J.  A dual component analysis explains the distinctive kinetics of cAMP accumulation in brown adipocytes.  J Biol Chem  274:  37770–37780, 1999.  [Crossref]

Bruheim,  S, Krobert,  KA, Andressen,  KW & Levy,  FO.  Unaltered agonist potency upon inducible 5-HT7(a) but not 5-HT4(b) receptor expression indicates agonist-independent association of 5-HT7(a) receptor and Gs.  Receptors Channels  9:  107–116, 2003.  [Crossref]

Buck,  M, Xu,  W & Rosen,  MK.  A two-state allosteric model for autoinhibition rationalizes WASP signal integration and targeting.  J Mol Biol  338:  271–285, 2004.  [Crossref]

Burzomato,  V, Beato,  M, Groot-Kormelink,  PJ, Colquhoun,  D & Sivilotti,  LG.  Single-channel behavior of heteromeric alpha1beta glycine receptors: an attempt to detect a conformational change before the channel opens.  J Neurosci  24:  10924–10940, 2004.  [Crossref]

Campo-Soria,  C, Chang,  Y & Weiss,  DS.  Mechanism of action of benzodiazepines on GABAA receptors.  Br J Pharmacol  148:  984–990, 2006.  [Crossref]

Carrijo,  JB, Antonio,  A & Rocha e Silva,  M.  On the nature of the antagonism atropine-acetylcholine on the guinea pig heart.  Acta Physiol Lat Am  27:  207–214, 1977.

Casadó,  V, Cortés,  A, Ciruela,  F, Mallol,  J, Ferré,  S, Lluis,  C, Canela,  EI & Franco,  R.  Old and new ways to calculate the affinity of agonists and antagonists interacting with G-protein-coupled monomeric and dimeric receptors: the receptor-dimer cooperativity index.  Pharmacol Ther  116:  343–354, 2007.  [Crossref]

Chang,  Y & Weiss,  DS.  Allosteric activation mechanism of the alpha1beta2gamma2 gamma-aminobutyric acid type A receptor revealed by mutation of the conserved M2 leucine.  Biophys J  77:  2542–2551, 1999.

Chen,  Y, Nong,  Y, Goudet,  C, Hemstapat,  K, de Paulis,  T, Pin,  JP & Conn,  PJ.  Interaction of novel positive allosteric modulators of metabotropic glutamate receptor 5 with the negative allosteric antagonist site is required for potentiation of receptor responses.  Mol Pharmacol  71:  1389–1398, 2007.  [Crossref]

Cherezov,  V, Rosenbaum,  DM, Hanson,  MA, Rasmussen,  SG, Thian,  FS, Kobilka,  TS, Choi,  HJ, Kuhn,  P, Weis,  WI, Kobilka,  BK & Stevens,  RC.   High-resolution srystal structure of an engineered human beta2-adrenergic G protein-coupled receptor science  318:  1258–1265, 2007.

Chidiac,  P.  Rethinking receptor-G protein-effector interactions.  Biochem Pharmacol  55:  549–556, 1998.  [Crossref]

Chidiac,  P, Nouet,  S & Bouvier,  M.  Agonist-induced modulation of inverse agonist efficacy at the beta 2-adrenergic receptor.  Mol Pharmacol  50:  662–669, 1996.

Childers,  SR, Li,  X, Xiao,  R & Eisenach,  JC.  Allosteric modulation of adenosine A1 receptor coupling to G-proteins in brain.  J Neurochem  93:  715–723, 2006.  [Crossref]

Christopoulos,  A & El-Fakahany,  EE.  Qualitative and quantitative assessment of relative agonist efficacy.  Biochem Pharmacol  58:  735–748, 1999.  [Crossref]

Christopoulos,  A & Kenakin,  T.  G protein coupled receptor allosterism and complexing.  Pharmacol Rev  54:  323–374, 2002.  [Crossref]

Christopoulos,  A, Lanzafame,  A & Mitchelson,  F.  Allosteric interactions at muscarinic cholinoceptors.  Clin Exp Pharmacol Physiol  25:  185–194, 1998.  [Crossref]

Christopoulos,  A, Grant,  MK & El-Fakahany,  EE.  Transducer abstraction: a novel approach to the detection of partial agonist efficacy in radioligand binding studies.  J Pharmacol Toxicol Methods  43:  55–67, 2000.  [Crossref]

Christopoulos,  A, Grant,  MK, Ayoubzadeh,  N, Kim,  ON, Sauerberg,  P, Jeppesen,  L & El-Fakahany,  EE.  Synthesis and pharmacological evaluation of dimeric muscarinic acetylcholine receptor agonists.  J Pharmacol Exp Ther  298:  1260–1268, 2001.

Colquhoun,  D.  Agonist-activated ion channels.  Br J Pharmacol  147:  S17–S26, 2006.  [Crossref]

Cornish-Bowden,  A    Fundamentals of Enzyme Kinetics .  London:  Portland Press, 1995.

Cornish-Bowden,  A.    Fundamentals of Enzyme Kinetics .  Colchester:  Portland Press Ltd., 2004.

Corsi,  M, Fina,  P & Trist,  DG.  Co-agonism in drug-receptor interaction: illustrated by the NMDA receptors.  Trends Pharmacol Sci  17:  220–222, 1996.  [Crossref]

Costa,  T, Ogino,  Y, Munson,  PJ, Onaran,  HO & Rodbard,  D.  Drug efficacy at guanine nucleotide-binding regulatory protein-linked receptors: thermodynamic interpretation of negative antagonism and of receptor activity in the absence of ligand.  Mol Pharmacol  41:  549–560, 1992.

Davila,  DF, Donis,  JH, Davila,  LA, Odreman,  WA, de Bellabarba,  GA & Villarroel,  V.  Anti-muscarinic autoantibodies and vagal modulation in Chagas disease: positive allosteric modulators vs desensitization and downregulation of M2 cardiac acetylcholine receptors.  Int J Cardiol  123:  328–329, 2008.  [Crossref]

Dawson,  RJP & Locher,  KP.  Structure of a bacterial multidrug ABC transporter.  Nature  443:  180–185, 2006.  [Crossref]

De Haën,  C.  The nonstoichiometric floating receptor model for hormone sensitive adenylyl cyclase.  J Theor Biol  58:  383–400, 1976.  [Crossref]

De Lean,  A, Stadel,  JM & Lefkowitz,  RJ.  A ternary complex model explains the agonist specific binding properties of the adenylate cyclase coupled beta-adrenergic receptor.  J Biol Chem  255:  7108–7117, 1980.

Dingledine,  R, Kleckner,  NW & McBain,  CJ.  The glycine coagonist site of the NMDA receptor.  Adv Exp Med Biol  268:  17–22, 1990.

Downing,  SS, Lee,  YT, Farb,  DH & Gibbs,  TT.  Benzodiazepine modulation of partial agonist efficacy and spontaneously active GABA(A) receptors supports an allosteric model of modulation.  Br J Pharmacol  145:  894–906, 2005.  [Crossref]

Dupre,  ML, Broyles,  JM & Mihic,  SJ.  Effects of a mutation in the TM2-TM3 linker region of the glycine receptor alpha1 subunit on gating and allosteric modulation.  Brain Res  1152C:  1–9, 2007.  [Crossref]

Durroux,  T.  Principles: a model for the allosteric interactions between ligand binding sites within a dimeric GPCR.  Trends Pharmacol Sci  26:  376–384, 2005.  [Crossref]

Dutzler,  R, Campbell,  EB, Chait,  BT & MacKinnon,  R.  X-ray structure of a ClC chloride channel at 3.0 A reveals the molecular basis of anion selectivity.  Nature  415:  287–294, 2002.  [Crossref]

Dutzler,  R, Campbell,  EB & MacKinnon,  R.  Gating the selectivity filter in ClC chloride channels.  Science  300:  108–112, 2003.  [Crossref]

Ehlert,  FJ.  Estimation of the affinities of allosteric ligands using radioligand binding and pharmacological null methods.  Mol Pharmacol  33:  187–194, 1988.

Ehlert,  FJ.  The ternary complex model.   Biomedical Application of Computer Modeling .  Boca Raton:  CRC Press, 2001.

Ehlert,  FJ.  Analysis of allosterism in functional assasys.  J Pharmacal ExpTher  315:  740–754, 2005.

Ehlert,  FJ.  On the analysis of ligand-directed signaling at G protein coupled receptors.   Naunyn Schmiedebergs Arch Pharmacol . 2008.

Ehlert,  FJ & Rathbun,  BE.  Signaling through the muscarinic receptor-adenylate cyclase system of the heart is buffered against GTP over a range of concentrations.  Mol Pharmacol  38:  148–158, 1990.

Ehlert,  FJ & &Griffin,  MT  Two-state Models and the Analysis of the Allosteric Effect of Gallamine at the M2 Muscarinic Receptor.   J Pharmacol Exp Ther . 2008.

Ferguson,  SS.  Phosphorylation-independent attenuation of GPCR signalling.  Trends Pharmacol Sci  28:  173–179, 2007.  [Crossref]

Finch,  EA, Turner,  TJ & Goldin,  SM.  Calcium as a coagonist of inositol 1,4,5-trisphosphate-induced calcium release.  Science  252:  443–446, 1991.  [Crossref]

Fisas,  A, Codony,  X, Romero,  G, Dordal,  A, Giraldo,  J, Merce,  R, Holenz,  J, Heal,  D, Buschmann,  H & Pauwels,  PJ.  Chronic 5-HT(6) receptor modulation by E-6837 induces hypophagia and sustained weight loss in diet-induced obese rats.  Br J Pharmacol  148:  973–983, 2006.  [Crossref]

Fitzsimons,  CP, Monczor,  F, Fernandez,  N, Shayo,  C & Davio,  C.  Mepyramine, a histamine H1 receptor inverse agonist, binds preferentially to a G protein-coupled form of the receptor and sequesters G protein.  J Biol Chem  279:  34431–34439, 2004.  [Crossref]

Fong,  TM.  Mechanistic hypotheses for the activation of G-protein-coupled receptors.  Cell Signal  8:  217–224, 1996.  [Crossref]

Foster,  AC & Kemp,  JA.  Glutamate- and GABA-based CNS therapeutics.  Curr Opin Pharmacol  6:  7–17, 2006.  [Crossref]

Franco,  R, Casado,  V, Mallol,  J, Ferre,  S, Fuxe,  K, Cortes,  A, Ciruela,  F, Lluis,  C & Canela,  EI.  Dimer-based model for heptaspanning membrane receptors.  Trends Biochem Sci  30:  360–366, 2005.  [Crossref]

Franco,  R, Casado,  V, Mallol,  J, Ferrada,  C, Ferre,  S, Fuxe,  K, Cortes,  A, Ciruela,  F, Lluis,  C & Canela,  EI.  The two-state dimer receptor model: a general model for receptor dimers.  Mol Pharmacol : –1912  69:  2006

1905.

Franco,  R, Casado,  V, Mallol,  J, Ferrada,  C, Ferre,  S, Fuxe,  K, Cortes,  A, Ciruela,  F, Lluis,  C & Canela,  EI.  G-protein-coupled receptor heteromers: function and ligand pharmacology.  Br J Pharmacol  153:  S90–S98, 2007a.  [Crossref]

Franco,  R, Casadó,  V, Cortés,  A, Ferrada,  C, Mallol,  J, Woods,  A, Lluis,  C, Canela,  EI & Ferré,  S.   Basic concepts in G-Protein-coupled receptor homo- and heterodimerization ScientificWorldJournal  7:  47–57, 2007b.

Francotte,  P, Tullio,  P, Goffin,  E, Dintilhac,  G, Graindorge,  E, Fraikin,  P, Lestage,  P, Danober,  L, Thomas,  JY, Caignard,  DH & Pirotte,  B.  Design, synthesis, and pharmacology of novel 7-substituted 3,4-dihydro-2H-1,2,4-benzothiadiazine 1,1-dioxides as positive allosteric modulators of AMPA receptors.  J Med Chem  50:  3153–3157, 2007.  [Crossref]

Frieden,  C.  Kinetic aspects of regulation of metabolic processes. The hysteretic enzyme concept.  J Biol Chem  245:  5788–5799, 1970.

Gao,  ZG & Jacobson,  KA.  Keynote review: allosterism in membrane receptors.  Drug Discov Today  11:  191–202, 2006.  [Crossref]

Ghosh,  P, Moitra,  K, Maki,  N & Dey,  S.  Allosteric modulation of the human P-glycoprotein involves conformational changes mimicking catalytic transition intermediates.  Arch Biochem Biophys  450:  100–112, 2006.  [Crossref]

Giraldo,  J.  Agonist induction, conformational selection, and mutant receptors.  FEBS Lett  256:  13–18, 2004.  [Crossref]

Giraldo,  J, Roche,  D, Rovira,  X & Serra,  J.  The catalytic power of enzymes: conformational selection or transition state stabilization?.  FEBS Lett  580:  2170–2177, 2006.  [Crossref]

Giraldo,  J, Serra,  J, Roche,  D & Rovira,  X.  Assessing receptor affinity for inverse agonists: Schild and Cheng-Prusoff methods revisited.  Curr Drug Targets  8:  197–202, 2007.  [Crossref]

Gouaux,  E & MacKinnon,  R.  Principles of selective ion transport in channels and pumps.  Science  310:  1461–1465, 2005.  [Crossref]

Gronlien,  JH, Haakerud,  M, Ween,  H, Thorin-Hagene,  K, Briggs,  CA & Gopalakrishnan,  M  Distinct profiles of {alpha}7 nAChR positive allosteric modulation revealed by structurally diverse chemotypes.  Mol Pharmacol  72:  715–724, 2007.  [Crossref]

Haber,  JE & Koshland,  DE jr  Relation of protein subunit interactions to the molecular species observed during cooperative binding of ligands.  Proc Natl Acad Sci USA  58:  2087–2093, 1967.  [Crossref]

Hall,  DA.  Modeling the functional effects of allosteric modulators at pharmacological receptors: an extension of the two-state model of receptor activation.  Mol Pharmacol  58:  1412–1423, 2000.

Hay,  DL, Christopoulos,  G, Christopoulos,  A, Poyner,  DR & Sexton,  PM.  Pharmacological discrimination of calcitonin receptor: receptor activity-modifying protein complexes.  Mol Pharmacol  67:  1655–1665, 2005.  [Crossref]

Hemstapat,  K, Da Costa,  H, Nong,  Y, Brady,  AE, Luo,  Q & Niswender,  CM  A novel family of potent negative allosteric modulators of group II metabotropic glutamate receptors.  J Pharmacol Exp Ther  322:  254–264, 2007.  [Crossref]

Himoe,  A.  Parks PG & Hess Gp. Investigations of the chymotrypsincatalyzed hydrolysis of specific substrates I. @ The pH dependence of the catalytic hydrolysis of N-acetyl-L-tryptophanamide by three forms of the enzyme at alkaline pH.  J Biol Chem  242:  919–929, 1967.

Hoare,  SR, Fleck,  BA, Gross,  RS, Crowe,  PD, Williams,  JP & Grigoriadis,  DE.  Allosteric ligands for the corticotropin releasing factor type 1 receptor modulate conformational states involved in receptor activation.   Mol Pharmacol . 2008.

Hogenkamp,  , Johnstone,  TB, Huang,  JC, Li,  WY, Tran,  M, Whittemore,  ER, Bagnera,  RE & Gee,  KW.  Enaminone amides as novel orally active GABA(A) receptor modulators.  J Med Chem  50:  3369–3379, 2007.  [Crossref]

Hogg,  RC, Buisson,  B & Bertrand,  D.  Allosteric modulation of ligand-gated ion channels.  Biochem Pharmacol  70:  1267–1276, 2005.  [Crossref]

Holst,  B, Brandt,  E, Bach,  A, Heding,  A & Schwartz,  TW.  Nonpeptide and peptide growth hormone secretagogues act both as ghrelin receptor agonist and as positive or negative allosteric modulators of ghrelin signaling.  Mol Endocrinol  19:  2400–2411, 2005.  [Crossref]

Holzgrabe,  U, De Amici,  M & Mohr,  K.  Allosteric modulators and selective agonists of muscarinic receptors.  J Mol Neurosci  30:  165–168, 2006.  [Crossref]

Hornigold,  DC, Mistry,  R, Raymond,  PD, Blank,  JL & Challiss,  RA.  Evidence for cross-talk between M2 and M3 muscarinic acetylcholine receptors in the regulation of second messenger and extracellular signal-regulated kinase signalling pathways in Chinese hamster ovary cells.  Br J Pharmacol  138:  1340–1350, 2003.  [Crossref]

Huang,  XP & Ellis,  J.  Mutational disruption of a conserved disulfide bond in muscarinic acetylcholine receptors attenuates positive homotropic cooperativity between multiple allosteric sites and has subtype-dependent effects on the affinities of muscarinic allosteric ligands.  Mol Pharmacol  71:  759–768, 2007.  [Crossref]

Hulme,  EC, Birdsall,  NJ, Burgen,  AS & Mehta,  P.  The binding of antagonists to brain muscarinic receptors.  Mol Pharmacol  14:  737–750, 1978.

Hunte,  C, Screpanti,  E, Venturi,  M, Rimon,  A, Padan,  E & Michel,  H.  Structure of a Na + /H+ antiporter and insights into mechanism of action and regulation by pH.  Nature  435:  1197–1202, 2005.  [Crossref]

Indarte,  M, Madura,  JD & Surratt,  CK.  Dopamine transporter comparative molecular modeling and binding site prediction using the LeuT(Aa) leucine transporter as a template.  Proteins  70:  1033–1046, 2007.  [Crossref]

Jacobs,  S & Cuatrecasas,  P.  The mobile receptor hypothesis and ‘cooperativity’ of hormone binding. Application to insulin.  Biochim Biophys Acta  433:  482–495, 1976.  [Crossref]

Jacobson,  KA & Gao,  ZG.  Adenosine receptors as therapeutic targets.  Nat Rev Drug Discov  5:  247–264, 2006.  [Crossref]

Jakubik,  J, Bacakova,  L, Lisa,  V, El-Fakahany,  EE & Tucek,  S.  Activation of muscarinic acetylcholine receptors via their allosteric binding sites.  Proc Natl Acad Sci USA  93:  8705–8709, 1996.  [Crossref]

Jasti,  J, Furukawa,  H, Gonzales,  EB & Gouaux,  E.  Structure of acid-sensing ion channel 1 at 1.9 A resolution and low pH.  Nature  449:  316–323, 2007.  [Crossref]

Jensen,  AA & Spalding,  TA.  Allosteric modulation of G-protein coupled receptors.  Eur J Pharm Sci  21:  407–420, 2004.  [Crossref]

Jensen,  AM, Sorensen,  TL, Olesen,  C, Moller,  JV & Nissen,  P.  Modulatory and catalytic modes of ATP binding by the calcium pump.  EMBO J  25:  2305–2314, 2006.  [Crossref]

Jentsch,  TJ, Neagoe,  I & Scheel,  O.  CLC chloride channels and transporters.  Curr Opin Neurobiol  15:  319–325, 2005a.  [Crossref]

Jentsch,  TJ, Poet,  M, Fuhrmann,  JC & Zdebik,  AA.  Physiological functions of CLC Cl-channels gleaned from human genetic disease and mouse models.  Annu Rev Physiol  67:  779–807, 2005b.  [Crossref]

Johnson,  JW & Ascher,  P.  Glycine potentiates the NMDA response in cultured mouse brain neurons.  Nature  325:  529–531, 1987.  [Crossref]

Jones,  SV, Heilman,  CJ & Brann,  MR.  Functional responses of cloned muscarinic receptors expressed in CHO-K1 cells.  Mol Pharmacol  40:  242–247, 1991.

Kashihara,  K, Varga,  EV, Waite,  SL, Roeske,  WR & Yamamura,  HI.  Cloning of the rat M3, M4 and M5 muscarinic acetylcholine receptor genes by the polymerase chain reaction (PCR) and the pharmacological characterization of the expressed genes.  Life Sci  51:  955–971, 1992.  [Crossref]

Katz,  B & Thesleff,  S.  A study of the desensitization produced by acetylcholine at the motor end-plate.  J Physiol  138:  63–80, 1957.

Kenakin,  T.  Differences between natural and recombinant G protein-coupled receptor systems with varying receptor/G protein stoichiometry.  Trends Pharmacol Sci  18:  456–464, 1997.

Kew,  JN.  Positive and negative allosteric modulation of metabotropic glutamate receptors: emerging therapeutic potential.  Pharmacol Ther  104:  233–244, 2004.  [Crossref]

Knudsen,  LB, Kiel,  D, Teng,  M, Behrens,  C, Bhumralkar,  D, Kodra,  JT, Holst,  JJ, Jeppesen,  CB, Johnson,  MD, de Jong,  JC, Jorgensen,  AS, Kercher,  T, Kostrowicki,  J, Madsen,  P, Olesen,  PH, Petersen,  JS, Poulsen,  F, Sidelmann,  UG, Sturis,  J, Truesdale,  L, May,  J & Lau,  J.  Small-molecule agonists for the glucagon-like peptide 1 receptor.  Proc Natl Acad Sci USA  104:  937–942, 2007.  [Crossref]

Koshland,  DE Jr, Nemethy,  G & Filmer,  D.  Comparison of experimental binding data and theoretical models in proteins containing subunits.  Biochemistry  5:  365–385, 1966.  [Crossref]

Kukkonen,  JP.  Explicit formulation of different receptor-G-protein interactions and effector regulation.  Bioinformatics  20:  2411–2420, 2004a.  [Crossref]

Kukkonen,  JP.  Regulation of receptor-coupling to (multiple) G proteins. A challenge for basic research and drug discovery.  Receptors Channels  10:  167–183, 2004b.

Kukkonen,  JP, Nasman,  J & Akerman,  KE.  Modelling of promiscuous receptor-Gi/Gs-protein coupling and effector response.  Trends Pharmacol Sci  22:  616–622, 2001.  [Crossref]

Kurganov,  BI.    Allosteric enzymes. Kinetic Behaviour .  Chichester:  Wiley & Sons, 1982.

Langmead,  CJ.  Screening for positive allosteric modulators: assessment of modulator concentration-response curves as a screening paradigm.  J Biomol Screen  12:  668–676, 2007.  [Crossref]

Langmead,  CJ & Christopoulos,  A.  Allosteric agonists of 7TM receptors: expanding the pharmacological toolbox.  Trends Pharmacol Sci  27:  475–581, 2006.  [Crossref]

Langmead,  , Fry,  VA, Forbes,  IT, Branch,  CL, Christopoulos,  A, Wood,  MD & Herdon,  HJ.  Probing the molecular mechanism of interaction between 4-n-butyl-1-[4-(2-methylphenyl)-4-oxo-1-butyl]-piperidine (AC-42) and the muscarinic M(1) receptor: direct pharmacological evidence that AC-42 is an allosteric agonist.  Mol Pharmacol  69:  236–246, 2006.

Lanzafame,  AA, Sexton,  PM & Christopoulos,  A.  Interaction studies of multiple binding sites on M4 muscarinic acetylcholine receptors.  Mol Pharmacol  70:  736–746, 2006.  [Crossref]

Lazareno,  S & Birdsall,  NJ.  Detection, quantitation, and verification of allosteric interactions of agents with labeled and unlabeled ligands at G protein-coupled receptors: interactions of strychnine and acetylcholine at muscarinic receptors.  Mol Pharmacol  48:  362–378, 1995.

Lazareno,  S, Dolezal,  V, Popham,  A & Birdsall,  NJ.  Thiochrome enhances acetylcholine affinity at muscarinic M4 receptors: receptor subtype selectivity via cooperativity rather than affinity.  Mol Pharmacol  65:  257–266, 2004.  [Crossref]

Lecourtier,  L, Homayoun,  H, Tamagnan,  G & Moghaddam,  B.  Positive allosteric modulation of metabotropic glutamate 5 (mGlu5) receptors reverses N-methyl-d-aspartate antagonist-induced alteration of neuronal firing in prefrontal cortex.  Biol Psychiatry  62:  739–746, 2007.  [Crossref]

Leff,  P.  The two state model of receptor activation.  Trends Pharmacol Sci  16:  89–97, 1995.  [Crossref]

Leff,  P, Scaramellini,  C, Law,  C & McKechnie,  K.  A three-state receptor model of agonist action.  Trends Pharmacol Sci  18:  355–362, 1997.

Lefkowitz,  RJ, Cotecchia,  S, Samama,  P & Costa,  T.  Constitutive activity of receptors coupled to guanine nucleotide regulatory proteins.  Trends Pharmacol Sci  14:  303–307, 1993.  [Crossref]

Litschig,  S, Gasparini,  F, Rueegg,  D, Stoehr,  N, Flor,  PJ, Vranesic,  I, Prezeau,  L, Pin,  JP, Thomsen,  C & Kuhn,  R.  CPCCOEt, a noncompetitive metabotropic glutamate receptor 1 antagonist, inhibits receptor signaling without affecting glutamate binding.  Mol Pharmacol  55:  453–461, 1999.

MacKinnon,  R.  Potassium channels.  FEBS Lett  555:  62–65, 2003.  [Crossref]

MacKinnon,  R.  Nobel Lecture. Potassium channels and the atomic basis of selective ion conduction.  Biosci Rep  24:  75–100, 2004.  [Crossref]

Maillet,  EL, Pellegrini,  N, Valant,  C, Bucher,  B, Hibert,  M, Bourguignon,  JJ & Galzi,  JL.  A novel, conformation-specific allosteric inhibitor of the tachykinin NK2 receptor (NK2R) with functionally selective properties.  FASEB J  21:  2124–2134, 2007.  [Crossref]

Mak,  DO, McBride,  SM & Foskett,  JK.  Spontaneous channel activity of the inositol 1,4,5-trisphosphate (InsP3) receptor (InsP3R). Application of allosteric modeling to calcium and InsP3 regulation of InsP3R single-channel gating.  J Gen Physiol  122:  583–603, 2003.  [Crossref]

Maki,  N, Moitra,  K, Ghosh,  P & Dey,  S.  Allosteric modulation bypasses the requirement for ATP hydrolysis in regenerating low affinity transition state conformation of human P-glycoprotein.  J Biol Chem  281:  10769–10777, 2006.  [Crossref]

Marvizon,  JC & Baudry,  M.  Receptor activation by two agonists: analysis by nonlinear regression and application to N-methyl-d-aspartate receptors.  Anal Biochem  213:  3–11, 1993.  [Crossref]

May,  LT, Avlani,  VA, Sexton,  PM & Christopoulos,  A.  Allosteric modulation of G protein-coupled receptors.  Curr Pharm Des : –2013  10:  2004

2003.

May,  LT, Leach,  K, Sexton,  PM & Christopoulos,  A.  Allosteric modulation of G protein-coupled receptors.  Annu Rev Pharmacol Toxicol  47:  1–51, 2007.  [Crossref]

Mehta,  AK & Ticku,  MK.  An update on GABAA receptors.  Brain Res Brain Res Rev  29:  196–217, 1999.  [Crossref]

Michal,  P, Lysikova,  M & Tucek,  S.  Dual effects of muscarinic M(2) acetylcholine receptors on the synthesis of cyclic AMP in CHO cells: dependence on time, receptor density and receptor agonists.  Br J Pharmacol  132:  1217–1228, 2001.  [Crossref]

Michelsen,  S, Sanchez,  C & Ebert,  B.  Lack of generalisation between the GABA(A) receptor agonist, gaboxadol, and allosteric modulators of the benzodiazepine binding site in the rat drug discrimination procedure.  Psychopharmacology (Berl)  193:  151–157, 2007.  [Crossref]

Migeon,  JC & Nathanson,  NM.  Differential regulation of cAMP-mediated gene transcription by m1 and m4 muscarinic acetylcholine receptors. Preferential coupling of m4 receptors to Gi alpha-2.  J Biol Chem  269:  9767–9773, 1994.

Miller,  C.  ClC chloride channels viewed through a transporter lens.  Nature  440:  484–489, 2006.  [Crossref]

Monczor,  F, Fernandez,  N, Legnazzi,  BL, Riveiro,  ME, Baldi,  A, Shayo,  C & Davio,  C.  Tiotidine, a histamine H2 receptor inverse agonist that binds with high affinity to an inactive G-protein-coupled form of the receptor. Experimental support for the cubic ternary complex model.  Mol Pharmacol  64:  512–520, 2003.  [Crossref]

Monod,  J, Wyman,  J & Changeux,  J-P.  On the nature of allosteric transitions: a plausible model.  J Mol Biol  12:  88–118, 1965.

Murakami,  S, Nakashima,  R, Yamashita,  E, Matsumoto,  T & Yamaguchi,  A.  Crystal structures of a multidrug transporter reveal a functionally rotating mechanism.  Nature  443:  173–179, 2006.  [Crossref]

Nasman,  J, Kukkonen,  JP, Ammoun,  S & Akerman,  KE.  Role of G-protein availability in differential signaling by alpha 2-adrenoceptors.  Biochem Pharmacol  62:  913–922, 2001.  [Crossref]

Neubig,  RR.  Membrane organization in G-protein mechanisms.  FASEB J  8:  939–946, 1994.

Nguitragool,  W & Miller,  C.  Uncoupling of a CLC Cl(–)/H(+) exchange transporter by polyatomic anions.  J Mol Biol  362:  682–690, 2006.  [Crossref]

Nicoll,  DA, Sawaya,  MR, Kwon,  S, Cascio,  D, Philipson,  KD & Abramson,  J.  The crystal structure of the primary Ca2+ sensor of the Na + /Ca2+ exchanger reveals a novel Ca2+ binding motif.  J Biol Chem  281:  21577–21581, 2006.  [Crossref]

O'Brien,  JA, Lemaire,  W, Wittmann,  M, Jacobson,  MA, Ha,  SN, Wisnoski,  DD, Lindsley,  CW, Schaffhauser,  HJ, Rowe,  B, Sur,  C, Duggan,  ME, Pettibone,  DJ, Conn,  PJ & Williams,  DL.  A novel selective allosteric modulator potentiates the activity of native metabotropic glutamate receptor subtype 5 in rat forebrain.  J Pharmacol Exp Ther  309:  568–577, 2004.  [Crossref]

Obara,  K, Miyashita,  N, Xu,  C, Toyoshima,  I, Sugita,  Y & Inesi,  G  Structural role of countertransport revealed in Ca(2 + ) pump crystal structure in the absence of Ca(2 + ).  Proc Natl Acad Sci USA  102:  14489–14496, 2005.  [Crossref]

Okada,  T, Fujiyoshi,  Y, Silow,  M, Navarro,  J, Landau,  EM & Shichida,  Y.  Functional role of internal water molecules in rhodopsin revealed by X-ray crystallography.  Proc Natl Acad Sci USA  99:  5982–5987, 2002.  [Crossref]

Okada,  T, Sugihara,  M, Bondar,  AN, Elstner,  M, Entel,  P & Buss,  V.  The retinal conformation and its environment in rhodopsin in light of a new 2.2 A crystal structure.  J Mol Biol  342:  571–583, 2004.  [Crossref]

Oldham,  WM, Van Eps,  N, Preininger,  AM, Hubbell,  WL & Hamm,  HE.  Mechanism of the receptor-catalyzed activation of heterotrimeric G proteins.  Nat Struct Mol Biol  13:  772–777, 2006.  [Crossref]

Olsen,  RW, Chang,  CS, Li,  G, Hanchar,  HJ & Wallner,  M.  Fishing for allosteric sites on GABA(A) receptors.  Biochem Pharmacol  68:  1675–1684, 2004.  [Crossref]

Onaran,  HO, Costa,  T & Rodbard,  D.  Beta gamma subunits of guanine nucleotide-binding proteins and regulation of spontaneous receptor activity: thermodynamic model for the interaction between receptors and guanine nucleotide-binding protein subunits.  Mol Pharmacol  43:  245–256, 1993.

Ong,  J & Kerr,  DI.  Clinical potential of GABAB receptor modulators.  CNS Drug Rev  11:  317–334, 2005.

Palczewski,  K.  G protein-coupled receptor rhodopsin.  Annu Rev Biochem  75:  743–767, 2006.  [Crossref]

Palczewski,  K, Kumasaka,  T, Hori,  T, Behnke,  CA, Motoshima,  H, Fox,  BA, Le Trong,  I, Teller,  DC, Okada,  T, Stenkamp,  RE, Yamamoto,  M & Miyano,  M.  Crystal structure of rhodopsin: A G protein-coupled receptor.  Science  289:  739–745, 2000.  [Crossref]

Parmentier,  ML, Prezeau,  L, Bockaert,  J & Pin,  JP.  A model for the functioning of family 3 GPCRs.  Trends Pharmacol Sci  23:  268–274, 2002.  [Crossref]

Pin,  JP, Kniazeff,  J, Liu,  J, Binet,  V, Goudet,  C, Rondard,  P & Prezeau,  L.  Allosteric functioning of dimeric class C G-protein-coupled receptors.  FEBS J  272:  2947–2955, 2005.  [Crossref]

Pineyro,  G, Azzi,  M, deLean,  A, Schiller,  PW & Bouvier,  M.  Reciprocal regulation of agonist and inverse agonist signaling efficacy upon short-term treatment of the human delta-opioid receptor with an inverse agonist.  Mol Pharmacol  67:  336–348, 2005.  [Crossref]

Plested,  AJ, Groot-Kormelink,  PJ, Colquhoun,  D & Sivilotti,  LG.  Single-channel study of the spasmodic mutation alpha1A52S in recombinant rat glycine receptors.  J Physiol  581:  51–73, 2007.  [Crossref]

Proska,  J & Tucek,  S.  Mechanisms of steric and cooperative actions of alcuronium on cardiac muscarinic acetylcholine receptors.  Mol Pharmacol  45:  709–717, 1994.

Proska,  J & Tucek,  S.  Competition between positive and negative allosteric effectors on muscarinic receptors.  Mol Pharmacol  48:  696–702, 1995.

Raddatz,  R, Schaffhauser,  H & Marino,  MJ.  Allosteric approaches to the targeting of G-protein-coupled receptors for novel drug discovery: a critical assessment.  Biochem Pharmacol  74:  383–391, 2007.  [Crossref]

Ramjeesingh,  M, Li,  C, She,  YM & Bear,  CE.  Evaluation of the membrane-spanning domain of ClC-2.  Biochem J  396:  449–460, 2006.  [Crossref]

Ritzen,  A, Mathiesen,  JM & Thomsen,  C.  Molecular pharmacology and therapeutic prospects of metabotropic glutamate receptor allosteric modulators.  Basic Clin Pharmacol Toxicol  97:  202–213, 2005.  [Crossref]

Ross,  RA.  Tuning the endocannabinoid system: allosteric modulators of the CB(1) receptor.  Br J Pharmacol  152:  565–566, 2007.  [Crossref]

Ross,  RA.  Allosterism and cannabinoid CB(1) receptors: the shape of things to come.  Trends Pharmacol Sci  28:  567–572, 2007.  [Crossref]

Rovira,  X, Roche,  D, Serra,  J, Kniazeff,  J, Pin,  JP & Giraldo,  J.  Modeling the Binding and Function of Metabotropic Glutamate Receptors.   J Pharmacol Exp Ther Feb 20; [Epub ahead of print]2008.

Rusch,  D & Forman,  SA.  Classic benzodiazepines modulate the open-close equilibrium in alpha1beta2gamma2L gamma-aminobutyric acid type A receptors.  Anesthesiology  102:  783–792, 2005.  [Crossref]

Rusch,  D, Zhong,  H & Forman,  SA.  Gating allosterism at a single class of etomidate sites on alpha1beta2gamma2L GABA A receptors accounts for both direct activation and agonist modulation.  J Biol Chem  279:  20982–20992, 2004.  [Crossref]

Ryder,  JW, Falcone,  JF, Manro,  JR, Svensson,  KA & Merchant,  KM.  Pharmacological characterization of cGMP regulation by the biarylpropylsulfonamide class of positive, allosteric modulators of alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors.  J Pharmacol Exp Ther  319:  293–298, 2006.  [Crossref]

Samama,  P, Cotecchia,  S, Costa,  T & Lefkowitz,  RJ.  A mutation induced activated state of the beta 2adrenergic receptor. Extending the ternary complex model.  J Biol Chem  268:  4625–4636, 1993.

Scaramellini,  C & Leff,  P.  A three-state receptor model: predictions of multiple agonist pharmacology for the same receptor type.  Ann N Y Acad Sci  861:  97–103, 1998.  [Crossref]

Scaramellini,  C & Leff,  P.  Theoretical implications of receptor coupling to multiple G proteins based on analysis of a three-state model.  Methods Enzymol  343:  17–29, 2002.

Scheer,  A, Fanelli,  F, Costa,  T, De Benedetti,  PG & Cotecchia,  S.  Constitutively active mutants of the alpha 1B-adrenergic receptor: role of highly conserved polar amino acids in receptor activation.  EMBO J  15:  3566–3578, 1996.

Schetz,  JA.  Allosteric modulation of dopamine receptors.  Mini Rev Med Chem  5:  555–561, 2005.  [Crossref]

Shipe,  WD, Wolkenberg,  SE, Williams,  DL Jr & Lindsley,  CW.  Recent advances in positive allosteric modulators of metabotropic glutamate receptors.  Curr Opin Drug Discov Dev  8:  449–457, 2005.

Solt,  K, Ruesch,  D, Forman,  SA, Davies,  PA & Raines,  DE.  Differential effects of serotonin and dopamine on human 5-HT3A receptor kinetics: interpretation within an allosteric kinetic model.  J Neurosci  27:  13151–13160, 2007.  [Crossref]

Spalding,  TA, Trotter,  C, Skjaerbaek,  N, Messier,  TL, Currier,  EA, Burstein,  ES, Li,  D, Hacksell,  U & Brann,  MR.  Discovery of an ectopic activation site on the M(1) muscarinic receptor.  Mol Pharmacol  61:  1297–1302, 2002.  [Crossref]

Teng,  M, Johnson,  C, Thomas,  MD, Kiel,  D, Lakis,  JN, Kercher,  T, Aytes,  S, Kostrowicki,  J, Bhumralkar,  D, Truesdale,  L, May,  J, Sidelman,  U, Kodra,  JT, Jørgensen,  AS, Olesen,  PH, de Jong,  JC, Madsen,  P, Behrens,  C, Pettersson,  I, Knudsen,  LB, Holst,  JJ & Lau,  J.  Small molecule ago-allosteric modulators of the human glucagon-like peptide-1 (hGLP-1) receptor.  Bioorg Med Chem Lett  17:  5472–5478, 2007.  [Crossref]

Timmermann,  DB, Grønlien,  JH, Kohlhaas,  KL, Nielsen,  EØ, Dam,  E, Jørgensen,  TD, Ahring,  PK, Peters,  D, Holst,  D, Chrsitensen,  JK, Malysz,  J, Briggs,  CA, Gopalakrishnan,  M & Olsen,  GM.  An allosteric modulator of the alpha7 nicotinic acetylcholine receptor possessing cognition-enhancing properties in vivo.  J Pharmacol Exp Ther  323:  294–307, 2007.  [Crossref]

Tolkovsky,  AM & Levitzki,  A.  Mode of coupling between the betaadrenergic receptor and adenylate cyclase in turkey erythrocytes.  Biochemistry  17:  3795–3810, 1978.  [Crossref]

Tolkovsky,  AM & Levitzki,  A.  Theories and predictions of models describing sequential interactions between the receptor, the GTP regulatory unit, and the catalytic unit of hormone dependent adenylate cyclases.  J Cyclic Nucleotide Res  7:  139–150, 1981.

Toyoshima,  C.  Ion pumping by calcium ATPase of sarcoplasmic reticulum..   Regulatory mechanisms of striated muscle contraction vol 592, Adv Exp Med Biol Ebashi S Tokoyo, Springer Verlag2007.

Toyoshima,  C, Nakasako,  M, Nomura,  H & Ogawa,  H.  Crystal structure of the calcium pump of sarcoplasmic reticulum at 2.6 A resolution.  Nature  405:  647–655, 2000.  [Crossref]

Tucek,  S & Proska,  J.  Allosteric modulation of muscarinic acetylcholine receptors.  Trends Pharmacol Sci  16:  205–212, 1995.  [Crossref]

Tucek,  S, Michal,  P & Vlachova,  V.  Dual effects of muscarinic M2 receptors on the synthesis of cyclic AMP in CHO cells: background and model.  Life Sci  68:  2501–2510, 2001.  [Crossref]

Tucek,  S, Michal,  P & Vlachova,  V.  Modelling the consequences of receptor-G-protein promiscuity.  Trends Pharmacol Sci  23:  171–176, 2002.  [Crossref]

Unwin,  N.  Refined structure of the nicotinic acetylcholine receptor at 4A resolution.  J Mol Biol  346:  967–989, 2005.  [Crossref]

van Rijn,  CM & Willems-van Bree,  E.  A four-ligand hypercube model to quantify allosteric interactions within the GABAA receptor complex.  Eur J Pharmacol  485:  43–51, 2004.  [Crossref]

Vauquelin,  G & Van Liefde,  I.  G protein-coupled receptors: a count of 1001 conformations.  Fundam Clin Pharmacol  19:  45–56, 2005.  [Crossref]

Vieira,  E, Huwyler,  J, Jolidon,  S, Knoflach,  F, Mutel,  V & Wichmann,  J.  9H-Xanthene-9-carboxylic acid [1,2,4]oxadiazol-3-yl- and (2H-tetrazol-5-yl)-amides as potent, orally available mGlu1 receptor enhancers.  Bioorg Med Chem Lett  15:  4628–4631, 2005.  [Crossref]

Vogel,  WK, Mosser,  VA, Bulseco,  DA & Schimerlik,  MI.  Porcine m2 muscarinic acetylcholine receptor-effector coupling in Chinese hamster ovary cells.  J Biol Chem  270:  15485–15493, 1995.  [Crossref]

Waelbroeck,  M.  Identification of drugs competing with d-tubocurarine for an allosteric site on cardiac muscarinic receptors.  Mol Pharmacol  46:  685–692, 1994.

Weiss,  JM, Morgan,  PH, Lutz,  MW & Kenakin,  TP.  The cubic ternary complex receptor occupancy model. III. Resurrecting efficacy.  J Theor Biol  181:  381–397, 1996a.  [Crossref]

Weiss,  JM, Morgan,  PH, Lutz,  MW & Kenakin,  TP.  The cubic ternary complex receptor-occupancy model II. Understanding apparent affinity.  J Theor Biol  178:  169–182, 1996b.  [Crossref]

Weiss,  JM, Morgan,  PH, Lutz,  MW & Kenakin,  TP.  The cubic ternary complex receptor-occupancy model I. Model description.  J Theor Biol  178:  151–167, 1996c.  [Crossref]

Wells,  JW.  Analysis and interpretation of binding at equilibrium.   Receptor–ligand Interactions. A Practical Approach .  Oxford:  IRL Press at Oxford University Press, 1992.

Winding,  B & Bindslev,  N.  Desensitization and reactivation of ACh-regulated exocrine secretion in hen tracheal epithelium.  Am J Physiol  264:  C342–C351, 1993.

Wreggett,  KA & Wells,  JW.  Cooperativity manifest in the binding properties of purified cardiac muscarinic receptors.  J Biol Chem  270:  22488–22499, 1995.  [Crossref]

Yamashita,  A, Singh,  SK, Kawate,  T, Jin,  Y & Gouaux,  E.  Crystal structure of a bacterial homologue of Na + /Cl–-dependent neurotransmitter transporters.  Nature  437:  215–223, 2005.  [Crossref]

Yernool,  D, Boudker,  O, Jin,  Y & Gouaux,  E.  Structure of a glutamate transporter homologue from Pyrococcus horikoshii.  Nature  431:  811–818, 2004.  [Crossref]

Yoshimura,  RF, Hogenkamp,  DJ, Li,  WY, Tran,  MB, Belluzzi,  JD, Whittemore,  ER, Leslie,  FM & Gee,  KW.  Negative allosteric modulation of nicotinic acetylcholine receptors blocks nicotine self-administration in rats.  J Pharmacol Exp Ther  323:  907–915, 2007.  [Crossref]

Yu,  EW, Aires,  JR & Nikaido,  H.  AcrB multidrug efflux pump of Escherichia coli: composite substrate-binding cavity of exceptional flexibility generates its extremely wide substrate specificity.  J Bacteriol  185:  5657–5664, 2003.  [Crossref]



Drug-acceptor interactions ISBN 978-91-977071-0-7 (print), 978-91-977071-1-4 (online)

This book is free to use under a Creative Commons Attribution-Noncommercial 3.0 Unported License.