Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET

User menu

  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Molecular Pharmacology
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET
  • My alerts
  • Log in
  • My Cart
Molecular Pharmacology

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Visit molpharm on Facebook
  • Follow molpharm on Twitter
  • Follow molpharm on LinkedIn
Research ArticleArticle

Evaluation of Operational Models of Agonism and Allosterism at Receptors with Multiple Orthosteric Binding Sites

Karen J. Gregory, Jesús Giraldo, Jiayin Diao, Arthur Christopoulos and Katie Leach
Molecular Pharmacology January 2020, 97 (1) 35-47; DOI: https://doi.org/10.1124/mol.119.118091
Karen J. Gregory
Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Melbourne, Australia (K.J.G., J.D., A.C., K.L.); Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Barcelona, Spain (J.G.); Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, Bellaterra, Spain (J.G.); and Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí and Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra, Spain (J.G.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jesús Giraldo
Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Melbourne, Australia (K.J.G., J.D., A.C., K.L.); Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Barcelona, Spain (J.G.); Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, Bellaterra, Spain (J.G.); and Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí and Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra, Spain (J.G.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jiayin Diao
Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Melbourne, Australia (K.J.G., J.D., A.C., K.L.); Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Barcelona, Spain (J.G.); Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, Bellaterra, Spain (J.G.); and Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí and Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra, Spain (J.G.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Arthur Christopoulos
Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Melbourne, Australia (K.J.G., J.D., A.C., K.L.); Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Barcelona, Spain (J.G.); Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, Bellaterra, Spain (J.G.); and Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí and Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra, Spain (J.G.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Katie Leach
Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Melbourne, Australia (K.J.G., J.D., A.C., K.L.); Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Barcelona, Spain (J.G.); Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, Bellaterra, Spain (J.G.); and Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí and Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra, Spain (J.G.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF
Loading

This article has a correction. Please see:

  • Correction to “Evaluation of Operational Models of Agonism and Allosterism at Receptors with Multiple orthosteric Binding Sites” - January 01, 2020

Abstract

Current operational models of agonism and allosterism quantify ligand actions at receptors where agonist concentration-response relationships are nonhyperbolic by introduction of a transducer slope that relates receptor occupancy to response. However, for some receptors nonhyperbolic concentration-response relationships arise from multiple endogenous agonist molecules binding to a receptor in a cooperative manner. Thus, we developed operational models of agonism in systems with cooperative agonist binding and evaluated the models by simulating data describing agonist effects. The models were validated by analyzing experimental data demonstrating the effects of agonists and allosteric modulators at receptors where agonist binding follows hyperbolic (M4 muscarinic acetylcholine receptors) or nonhyperbolic relationships (metabotropic glutamate receptor 5 and calcium-sensing receptor). For hyperbolic agonist concentration-response relationships, no differences in estimates of ligand affinity, efficacy, or cooperativity were observed when the slope was assigned to either a transducer slope or agonist binding slope. In contrast, for receptors with nonhyperbolic agonist concentration-response relationships, estimates of ligand affinity, efficacy, or cooperativity varied depending on the assignment of the slope. The extent of this variation depended on the magnitude of the slope value and agonist efficacy, and for allosteric modulators on the magnitude of cooperativity. The modified operational models described herein are well suited to analyzing agonist and modulator interactions at receptors that bind multiple orthosteric agonists in a cooperative manner. Accounting for cooperative agonist binding is essential to accurately quantify agonist and drug actions.

SIGNIFICANCE STATEMENT Some orthosteric agonists bind to multiple sites on a receptor, but current analytical methods to characterize such interactions are limited. Herein, we develop and validate operational models of agonism and allosterism for receptors with multiple orthosteric binding sites, and demonstrate that such models are essential to accurately quantify agonist and drug actions. These findings have important implications for the discovery and development of drugs targeting receptors such as the calcium-sensing receptor, which binds at least five calcium ions.

Introduction

The past 30 years have seen major advances in quantifying the relationship between receptor occupancy and response, with the operational model of agonism (Black and Leff, 1983) representing one of the most common analytical approaches. The operational model of agonism describes agonist effects based on agonist affinity [equilibrium dissociation constant of the orthosteric agonist (KA)] and observed efficacy in a given test system. The latter is defined by a transducer ratio (τ), which is a function of both tissue- and agonist-specific components; it is the ratio of the total receptor number and a transducer parameter that defines the avidity with which a given agonist-occupied receptor complex promotes the final observed pharmacological effect. As such, the operational model of agonism is a useful tool for quantifying agonism in a comparable manner across different test systems (Black and Leff, 1983), and has subsequently been extended or modified to also quantify effects of allosteric modulators and biased agonists (Leach et al., 2007, 2010; Kenakin, 2012).

The operational model of agonism has been most commonly applied to characterize the activity of agonists that display both rectangular hyperbolic or nonhyperbolic concentration-response curves, i.e., normally empirically characterized by Hill slopes that are equal to or different from unity, respectively. The key underlying assumption in the majority of instances to date where an agonist concentration-response curve displays a Hill slope significantly different from 1has been ascribed in the most common form of the operational model to differences in the postreceptor machinery that transduces occupancy to response, i.e., through introduction of a so-called transducer slope (n) (Black et al., 1985). For instance, steep or shallow Hill slopes could arise due to changes in the sensitivity of one or more steps in a receptor’s signal transduction mechanism, while the initial agonist-receptor binding event is assumed to be characterized by a simple hyperbolic one-to-one relationship. However, for ion channels and a number of G protein–coupled receptors (GPCRs), particularly the class C GPCR subfamily, nonhyperbolic concentration-response relationships can also arise from cooperative binding of multiple equivalents of the same endogenous agonist molecule prior to any subsequent processing of the stimulus by the cellular transduction machinery. For example, while a number of small molecule calcium-sensing receptor (CaSR) agonists produce responses characterized by Hill slopes close to unity (Cook et al., 2015; Keller et al., 2018), indicating a transducer slope equal to unity, it is also well established that CaSR responses to its endogenous activator, extracellular calcium (Embedded Image), and other divalent cations are characterized by extremely high Hill slopes, ranging from 2 to 4 (Brown, 1983; Davey et al., 2012; Leach et al., 2015). The most parsimonious explanation to account for these disparate observations is that the operational transducer slope linking CaSR agonist occupancy to response can adequately be described by a transducer slope equal to unity, which suggests that the cooperativity observed in response to activators such as Embedded Image ions arises at the level of binding, not function. This is also in accord with known pharmacological and structural studies of the CaSR that have identified multiple binding sites for Embedded Image ions (Geng et al., 2016). As a consequence, the classic operational model of agonism as applied to concentration-response curves of nonunit Hill slopes is suboptimal for such situations.

Herein, we sought to develop and evaluate an operational model of agonism that describes orthosteric agonist binding to multiple sites in a cooperative manner, referred to as the cooperative agonist operational model. The cooperative agonist operational model was superior to the original Black-Leff operational model of agonism in fitting Embedded Image-CaSR concentration-response data. We also extended this cooperative agonist operational model to incorporate allosteric modulation of the affinity and efficacy of an agonist that binds cooperatively to multiple sites. This operational model of allosterism with cooperative agonist binding was fitted to data describing the actions of CaSR positive allosteric modulators (PAMs) and negative allosteric modulators (NAMs), and revealed that if cooperative agonist binding is not taken into consideration, under- or overestimates of PAM and NAM affinity and cooperativity can occur.

Materials and Methods

Materials.

Dulbecco’s modified Eagle’s medium (DMEM), Flp-In human embryonic kidney (HEK) T-REx cells, blasticidin S HCl, and FBS were obtained from Invitrogen (Carlsbad, CA), while hygromycin B was obtained from Roche (Mannheim, Germany). Fluo-8 AM was obtained from Abcam (Cambridge, MA).

CaSR-Expressing HEK293 Cell Lines.

The generation of DNA and Flp-In HEK T-REx cells stably expressing c-myc-tagged wild-type (WT) CaSR in pcDNA5/frt/TO have been described previously (Davey et al., 2012; Leach et al., 2016). Flp-In HEK T-REx CaSR cells were maintained in DMEM containing 5% FBS, 200 μg/ml hygromycin B, and 5 μg/ml blasticidin S HCl. To generate a tetracycline-inducible Flp-In HEK cell line stably expressing an N-terminally truncated CaSR, N-terminally truncated CaSR corresponding to amino acids 600–903 with an N-terminal influenza hemagglutinin signal peptide followed by a c-myc epitope and rhodopsin signal peptide in pcDNA3.1+ (Leach et al., 2016) was transferred to pcDNA5/frt/TO using BamHI and NotI restriction sites. Flp-In HEK T-REx cells were seeded into 25 cm2 flasks in DMEM containing 5% FBS and allowed to reach 80% confluency. Cells were transfected with 0.5 µg pcDNA5/frt/TO containing the N-terminally truncated CaSR plus 5 µg POG44 with lipofectamine 2000 (1:4 DNA:lipofectamine 2000) according to the manufacturer’s instructions. The following day, cells were transferred to a T75 cm2 flask, and 24 hours later DMEM was replaced with DMEM containing 5% FBS, 200 µg/ml hygromycin, and 5 µg/ml blasticidin S HCl. The selection DMEM was replaced every 3 days until untransfected cells had died (∼10 days), and antibiotic-resistant cells were expanded and maintained in DMEM containing 5% FBS, 200 µg/ml hygromycin, and 5 µg/ml blasticidin S HCl. All cell lines were routinely tested for mycoplasma contamination using the Lonza MycoAlert mycoplasma detection kit.

Determination of WT and N-Terminally Truncated CaSR Cell Surface Expression Using Fluorescence-Activated Cytometry.

Flp-In HEK T-REx WT and N-terminally truncated CaSR-expressing cells were seeded at 80,000 cells per well into a 96-well plate and expression was induced with 100 ng/ml tetracycline overnight at 37°C. The following day, cells were harvested and washed in 1X PBS with 0.1% bovine serum albumin and 2 mM EDTA (wash buffer) by centrifugation (350g, 4°C for 3 minutes) before resuspension and 30-minute incubation in 100 μl blocking buffer (1X PBS, 5% bovine serum albumin, and 2 mM EDTA) containing 1 μg/ml AF647-conjugated 9E10 made in-house as previously described (Cook et al., 2015). Cells were washed as previously described and resuspended in wash buffer containing propidium iodide. Live cell fluorescence was measured using a BD FACS Canto analyzer (Becton Dickinson).

Calcium Mobilization Assays.

Cells were seeded in clear 96-well plates coated with poly-D-lysine (50 μg/ml−1) at 80,000 cells per well and incubated overnight in the presence of 0 or 100 ng/ml−1 tetracycline. The following day, cells were washed with assay buffer (150 mM NaCl, 2.6 mM KCl, 1.18 mM MgCl2, 10 mM D-Glucose, 10 mM HEPES, 0.1 mM CaCl2, 0.5% bovine serum albumin, and 4 mM probenecid at pH 7.4) and loaded with Fluo-8 AM (1 μM in assay buffer) for 1 hour at 37°C. Cells were washed with assay buffer prior to the addition of fresh assay buffer.

For all studies, each well was treated with a single agonist concentration. The release of intracellular calcium (Ca2+i) was measured at 37°C using FlexStation 1 or 3 (Molecular Devices, Sunnyvale, CA). Fluorescence was detected for 60 seconds at 490 nm excitation and 520 nm emission and the peak Ca2+i mobilization response (approximately 12 seconds after agonist addition) was used for the subsequent determination of the agonist response. Relative peak fluorescence units were normalized to the fluorescence stimulated by 1 μM ionomycin to account for differences in cell number and loading efficiency.

Model Derivation and Data Analysis.

The derivation of operational models describing the effect of an agonist in the absence or presence of an allosteric modulator at a receptor with multiple agonist binding sites is presented in the Supplemental Appendix. The script input to use the two equations in the GraphPad Prism program is also presented in the Supplemental Appendix.

Data simulations were performed using the original Black-Leff operational model of agonism, referred to hereinafter as the Black-Leff model (eq. 1), or a modified cooperative agonist operational model (eq. 2), where an additional slope [i.e., the binding slope linking agonist concentration to receptor occupancy (nB)] was incorporated to take into account multiple agonist binding sites, and thus the steepness of the slope describing the agonist concentration-occupancy relationship:Embedded Image(1)Embedded Image(2)where [A] is the agonist concentration, KA is the agonist equilibrium dissociation constant; τA is an operational measure of agonist efficacy; Em is the maximal system response; and nT is the transducer slope linking agonist concentration to response.

For simplicity, the modified cooperative operational model of agonism (eq. 2) makes the following assumptions:

  1. The receptor is either empty or fully occupied, and only the fully occupied receptor exerts an effect. The lack of partially occupied receptor molecules could arise if multiple agonist molecules bind simultaneously to the receptor, or if agonist molecules bind sequentially with high positive cooperativity. Thus, once one binding site is occupied, positive cooperativity drives occupancy of all other sites. The latter scenario is most likely; hence, this is why we have called this model the cooperative agonist operational model.

  2. The model cannot discern cooperativity between the multiple binding sites, thus nB may not be the true number of binding sites; therefore, nB is a binding slope coefficient.

  3. The KA value is a geometric mean of the microscopic dissociation constants for each binding site.

The Eqs. 1 or 2 were fitted to agonist concentration response data in order to quantify agonist affinity and efficacy. When fitting eq. 2 to experimental data, the transducer slope was constrained to unity (see Results for validation of this assumption).

Data describing the interaction between glutamate and PAMs and NAMs at metabotropic glutamate receptor (mGlu) subtype 5, or between acetylcholine (ACh) and 3-amino-5-chloro-6-methoxy-4-methyl-thieno[2,3-b]pyridine-2-carboxylic acid cyclopropylamide (LY2033298) at the M4 muscarinic acetylcholine receptor (mAChR), were fitted to our original operational model of allosterism (eq. 3) or to the new operational model of allosterism with cooperative agonist binding (eq. 4). For simplicity and for the purpose of fitting experimental data, eqs. 3–8 assume a single allosteric modulator binding site; therefore, they do not take into account modulator cooperative binding:Embedded Image(3)Embedded Image(4)where KA is the equilibrium dissociation constant of the orthosteric agonist, which was fixed in some instances to the affinity determined in radioligand binding assays (Mutel et al., 2000; Leach et al., 2010); KB is the equilibrium dissociation constant of the allosteric ligand; τA and τB are the operational efficacies of the orthosteric agonist and allosteric ligand, respectively; α and β are the allosteric effects on orthosteric agonist affinity and efficacy, respectively (it should be noted that β is not a reciprocal efficacy cooperativity factor) (Leach et al., 2007; Giraldo, 2015); and [A] and [B] are the orthosteric agonist and allosteric ligand concentrations, respectively.

To fit the operational model of allosterism to data describing the interaction between Embedded Image and cinacalcet at the CaSR, the original operational model of allosterism shown by eq. 3 was simplified, because for a full agonist like Embedded Image (i.e., one that generates the maximal system response at submaximal receptor occupancies) KA >> [A]. Furthermore, because the CaSR’s orthosteric agonist, Embedded Image, was present in the assay buffer, the contaminating agonist was included in the equations used to analyze CaSR PAM (cinacalcet) and NAM (2-chloro-6-[(2R)-2-hydroxy-3-[(2-methyl-1-naphthalen-2-ylpropan-2-yl)amino]propoxy]benzonitrile [NPS2143]) data (Keller et al., 2018). Therefore, data describing the interaction between Embedded Image and cinacalcet or NPS2143 at the CaSR were fitted to the original operational model of allosterism with the contaminating agonist (eqs. 5 and 6, respectively) or to an operational model of allosterism with cooperative agonist binding and contaminating agonist (eqs. 7 and 8, respectively):Embedded Image(5)where EC50 is the agonist concentration that elicits a half-maximal response, in which it should be noted that inclusion of [EC50] involves some simplifying assumptions that facilitate data fitting (Aurelio et al., 2009; Giraldo, 2015); [C] is the contaminating agonist concentration; and all other parameters are as described for eq. 3:Embedded Image(6)Embedded Image(7)where EC50, KB, τB, α, β, [A], [B], [C], and Em are as described for eq. 5:Embedded Image(8)where all other parameters are as described for eq. 4.

The Hill equation (Eq 9) was fitted to simulated data:Embedded Image(9)where [A] is the agonist concentration; Emax is the maximum agonist effect; and nH is the Hill slope.

Nonlinear regression analysis was performed in GraphPad Prism 7 or 8. Potency, affinity, cooperativity, and efficacy parameters were estimated as logarithms (Christopoulos, 1998). An extra sum of squares F test was used to determine whether data were fitted best when the Hill slope, binding slope, or transducer slope (nH, nB, or nT, respectively) were significantly different from unity, where P < 0.05 was considered significant.

Results

The Contribution of Slope Factors to Agonist Concentration-Response Relationships.

We first evaluated the contribution of the agonist binding slope (nB) or transducer slope (nT) to the concentration-response curve of two agonists with different efficacies by simulating variations in nB or nT using the cooperative agonist operational model (eq. 2). We specifically wanted to evaluate a system with cooperative agonist binding; therefore, we based our simulations on Embedded Image activation of the CaSR. The Embedded Image concentration-response relationship for the CaSR’s best characterized physiologic role, inhibition of parathyroid hormone secretion, occurs over a Embedded Image concentration range of 0.8–1.5 mM with an EC50 value of ∼1.2 mM (Brown, 1983). Thus, for these simulations, the affinity of the agonist (KA) was assumed to be 1.2 mM and nB or nT were assumed to be between 1 and 3. Simulated data were subsequently fitted to a Hill equation (eq. 9). Unsurprisingly, increasing nB or nT increases the Hill slope of the agonist concentration-response curve (Fig. 1; Table 1). Furthermore, increasing nT decreases agonist potency. Interestingly, the effect of nB on agonist potency depends on the magnitude of nT and τA. For instance, increasing nB decreases agonist potency for higher efficacy agonists (τA 3). However, for lower efficacy agonists (τA 1), increasing nB decreases agonist potency when nT ≤ 1, but increases agonist potency when nT ≥ 2.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

The nT and nB values differentially contribute to agonist concentration-response relationships. Simulations demonstrating the influence of nT or nB on concentration-response relationships for agonists with different efficacies (τA). Data were simulated using the cooperative agonist operational model (eq. 2), where the affinity of the agonist (KA) is 1.2 mM and nT and nB are between 1 and 3. Curves through the data are the fits to a four-parameter Hill equation (eq. 9), where the parameters describing the fits are given in Table 1.

View this table:
  • View inline
  • View popup
TABLE 1

Simulation of agonist concentration-response relationships upon changes in binding or transducer slopes and τA

Data were simulated using the cooperative agonist operational model (eq. 2) and a Hill equation (eq. 9) was fitted to simulated data to determine agonist potency, maximum effect, and Hill slope.

We next sought to directly compare the influence of the binding or transducer slope by simulating concentration-response curves for agonists with varying efficacies using the Black-Leff operational model (eq. 1, which contains a transducer slope, nT) and the cooperative agonist operational model (eq. 2, which contains transducer and binding slopes, nT and nB, respectively). As can be seen in Fig. 2 and Table 2, when nT or nB are equal to 1, variations in τA have an identical effect on empirical agonist concentration-response parameters (potency, Hill slope, or Emax) regardless of the model. In contrast, when nT or nB are greater than 1, variations in τA result in major differences in the agonist concentration-response profile predicted with the two different operational models of agonism. Specifically, the Black-Leff model predicts that high-efficacy agonists have greater potency relative to affinity (due to amplification of the steps between agonist binding and response), while for low-efficacy agonists, the EC50 value may be less than the KA value for curves that possess nonunity Hill slopes. The latter was previously noted by Black et al. (1985). Furthermore, the Hill slope decreases alongside decreases in τA. In comparison, the cooperative agonist operational model predicts that when nT = 1, agonist EC50 may approach but not be less than its KA value, regardless of whether nB > 1, and there is no effect of τA on the Hill slope of the agonist concentration-response curve (Fig. 2; Table 2).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Cooperative agonist binding influences agonist concentration-response relationships. Simulations demonstrating the influence of agonist efficacy (τA) on agonist concentration-response relationships when the slope is governed by the transducer slope (nT, Black-Leff operational model of agonism) or the agonist binding slope (nB, cooperative agonist operational model). Data were simulated using the Black-Leff model (eq. 1) or the cooperative agonist operational model (eq. 2), where the affinity of the agonist (KA) is 1.2 mM and nT and nB are between 1 and 3. Curves through the data are the fits to a Hill equation (eq. 9), where the parameters describing the fits are given in Table 2.

View this table:
  • View inline
  • View popup
TABLE 2

Simulation of agonist concentration-response relationships upon changes in binding or transducer slopes and τA

Data were simulated using the Black-Leff model (eq. 1) or the cooperative agonist operational model (eq. 2) and a Hill equation (eq. 9) was fitted to simulated data to determine agonist potency, maximum effect, and Hill slope.

Quantification of Experimentally Derived Agonist Concentration-Response Data.

We next tested whether our simulations were recapitulated in a functional assay that measures CaSR activation. To do so, we measured Embedded Image-mediated Ca2+i mobilization following titration of CaSR expression using a tetracycline inducible system. In the absence of tetracycline, the maximal response to Embedded Image is approximately 50% of the maximal response obtained under full induction of receptor expression (100 ng/ml tetracycline). In this system, fitting a Hill equation (eq. 9) to both data sets indicated that the data were fitted best when the Hill slope was unchanged with different receptor expression levels, i.e., different magnitudes of τA (P < 0.05, extra sum-of-squares F test; data not shown). For the CaSR, the small molecule allosteric agonists 1-(1,3-benzothiazol-2-yl)-1-(2,4-dimethylphenyl)ethanol (AC265347) and R-1-(6-methoxy-4′-(trifluoromethyl)-3-biphenylyl)-N-(R)- 1-phenylethyl)ethanamine (calcimimetic B) activate the CaSR with a Hill slope of 1 (Cook et al., 2015; Keller et al., 2018). Similarly, when cooperative agonist binding is prevented by removal of the CaSR’s N-terminal domain, and consequently the primary Embedded Image binding sites, the Hill slope for Embedded Image is not significantly different from unity (as shown subsequently). This provides experimental evidence that the CaSR’s transducer slope is equal to unity, and that the steep Hill slopes observed for Embedded Image at the full-length CaSR thus likely arise from a binding slope greater than 1. Thus, when fitting CaSR experimental data to the cooperative agonist operational model, nT was constrained to unity.

When the data were fitted to the classic Black-Leff model, the estimated KA value was 0.2 mM (Fig. 3A; Table 3). In comparison, the cooperative agonist operational model yielded a KA estimate of 1.1 mM, which is in close agreement with the EC50 value (1.2 mM) of Embedded Image for suppressing parathyroid hormone release (Brown, 1983) and Embedded Image affinity estimates for the CaSR extracellular domain determined using spectroscopic approaches (Huang et al., 2009; Zhang et al., 2014). For both analyses, data were fitted best when the binding slope (cooperative agonist operational model) or transducer slope (Black-Leff model) were different from unity (P < 0.05, extra sum of squares F test).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Embedded Image-CaSR concentration-response relationships fit well to the cooperative agonist operational model. (A) Embedded Image-mediated Ca2+i mobilization at the WT CaSR following overnight receptor induction with 100 ng/ml tetracycline (tet) or in the absence of tetracycline. Data are mean + S.D. from four independent experiments performed in duplicate. Curves through the data are the fits to the Black-Leff model (blue line) or the cooperative agonist operational model (red line), where the parameters describing the fits are given in Table 3. Although both models fit the data, the cooperative agonist operational model more accurately predicts the expected affinity of Embedded Image at the CaSR (Table 3). (B) Extracellular trivalent gadolinium (Embedded Image) and Embedded Image concentration-response curves at the WT CaSR following overnight receptor induction with 100 ng/ml tetracycline. Data are mean + S.D. from four independent experiments performed in duplicate. Curves through the data are the fits to a four-parameter Hill equation. (C) N-terminally truncated CaSR snake diagram. (D) Embedded Image and Embedded Image concentration-response curves at an N-terminally truncated CaSR following overnight receptor induction with 100 ng/ml tetracycline. Data are mean + S.D. from five independent experiments performed in duplicate. Curves through the data are the fits to the cooperative agonist operational model, where the parameters describing the fits are given in Table 3.

View this table:
  • View inline
  • View popup
TABLE 3

Quantification of Embedded Image affinity and efficacy for CaSR-Ca2+i mobilization using the operational model of agonism

Embedded Image concentration-response curves were generated at the WT receptor following overnight incubation of cells with or without 100 ng/ml tetracycline to induce receptor expression. For N-terminally truncated CaSRs, Embedded Image and Embedded Image concentration-response curves were generated upon induction of expression with 100 ng/ml tetracycline, where Embedded Image is a partial agonist and Embedded Image is a full agonist. The Black-Leff model (eq. 1) or the cooperative agonist operational model (eq. 2) was fitted to data to determine agonist affinity, efficacy, and transducer or binding slope.

To further validate our simulations in a functional assay, we next sought to quantify the affinity and efficacy of a CaSR partial agonist. To do so, we took advantage of observations that in comparison with at the WT CaSR (Fig. 3B), Embedded Image acting via the 7 transmembrane domain is a partial agonist at an N-terminally truncated CaSR (depicted in Fig. 3C) relative to the extracellular trivalent gadolinium cation (Fig. 3D). The fluorescence-activated cytometry analysis confirmed cell surface expression of the WT and N-terminally truncated CaSR (Supplemental Fig. 1). We quantified Embedded Image affinity and efficacy at the N-terminally truncated CaSR using the original Black-Leff model or the cooperative agonist operational model (Table 3). In both instances, data were best fitted when the binding slope (the cooperative agonist operational model) or transducer slope (the Black-Leff model) were not different from unity (P < 0.05, extra sum of squares F test), which is consistent with a reduction in positively cooperative Embedded Image binding, as would be expected upon removal of the four primary Embedded Image binding sites located in the N-terminal domain (Geng et al., 2016). Thus, the parameters determined at the N-terminally truncated receptor were identical regardless of the equation used to analyze the data. The Embedded Image affinity and cooperativity estimates at the N-terminally truncated receptor were compared with those determined at the full-length WT CaSR, indicating a reduction in Embedded Image affinity at the N-terminally truncated receptor in comparison with WT (Table 3), again consistent with removal of the primary Embedded Image binding sites. However, due to a lower estimate of Embedded Image efficacy at the WT receptor when WT data were analyzed using the Black-Leff model, only the cooperative agonist operational model accurately quantified a reduction in Embedded Image efficacy at the N-terminally truncated receptor in comparison with at the WT receptor. This is consistent with a lower Embedded Image Emax value at the N-terminally truncated receptor (∼60% of the maximum response stimulated by extracellular trivalent gadolinium) in comparison with WT (∼100% extracellular trivalent gadolinium Emax) (Table 3). Thus, only the cooperative agonist operational model accurately estimated Embedded Image partial agonism at the N-terminally truncated receptor.

Quantifying Allosteric Interactions in Systems with Different Degrees of Cooperative Agonist Binding.

Having established that the cooperative agonist operational model best fitted our Embedded Image-WT CaSR concentration-response curves with Hill slopes greater than 1, we next extended this model to allow for quantification of allosteric modulation of an agonist response. The operational model of agonism and allosterism (Leach et al., 2007, 2010) (referred to herein as the original operational model of agonism and allosterism), which takes into account the allosteric effects on agonist affinity and efficacy, combines the allosteric ternary complex models developed by Stockton et al. (1983) and Ehlert (1988) and the Black-Leff operational model of agonism. In our original model (Leach et al., 2007), the allosteric modulator can also possess intrinsic efficacy. Introduction of a slope in that model once again assumed that the slope linked occupancy to response, not to the original binding events, which were assumed to be described as simple one-to-one hyperbolic functions. Therefore, we adapted this operational model of allosterism to account for cooperative agonist binding, referred to hereinafter as the operational model of allosterism with cooperative agonist binding. To validate this operational model of allosterism with cooperative agonist binding, we reanalyzed existing data demonstrating positive and negative allosteric modulation at three model GPCRs with different agonist Hill slopes: CaSR (a class C GPCR where the primary endogenous agonist, Embedded Image, has a Hill slope of 2–4), mGlu5 (a class C GPCR where the primary endogenous agonist, L-glutamate, has a Hill slope of ∼1.8) (Sengmany and Gregory, 2016), and M4 mAChR (a class A GPCR where the endogenous agonist, acetylcholine, has a Hill slope of 1) (Leach et al., 2010, 2011) (Supplemental Fig. 2). In all instances, nT was assumed to be unity and all allosteric modulators were assumed to bind to a single site (i.e., the modulator binding slope is unity).

For the CaSR, we analyzed allosteric modulation of Embedded Image by cinacalcet (PAM) or NPS2143 (NAM) (Leach et al., 2016) with the original operational model of agonism and allosterism with contaminating (i.e., ambient buffer) agonist (eqs. 5 or 6, respectively) and the newly derived operational model of allosterism with cooperative agonist binding and contaminating agonist (eqs. 7 or 8, respectively) (Fig. 4). Similar to our analysis of agonist concentration-response curves, data were fitted best when the binding slope (operational model of allosterism with cooperative agonist binding) or transducer slope (original operational model of agonism and allosterism) was different from unity (P < 0.05, extra sum of squares F test). Compared with the original operational model of agonism and allosterism, the estimated affinity for Embedded Image determined using the operational model of allosterism with cooperative agonist binding (1.4 mM) (Table 4) was once again closer to the assumed Embedded Image affinity based on its EC50 value for suppression of parathyroid hormone release (1.2 mM) and quantification of the Embedded Image KA at the extracellular domain using spectroscopic approaches (3–5 mM) (Huang et al., 2009; Zhang et al., 2014). Furthermore, the estimated affinity and negative cooperativity of NPS2143 were greater (5.5- and 35-fold, respectively) when cooperative agonist binding was factored into the analysis (Table 4). For the PAM, cinacalcet, the operational model of allosterism with cooperative agonist binding yielded a 3-fold lower affinity estimate but an 8-fold greater magnitude of positive cooperativity.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Allosteric modulation at the CaSR is fitted well by an operational model of allosterism with cooperative agonist binding. Allosteric modulation of Embedded Image-mediated Ca2+i mobilization at the CaSR by NPS2143 (NAM) or cinacalcet (PAM). Data were previously published (Leach et al., 2016) and are mean + S.D. from at least 11 independent experiments performed in duplicate. Curves through the data are the fits to the operational model of allosterism with cooperative agonist binding and contaminating agonist (eq. 7 for cinacalcet and eq. 8 for NPS2143), where the parameters describing the fits are given in Table 4.

View this table:
  • View inline
  • View popup
TABLE 4

Comparison of parameters describing CaSR allosteric interactions analyzed with different allosteric models

CaSR allosteric interactions with Embedded Image in Ca2+i mobilization assays were analyzed with the original operational model of agonism and allosterism (eqs. 5 or 6) vs. the operational model of allosterism with cooperative agonist binding (eqs. 7 or 8) to determine Embedded Image potency, efficacy, and affinity; modulator affinity, cooperativity, and transducer or binding slope; and maximum system response.

We next analyzed allosteric modulation of glutamate at mGlu5 (eqs. 3 or 4) by a representative full NAM (2-methyl-6-(2-phenylethynyl)pyridine [MPEP]) that completely inhibits glutamate-mediated activation of Ca2+i mobilization, a partial NAM (2-[2-(3-methoxyphenyl)ethynyl]-5-methylpyridine [M-5MPEP]) that only partially inhibits glutamate-mediated activation of Ca2+i mobilization (Sengmany et al., 2019), a pure PAM (N-(1,3-Diphenyl-1H-pyrazolo-5-yl)-4-nitrobenzamide [VU-29]), and a mixed PAM-agonist (1-(4-(2,4-difluorophenyl)piperazin-1-yl)-2-((4-fluorobenzyl)oxy)ethanone [DPFE]) (Sengmany et al., 2017). Similar to analyses at the CaSR, all data were fitted best when the binding slope (operational model of allosterism with cooperative agonist binding) or transducer slope (original operational model of agonism and allosterism) was different from unity (P < 0.05, extra sum of squares F test). However, for each modulator, the affinity and cooperativity estimates were similar (within 3-fold) irrespective of the analytical model applied (Fig. 5; Table 5). Therefore, although the glutamate-mGlu5 concentration-response relationship has a Hill slope greater than unity, quantification of allosteric interactions at mGlu5 is largely unaffected by whether the empirical slope is assumed to be determined by the transducer slope or the agonist binding slope.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Allosteric modulation at mGlu5 is fitted well by an operational model of allosterism with cooperative agonist binding. Allosteric modulation of glutamate-mediated Ca2+i mobilization at mGlu5 by MPEP (full NAM), M-5MPEP (partial NAM), VU-29 (PAM), and DPFE (ago-PAM). Data were previously published (Sengmany et al., 2017, 2019) and are mean + S.D. from at least three independent experiments performed in duplicate previously published. Curves through the data are the fits to the operational model of allosterism with cooperative agonist binding (eq. 4), where the parameters describing the fits are given in Table 5.

View this table:
  • View inline
  • View popup
TABLE 5

Comparison of parameters describing mGlu5 allosteric interactions analyzed with different allosteric models

mGlu5 allosteric interactions with glutamate in Ca2+i mobilization assays were analyzed with the original operational model of agonism and allosterism (eq. 3) or the operational model of allosterism with cooperative agonist binding (eq. 4) to determine modulator affinity, glutamate or modulator efficacy, respectively), cooperativity, transducer or binding slopes, and maximum system response.

For the M4 mAChR, we analyzed previously published positive allosteric modulation of ACh by the PAM agonist LY2033298 in guanosine 5′-O-(3-[35S]thio)triphosphate binding assays (Leach et al., 2010) (eqs. 3 or 4). As expected, in the absence of cooperative ACh binding at the M4 mAChR, data for the interaction between ACh and LY2033298 were fitted best by both operational models of agonism and allosterism when the slope was not different from unity (P < 0.05, extra sum of squares F test); therefore, both equations yielded identical estimates of affinity and cooperativity (Fig. 6; Table 6).

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Allosteric modulation at the M4 mAChR is fitted well by an operational model of allosterism with no cooperative agonist binding. Positive modulation of ACh-mediated guanosine 5′-O-(3-[35S]thio)triphosphate ([35S]GTPγS) binding by LY2033298 (ago-PAM) at the M4 mAChR. Data were previously published (Leach et al., 2010) and are mean + S.D. from at least three independent experiments performed in duplicate. Curves through the data are the fits to the operational model of allosterism with cooperative agonist binding (eq. 4), where the equation was fitted best when n = 1 (i.e., no cooperative agonist binding). The parameters describing the fits are given in Table 6.

View this table:
  • View inline
  • View popup
TABLE 6

Comparison of parameters describing M4 mAChR allosteric interactions analyzed with different allosteric models

M4 mAChR LY2033298 interactions with ACh in guanosine 5′-O-(3-[35S]thio)triphosphate binding assays were analyzed with the original operational model of agonism and allosterism (eq. 3) or the operational model of allosterism with cooperative agonist binding (eq. 4) to determine agonist or modulator affinity, agonist or modulator efficacy, cooperativity, transducer or binding slopes, and maximum system response. Values obtained from fitting eqs. 3 and 4 to the data were identical; therefore, they are only presented in a single column.

We next sought to establish why quantification of PAM and NAM affinity and cooperativity were not greatly affected by the assignment of the slope at mGlu5, where the glutamate Hill slope is greater than unity. To do so, we simulated the interaction between an orthosteric agonist and a NAM or PAM with the operational model of allosterism with cooperative agonist binding and analyzed the simulated data with the original operational model of agonism and allosterism. For these simulations, orthosteric agonist affinity (1 µM), τA (10), and modulator affinity (10 nM) were held constant, and different magnitudes of positive or negative cooperativity were examined alongside changes in the magnitude of cooperative agonist binding. Consistent with our analysis of mGlu5 allosteric interaction data, when the agonist binding slope ranged from 1 to 2 the affinity of NAMs and PAMs estimated using the original operational model of agonism and allosterism fell within ∼3-fold of the affinity simulated with the operational model of allosterism with cooperative agonist binding. Furthermore, only the affinity of the full NAM (where αβ is assumed to approach zero) fell outside this 3-fold window once the agonist binding slope exceeded 2 (Fig. 7A). Similarly, when analyzed with the original operational model of agonism and allosterism, the αβ values (0.1–10) were within a 3-fold range of simulated data when nT ranged from 1 to 2. However, where the magnitude of cooperativity exceeded 10 for a PAM, or 0.1 for a NAM, the influence of cooperative agonist binding becomes more pronounced (Fig. 7B). These simulations confirm that the assignment of n to a transducer or binding slope would not be expected to greatly influence quantification of mGlu5 allosteric modulation of glutamate for modulators with moderate cooperativity, which is consistent with our Ca2+i mobilization data.

Fig. 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 7.

Modulator affinity and cooperativity are influenced by assignment of the slope in the operational model of allosterism. Simulations demonstrating the influence of the agonist binding slope on the estimated modulator affinity and cooperativity when interaction data are analyzed using the original operational model of allosterism. Interaction data between an orthosteric agonist and a NAM or PAM were simulated with the operational model of allosterism with cooperative agonist binding. Orthosteric agonist affinity (1 µM), τA (10), and modulator affinity (10 nM) were held constant; different magnitudes of positive or negative cooperativity were examined alongside changes in the magnitude of cooperative agonist binding. The simulated data were analyzed with the original operational model of allosterism, and log KB (A) or log αβ (B) estimates were plotted against the agonist binding slope.

Discussion

In the current study, we have assessed operational models of agonism and allosterism that account for receptors whose agonists bind multiple binding sites in a cooperative manner. The modified models accurately fit experimental data at an exemplar GPCR, the CaSR, which has high sensitivity for Embedded Image due to multiple Embedded Image binding sites that are linked in a positively cooperative manner. We show that agonist Hill slopes that differ from unity and remain unchanged by alterations in receptor expression levels or cellular coupling efficiencies (i.e., where τA differs) may be indicative of cooperative agonist binding. We demonstrate that if a steep Hill slope such as that observed at the CaSR is attributed to the transducer slope rather than to the agonist binding slope, the Black-Leff operational model of agonism underestimates agonist efficacy and overestimates agonist affinity. Extension to allosteric interactions shows the importance of accounting for cooperative agonist binding, since different models fitted to the same allosteric interaction data yield divergent modulator affinity and cooperativity estimates. For instance, the original operational model of agonism and allosterism estimates lower CaSR PAM and NAM cooperativity values and higher or lower affinity values, respectively. Data simulations support these findings and demonstrate that the impact of cooperative binding on estimates of modulator affinity and cooperativity is more pronounced as the magnitude of modulator cooperativity or cooperative agonist binding is increased. This was evidenced by our demonstration that for mGlu5, where the glutamate Hill slope is ∼1.8, differences in affinity and cooperativity are within the margin of experimental error (∼3-fold range). Accordingly, for agonist-receptor concentration-response relationships with Hill slopes equal to unity, it does not matter whether the slope is governed by the transducer slope or the agonist binding slope.

The operational models of agonism and allosterism with cooperative agonist binding have important practical uses for analyzing data at receptors that possess multiple agonist binding sites. This is particularly true for the CaSR, but also for ion channels and other GPCRs where agonist binding coefficients differ from unity, such as GPR39, which binds at least two Zn2+ ions and responds to Zn2+ with a Hill slope of 2–3 (Storjohann et al., 2008; Sato et al., 2016). Similarly, cooperative binding can occur across a GPCR dimer, which may account for the steep Hill slope at mGlu5 demonstrated in the present study. For instance, the mGlu2 orthosteric agonist (1S,2S,5R,6S)-2-aminobicyclo[3.1.0]hexane-2,6-dicarboxylic acid (LY354740) stabilizes conformational rearrangements of a metabotropic glutamate receptor subtype 2 and 4 heterodimer (mGlu2-mGlu4) with a shallow Hill slope (Moreno Delgado et al., 2017), which is increased to unity when LY354740 is prevented from binding to the mGlu2 or mGlu4 orthosteric binding site in the dimer. In contrast, the glutamate Emax value is reduced when it can bind to only one of the orthosteric sites in the heterodimer (Moreno Delgado et al., 2017). These findings indicate negative (LY354740) and positive (glutamate) cooperativity across the dimer, respectively. Negative cooperative binding has also been demonstrated at several class A GPCRs dimers, including the 5-hydroxytryptamine2A, A3 adenosine, H3 histamine, and D2 dopamine receptors (Sinkins and Wells, 1993; Vivo et al., 2006; Brea et al., 2009; May et al., 2011). Accurate quantification of ligand affinity, efficacy, and cooperativity at such receptors using functional assays is critical, particularly where radioligand binding–based methods are not available. Indeed, there are no commercially available radioligands for the CaSR, thus pharmacological characterization, and indeed drug discovery at this receptor has generally relied upon functional measures of receptor activity to quantify drug actions. However, it must be noted that experimentally derived pharmacological data for agonists with steep Hill slopes will likely only be fitted to the cooperative agonist operational model and the operational model of allosterism with cooperative agonist binding when one of the binding or transducer slopes (nB or nT) is known and constrained as such. This was a key advantage of analyzing pharmacological data at the CaSR, where we showed experimentally that the transducer slope is 1.

Our findings have important implications for past and present drug discovery efforts at class C GPCRs and beyond. Establishing structure-activity relationship profiles that dictate drug affinity and cooperativity is essential for predicting drug efficacy in vivo. However, underestimates of cooperativity at class C GPCRs with cooperative agonist binding may explain previous observations that class C GPCR allosteric modulators have limited cooperativity when compared with their class A GPCR counterparts. For example, PAMs with αβ values >100 have been reported for class A GPCRs (Leach et al., 2010; Abdul-Ridha et al., 2014), whereas for the CaSR potentiation is at most 5-fold for many modulators (Cook et al., 2015; Leach et al., 2016; Diepenhorst et al., 2018). Thus, for GPCRs with cooperative agonist binding, larger differences between modulator cooperativities were likely previously unappreciated. This is important because allosteric modulator cooperativity can predict likely clinical efficacy or adverse effect liability. Inaccurate estimation of allosteric modulator affinity or cooperativity due to a failure to consider cooperative agonist binding likely also impacts interpretation of structure-function studies. If cooperativity values are narrowed, then more subtle effects of mutations on modulator cooperativity may have been missed.

In conclusion, we have validated a method for quantifying agonist and allosteric modulator actions at receptors possessing multiple agonist binding sites that interact in a cooperative manner. Our operational models of agonism and allosterism with cooperative agonist binding more accurately quantify the actions of both orthosteric and allosteric drugs acting at GPCRs with cooperative agonist binding and may be used for future drug discovery efforts at these important receptors.

Appendix

Cooperative agonist operational model

A model of signal transduction by a receptor with nB binding sites is displayed in the equilibrium:Embedded ImageWhere KE is the value of AnBR that elicits half the maximal system effect and agonist affinity for each site can be described by the equilibrium dissociation constants:Embedded Image

We see that Embedded Image, where KA is the geometric mean of the individual equilibrium dissociation constants.

For simplicity, the receptor is considered either empty (R) or fully occupied (AnBR): Embedded ImageEmbedded Image

The total receptor concentration can be expressed as:Embedded ImagewhereEmbedded ImagethereforeEmbedded ImageEmbedded ImageEmbedded Image

Receptor occupancy is thus denoted: Embedded Image

In accordance with the scheme of the operational model of agonism, the logistic function for the transduction of occupancy into response is: Embedded Image where nT is a logistic slope factor describing the transduction of agonist binding into a response (the transducer slope).

Using the previous expression of [AnBR] gives:Embedded Image

Multiplying numerator and denominator by (KAnB + [A]nB)nT gives:Embedded Image

Dividing through by KE and redefining Embedded Image as τA gives an operational model of agonism for a receptor with nB binding sites (Eq 2 in the main text):Embedded Image

For use in Graphpad Prism or similar software, the above equation is described by the following notations, where nT or nB will likely need to be fixed to a known or theoretical value to fit real experimental data to this equation, and where a “Basal” response parameter is introduced to accommodate ligand-independent effects that deviate from zero:Embedded ImageEmbedded Image

Operational model of allosterism with cooperative agonist binding

In a ternary complex consisting of a receptor (lacking constitutive activity), an orthosteric agonist that binds nB binding sites, and an allosteric modulator, the stimulus-generating species are AnBR, BR and AnBRB. A model of ligand-receptor interactions in this ternary complex is displayed in the equilibrium:Embedded ImageWhere Embedded Image, Embedded Image, and Embedded Image

From the latter expression we obtain:Embedded Image

Assuming the pharmacological effect (E) is a function of the total stimulus (ST) arising from the sum of the stimuli generated from each individual ligand-bound receptor species, then:Embedded Image

It is further assumed that the stimulus (S) generated from each individual ligand-bound receptor species is proportional to the product of the intrinsic efficacy of the ligand (ε) and the concentration of the ligand-bound receptor species, thus:Embedded Imagewhere Embedded Imagedenotes the intrinsic efficacy of the orthosteric agonist, εB denotes the intrinsic efficacy of the allosteric modulator, and β is a coupling factor that describes the effect of the allosteric modulator on the intrinsic efficacy of the orthosteric agonist when the two ligands are bound to the same receptor.

Thus, the effect (E) of an agonist in the presence of an allosteric modulator is processed through the following logistic equationEmbedded Imagewhere Em denotes the maximum possible response, KS denotes a constant that governs the efficiency of stimulus-response coupling and nT denotes a logistic slope factor (the transducer slope).

If we consider the total receptor concentration:Embedded ImagewhereEmbedded Imagesubstituting the above terms into Embedded Image gives the following operational model of allosterism at a receptor with nB orthosteric agonist binding sites (Eq 4 in the main text):Embedded ImagewhereEmbedded ImageandEmbedded Image

For use in Graphpad Prism or similar software, the above equation is described by the following notations, where nT or nB will likely need to be fixed to a known or theoretical value to fit real experimental data to this equation:Embedded ImageEmbedded ImageEmbedded ImageEmbedded ImageEmbedded Image

For the purposes of the current study, the model has been simplified further to enable analysis of data when the orthosteric agonist is a full agonist. Therefore, equation 4 in the manuscript reduces to:Embedded Image

Dividing through by Embedded Image, and defining [EC50]=KA/τA yields the following expression:Embedded Image

Authorship Contributions

Participated in research design: Gregory, Christopoulos, Leach.

Conducted experiments: Diao.

Contributed new reagents or analytic tools: Giraldo, Christopoulos.

Performed data analysis: Diao, Gregory, Giraldo, Leach.

Wrote or contributed to the writing of the manuscript: Gregory, Giraldo, Christopoulos, Leach.

Footnotes

    • Received August 12, 2019.
    • Accepted November 4, 2019.
  • This work was supported by the National Health and Medical Research Council of Australia [Project Grants APP1085143 (to K.L.), APP1138891 (to K.L. and K.J.G.), APP1084775 (to K.J.G.), and APP1127322 (to K.J.G.); Program Grant APP1055134 (to A.C.); and Senior Principal Research Fellowship APP1102950 (to A.C.)]; the Australian Research Council [Discovery Grant DP170104228 (to K.L. and K.J.G.), and Future Fellowships FT160100075 (to K.L.) and FT170100392 (to K.J.G.)]; and the Ministerio de Ciencia, Innovación y Universidades of Spain [Project Grant SAF2017-87199-R (to J.G.)].

  • https://doi.org/10.1124/mol.119.118091.

  • ↵Embedded ImageThis article has supplemental material available at molpharm.aspetjournals.org.

Abbreviations

AC265347
1-(1,3-benzothiazol-2-yl)-1-(2,4-dimethylphenyl)ethanol
ACh
acetylcholine
Ca2+i
intracellular calcium
Ca2+o
extracellular calcium
calcimimetic B
R-1-(6-methoxy-4′-(trifluoromethyl)-3-biphenylyl)-N-(R)- 1-phenylethyl)ethanamine
CaSR
calcium-sensing receptor
DMEM
Dulbecco’s modified Eagle’s medium
DPFE
1-(4-(2,4-difluorophenyl)piperazin-1-yl)-2-((4-fluorobenzyl)oxy)ethanone
Em
maximal system response
Emax
maximum effect
GPCR
G protein–coupled receptor
HEK
human embryonic kidney
KA
equilibrium dissociation constant of the orthosteric agonist
KB
equilibrium dissociation constant of the allosteric ligand
LY2033298
3-amino-5-chloro-6-methoxy-4-methyl-thieno[2,3-b]pyridine-2-carboxylic acid cyclopropylamide
M-5MPEP
2-[2-(3-methoxyphenyl)ethynyl]-5-methylpyridine
mAChR
muscarinic acetylcholine receptor
mGlu
metabotropic glutamate receptor
MPEP
2-methyl-6-(2-phenylethynyl)pyridine
NAM
negative allosteric modulator
nB
binding slope linking agonist concentration to receptor occupancy
nH
Hill slope
nT
transducer slope linking agonist concentration to response
NPS2143
2-chloro-6-[(2R)-2-hydroxy-3-[(2-methyl-1-naphthalen-2-ylpropan-2-yl)amino]propoxy]benzonitrile
PAM
positive allosteric modulator
WT
wild type
VE-29
2-chloro-6-[(2R)-2-hydroxy-3-[(2-methyl-1-naphthalen-2-ylpropan-2-yl)amino]propoxy]benzonitrile
  • Copyright © 2019 by The American Society for Pharmacology and Experimental Therapeutics

References

  1. ↵
    1. Abdul-Ridha A,
    2. Lane JR,
    3. Mistry SN,
    4. López L,
    5. Sexton PM,
    6. Scammells PJ,
    7. Christopoulos A, and
    8. Canals M
    (2014) Mechanistic insights into allosteric structure-function relationships at the M1 muscarinic acetylcholine receptor. J Biol Chem 289:33701–33711.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Aurelio L,
    2. Valant C,
    3. Flynn BL,
    4. Sexton PM,
    5. Christopoulos A, and
    6. Scammells PJ
    (2009) Allosteric modulators of the adenosine A1 receptor: synthesis and pharmacological evaluation of 4-substituted 2-amino-3-benzoylthiophenes. J Med Chem 52:4543–4547.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Black JW and
    2. Leff P
    (1983) Operational models of pharmacological agonism. Proc R Soc Lond B Biol Sci 220:141–162.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Black JW,
    2. Leff P,
    3. Shankley NP, and
    4. Wood J
    (1985) An operational model of pharmacological agonism: the effect of E/[A] curve shape on agonist dissociation constant estimation. Br J Pharmacol 84:561–571.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Brea J,
    2. Castro M,
    3. Giraldo J,
    4. López-Giménez JF,
    5. Padín JF,
    6. Quintián F,
    7. Cadavid MI,
    8. Vilaró MT,
    9. Mengod G,
    10. Berg KA, et al.
    (2009) Evidence for distinct antagonist-revealed functional states of 5-hydroxytryptamine2A receptor homodimers. Mol Pharmacol 75:1380–1391.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Brown EM
    (1983) Four-parameter model of the sigmoidal relationship between parathyroid hormone release and extracellular calcium concentration in normal and abnormal parathyroid tissue. J Clin Endocrinol Metab 56:572–581.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Christopoulos A
    (1998) Assessing the distribution of parameters in models of ligand-receptor interaction: to log or not to log. Trends Pharmacol Sci 19:351–357.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Cook AE,
    2. Mistry SN,
    3. Gregory KJ,
    4. Furness SG,
    5. Sexton PM,
    6. Scammells PJ,
    7. Conigrave AD,
    8. Christopoulos A, and
    9. Leach K
    (2015) Biased allosteric modulation at the CaS receptor engendered by structurally diverse calcimimetics. Br J Pharmacol 172:185–200.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Davey AE,
    2. Leach K,
    3. Valant C,
    4. Conigrave AD,
    5. Sexton PM, and
    6. Christopoulos A
    (2012) Positive and negative allosteric modulators promote biased signaling at the calcium-sensing receptor. Endocrinology 153:1232–1241.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Diepenhorst NA,
    2. Leach K,
    3. Keller AN,
    4. Rueda P,
    5. Cook AE,
    6. Pierce TL,
    7. Nowell C,
    8. Pastoureau P,
    9. Sabatini M,
    10. Summers RJ, et al.
    (2018) Divergent effects of strontium and calcium-sensing receptor positive allosteric modulators (calcimimetics) on human osteoclast activity. Br J Pharmacol 175:4095–4108.
    OpenUrl
  11. ↵
    1. Ehlert FJ
    (1988) Estimation of the affinities of allosteric ligands using radioligand binding and pharmacological null methods. Mol Pharmacol 33:187–194.
    OpenUrlAbstract
  12. ↵
    1. Geng Y,
    2. Mosyak L,
    3. Kurinov I,
    4. Zuo H,
    5. Sturchler E,
    6. Cheng TC,
    7. Subramanyam P,
    8. Brown AP,
    9. Brennan SC,
    10. Mun HC, et al.
    (2016) Structural mechanism of ligand activation in human calcium-sensing receptor. eLife 5:e13662.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Giraldo J
    (2015) Operational models of allosteric modulation: caution is needed. Trends Pharmacol Sci 36:1–2.
    OpenUrl
  14. ↵
    1. Huang Y,
    2. Zhou Y,
    3. Castiblanco A,
    4. Yang W,
    5. Brown EM, and
    6. Yang JJ
    (2009) Multiple Ca2+ binding sites in the extracellular domain of the Ca2+-sensing receptor corresponding to cooperative Ca2+ response. Biochemistry 48:388–398.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Keller AN,
    2. Kufareva I,
    3. Josephs TM,
    4. Diao J,
    5. Mai VT,
    6. Conigrave AD,
    7. Christopoulos A,
    8. Gregory KJ, and
    9. Leach K
    (2018) Identification of global and ligand-specific calcium sensing receptor activation mechanisms. Mol Pharmacol 93:619–630.
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Kenakin TP
    (2012) Biased signalling and allosteric machines: new vistas and challenges for drug discovery. Br J Pharmacol 165:1659–1669.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Leach K,
    2. Conigrave AD,
    3. Sexton PM, and
    4. Christopoulos A
    (2015) Towards tissue-specific pharmacology: insights from the calcium-sensing receptor as a paradigm for GPCR (patho)physiological bias. Trends Pharmacol Sci 36:215–225.
    OpenUrlCrossRef
  18. ↵
    1. Leach K,
    2. Davey AE,
    3. Felder CC,
    4. Sexton PM, and
    5. Christopoulos A
    (2011) The role of transmembrane domain 3 in the actions of orthosteric, allosteric, and atypical agonists of the M4 muscarinic acetylcholine receptor. Mol Pharmacol 79:855–865.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Leach K,
    2. Gregory KJ,
    3. Kufareva I,
    4. Khajehali E,
    5. Cook AE,
    6. Abagyan R,
    7. Conigrave AD,
    8. Sexton PM, and
    9. Christopoulos A
    (2016) Towards a structural understanding of allosteric drugs at the human calcium-sensing receptor. Cell Res 26:574–592.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Leach K,
    2. Loiacono RE,
    3. Felder CC,
    4. McKinzie DL,
    5. Mogg A,
    6. Shaw DB,
    7. Sexton PM, and
    8. Christopoulos A
    (2010) Molecular mechanisms of action and in vivo validation of an M4 muscarinic acetylcholine receptor allosteric modulator with potential antipsychotic properties. Neuropsychopharmacology 35:855–869.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Leach K,
    2. Sexton PM, and
    3. Christopoulos A
    (2007) Allosteric GPCR modulators: taking advantage of permissive receptor pharmacology. Trends Pharmacol Sci 28:382–389.
    OpenUrlCrossRefPubMed
  22. ↵
    1. May LT,
    2. Bridge LJ,
    3. Stoddart LA,
    4. Briddon SJ, and
    5. Hill SJ
    (2011) Allosteric interactions across native adenosine-A3 receptor homodimers: quantification using single-cell ligand-binding kinetics. FASEB J 25:3465–3476.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Moreno Delgado D,
    2. Møller TC,
    3. Ster J,
    4. Giraldo J,
    5. Maurel D,
    6. Rovira X,
    7. Scholler P,
    8. Zwier JM,
    9. Perroy J,
    10. Durroux T, et al.
    (2017) Pharmacological evidence for a metabotropic glutamate receptor heterodimer in neuronal cells. eLife 6:e25233.
    OpenUrl
  24. ↵
    1. Mutel V,
    2. Ellis GJ,
    3. Adam G,
    4. Chaboz S,
    5. Nilly A,
    6. Messer J,
    7. Bleuel Z,
    8. Metzler V,
    9. Malherbe P,
    10. Schlaeger EJ, et al.
    (2000) Characterization of [3H]quisqualate binding to recombinant rat metabotropic glutamate 1a and 5a receptors and to rat and human brain sections. J Neurochem 75:2590–2601.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Sato S,
    2. Huang XP,
    3. Kroeze WK, and
    4. Roth BL
    (2016) Discovery and characterization of novel GPR39 agonists allosterically modulated by zinc. Mol Pharmacol 90:726–737.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Sengmany K and
    2. Gregory KJ
    (2016) Metabotropic glutamate receptor subtype 5: molecular pharmacology, allosteric modulation and stimulus bias. Br J Pharmacol 173:3001–3017.
    OpenUrl
  27. ↵
    1. Sengmany K,
    2. Hellyer SD,
    3. Albold S,
    4. Wang T,
    5. Conn PJ,
    6. May LT,
    7. Christopoulos A,
    8. Leach K, and
    9. Gregory KJ
    (2019) Kinetic and system bias as drivers of metabotropic glutamate receptor 5 allosteric modulator pharmacology. Neuropharmacology 149:83–96.
    OpenUrl
  28. ↵
    1. Sengmany K,
    2. Singh J,
    3. Stewart GD,
    4. Conn PJ,
    5. Christopoulos A, and
    6. Gregory KJ
    (2017) Biased allosteric agonism and modulation of metabotropic glutamate receptor 5: implications for optimizing preclinical neuroscience drug discovery. Neuropharmacology 115:60–72.
    OpenUrl
  29. ↵
    1. Sinkins WG and
    2. Wells JW
    (1993) G protein-linked receptors labeled by [3H]histamine in guinea pig cerebral cortex. II. Mechanistic basis for multiple states of affinity [corrected] [published correction appears in Mol Pharmacol (1993) 44:1278]. Mol Pharmacol 43:583–594.
    OpenUrlAbstract
  30. ↵
    1. Stockton JM,
    2. Birdsall NJ,
    3. Burgen AS, and
    4. Hulme EC
    (1983) Modification of the binding properties of muscarinic receptors by gallamine. Mol Pharmacol 23:551–557.
    OpenUrlAbstract
  31. ↵
    1. Storjohann L,
    2. Holst B, and
    3. Schwartz TW
    (2008) Molecular mechanism of Zn2+ agonism in the extracellular domain of GPR39. FEBS Lett 582:2583–2588.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Vivo M,
    2. Lin H, and
    3. Strange PG
    (2006) Investigation of cooperativity in the binding of ligands to the D2 dopamine receptor. Mol Pharmacol 69:226–235.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Zhang C,
    2. Zhuo Y,
    3. Moniz HA,
    4. Wang S,
    5. Moremen KW,
    6. Prestegard JH,
    7. Brown EM, and
    8. Yang JJ
    (2014) Direct determination of multiple ligand interactions with the extracellular domain of the calcium-sensing receptor. J Biol Chem 289:33529–33542.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top

In this issue

Molecular Pharmacology: 97 (1)
Molecular Pharmacology
Vol. 97, Issue 1
1 Jan 2020
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Editorial Board (PDF)
  • Front Matter (PDF)
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Molecular Pharmacology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Evaluation of Operational Models of Agonism and Allosterism at Receptors with Multiple Orthosteric Binding Sites
(Your Name) has forwarded a page to you from Molecular Pharmacology
(Your Name) thought you would be interested in this article in Molecular Pharmacology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Research ArticleArticle

Operational Models for Receptors with Multiple Agonist Sites

Karen J. Gregory, Jesús Giraldo, Jiayin Diao, Arthur Christopoulos and Katie Leach
Molecular Pharmacology January 1, 2020, 97 (1) 35-47; DOI: https://doi.org/10.1124/mol.119.118091

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Share
Research ArticleArticle

Operational Models for Receptors with Multiple Agonist Sites

Karen J. Gregory, Jesús Giraldo, Jiayin Diao, Arthur Christopoulos and Katie Leach
Molecular Pharmacology January 1, 2020, 97 (1) 35-47; DOI: https://doi.org/10.1124/mol.119.118091
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Appendix
    • Authorship Contributions
    • Footnotes
    • Abbreviations
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • The binding site for KCI807 in the androgen receptor
  • Fatty acid amide hydrolase in cisplatin nephrotoxicity
  • eCB Signaling System in hiPSC-Derived Neuronal Cultures
Show more Articles

Similar Articles

Advertisement
  • Home
  • Alerts
Facebook   Twitter   LinkedIn   RSS

Navigate

  • Current Issue
  • Fast Forward by date
  • Fast Forward by section
  • Latest Articles
  • Archive
  • Search for Articles
  • Feedback
  • ASPET

More Information

  • About Molecular Pharmacology
  • Editorial Board
  • Instructions to Authors
  • Submit a Manuscript
  • Customized Alerts
  • RSS Feeds
  • Subscriptions
  • Permissions
  • Terms & Conditions of Use

ASPET's Other Journals

  • Drug Metabolism and Disposition
  • Journal of Pharmacology and Experimental Therapeutics
  • Pharmacological Reviews
  • Pharmacology Research & Perspectives
ISSN 1521-0111 (Online)

Copyright © 2023 by the American Society for Pharmacology and Experimental Therapeutics