Agonist-Specific Conformations of the M2 Muscarinic Acetylcholine Receptor Assessed by Molecular Dynamics.

Binding of muscarinic ligands, both antagonists and agonists, and their effects on the conformation of the M2 acetylcholine receptor were modeled in silico and compared to experimental data. After docking of antagonists to the M2 receptor in an inactive conformation (3UON, 5ZK3, 5ZKB, or 5ZKB) and agonists in an active conformation (4MQS), 100 ns of conventional molecular dynamics (MD) followed by 500 ns of accelerated MD was run. Conventional MD revealed ligand-specific interactions with the receptor. Antagonists stabilized the receptor in an inactive conformation during accelerated MD. The receptor in complex with various agonists attained different conformations specific to individual agonists. The magnitude of the TM6 movement correlated with agonist efficacy at the non-preferential Gs pathway. The shape of the intracellular opening where the receptor interacts with a G-protein was different for the classical agonist carbachol, super-agonist iperoxo, and Gi/o-biased partial agonists JR-6 and JR-7, being compatible with experimentally observed agonist bias at the G-protein level. Moreover, a wash-resistant binding of the unique agonist xanomeline associated with interactions with membrane lipids was formed during accelerated MD. Thus, accelerated MD is suitable for modeling of ligand-specific receptor binding and receptor conformations that is essential for the design of experiments aimed at identification of the secondary binding sites and understanding molecular mechanisms underlying receptor activation.

Structural Insight into Specificity of Interactions between Nonconventional Three-finger Weak Toxin from Naja kaouthia (WTX) and Muscarinic Acetylcholine Receptors.

Weak toxin from Naja kaouthia (WTX) belongs to the group of nonconventional "three-finger" snake neurotoxins. It irreversibly inhibits nicotinic acetylcholine receptors and allosterically interacts with muscarinic acetylcholine receptors (mAChRs). Using site-directed mutagenesis, NMR spectroscopy, and computer modeling, we investigated the recombinant mutant WTX analogue (rWTX) which, compared with the native toxin, has an additional N-terminal methionine residue. In comparison with the wild-type toxin, rWTX demonstrated an altered pharmacological profile, decreased binding of orthosteric antagonist N-methylscopolamine to human M1- and M2-mAChRs, and increased antagonist binding to M3-mAChR. Positively charged arginine residues located in the flexible loop II were found to be crucial for rWTX interactions with all types of mAChR. Computer modeling suggested that the rWTX loop II protrudes to the M1-mAChR allosteric ligand-binding site blocking the entrance to the orthosteric site. In contrast, toxin interacts with M3-mAChR by loop II without penetration into the allosteric site. Data obtained provide new structural insight into the target-specific allosteric regulation of mAChRs by "three-finger" snake neurotoxins.

Apolipoprotein E4 reduces evoked hippocampal acetylcholine release in adult mice.

Apolipoprotein E4 (apoE4) is the most prevalent genetic risk factor for Alzheimer's disease. We utilized apoE4-targeted replacement mice (approved by the Tel Aviv University Animal Care Committee) to investigate whether cholinergic dysfunction, which increases during aging and is a hallmark of Alzheimer's disease, is accentuated by apoE4. This revealed that levels of the pre-synaptic cholinergic marker, vesicular acetylcholine transporter in the hippocampus and the corresponding electrically evoked release of acetylcholine, are similar in 4-month-old apoE4 and apolipoprotein E3 (apoE3) mice. Both parameters decrease with age. This decrease is, however, significantly more pronounced in the apoE4 mice. The levels of cholinacetyltransferase (ChAT), acetylcholinesterase (AChE), and butyrylcholinesterase (BuChE) were similar in the hippocampus of young apoE4 and apoE3 mice and decreased during aging. For ChAT, this decrease was similar in the apoE4 and apoE3 mice, whereas it was more pronounced in the apoE4 mice, regarding their corresponding AChE and BuChE levels. The level of muscarinic receptors was higher in the apoE4 than in the apoE3 mice at 4 months and increased to similar levels with age. However, the relative representation of the M1 receptor subtype decreased during aging in apoE4 mice. These results demonstrate impairment of the evoked release of acetylcholine in hippocampus by apoE4 in 12-month-old mice but not in 4-month-old mice. The levels of ChAT and the extent of the M2 receptor-mediated autoregulation of ACh release were similar in the adult mice, suggesting that the apoE4-related inhibition of hippocampal ACh release in these mice is not driven by these parameters. Evoked ACh release from hippocampal and cortical slices is similar in 4-month-old apoE4 and apoE3 mice but is specifically and significantly reduced in hippocampus, but not cortex, of 12-month-old apoE4 mice. This effect is accompanied by decreased VAChT levels. These findings show that the hipocampal cholinergic nerve terminals are specifically affected by apoE4 and that this effect is age dependent.

Changes in Membrane Cholesterol Differentially Influence Preferential and Non-preferential Signaling of the M1 and M3 Muscarinic Acetylcholine Receptors.

We have found earlier that changes in membrane cholesterol content have distinct impact on signaling via the M1, M2, or M3 receptors expressed in CHO cells (CHO-M1 through CHO-M3). Now we investigated whether gradual changes in membrane cholesterol exerts differential effects on coupling of the M1 and M3 muscarinic receptors to preferential signaling pathways through Gq/11 and non-preferential Gs G-proteins signaling. Changes in membrane cholesterol resulted in only marginal alterations of antagonist and agonist affinity of the M1 and M3 receptors, and did not influence precoupling of either subtype. Changes in membrane cholesterol did not influence parameters of carbachol-stimulated GTP-γ(35)S binding in CHO-M1 membranes while reduction as well as augmentation of membrane cholesterol lowered the efficacy but increased the potency of carbachol in CHO-M3 membranes. Gradual increase or decrease in membrane cholesterol concentration dependently attenuated agonist-induced inositolphosphates release while only cholesterol depletion increased basal values in both cell lines. Similarly, membrane cholesterol manipulation modified basal and agonist-stimulated cAMP synthesis via Gs in the same way in both cell lines. These results demonstrate that changes in membrane cholesterol concentration differentially impact preferential and non-preferential M1 and M3 receptor signaling. They point to the activated G-protein/effector protein interaction as the main site of action in alterations of M1 receptor-mediated stimulation of second messenger pathways. On the other hand, modifications in agonist-stimulated GTP-γ(35)S binding in CHO-M3 membranes indicate that in this case changes in ligand-activated receptor/G-protein interaction may also play a role.

Classical and atypical agonists activate M1 muscarinic acetylcholine receptors through common mechanisms.

We mutated key amino acids of the human variant of the M1 muscarinic receptor that target ligand binding, receptor activation, and receptor-G protein interaction. We compared the effects of these mutations on the action of two atypical M1 functionally preferring agonists (N-desmethylclozapine and xanomeline) and two classical non-selective orthosteric agonists (carbachol and oxotremorine). Mutations of D105 in the orthosteric binding site and mutation of D99 located out of the orthosteric binding site decreased affinity of all tested agonists that was translated as a decrease in potency in accumulation of inositol phosphates and intracellular calcium mobilization. Mutation of D105 decreased the potency of the atypical agonist xanomeline more than that of the classical agonists carbachol and oxotremorine. Mutation of the residues involved in receptor activation (D71) and coupling to G-proteins (R123) completely abolished the functional responses to both classical and atypical agonists. Our data show that both classical and atypical agonists activate hM1 receptors by the same molecular switch that involves D71 in the second transmembrane helix. The principal difference among the studied agonists is rather in the way they interact with D105 in the orthosteric binding site. Furthermore, our data demonstrate a key role of D105 in xanomeline wash-resistant binding and persistent activation of hM1 by wash-resistant xanomeline.

Molecular mechanisms of methoctramine binding and selectivity at muscarinic acetylcholine receptors.

Methoctramine (N,N'-bis[6-[[(2-methoxyphenyl)-methyl]hexyl]-1,8-octane] diamine) is an M(2)-selective competitive antagonist of muscarinic acetylcholine receptors and exhibits allosteric properties at high concentrations. To reveal the molecular mechanisms of methoctramine binding and selectivity we took advantage of reciprocal mutations of the M(2) and M(3) receptors in the second and third extracellular loops that are involved in the binding of allosteric ligands. To this end we performed measurements of kinetics of the radiolabeled antagonists N-methylscopolamine (NMS) in the presence of methoctramine and its precursors, fluorescence energy transfer between green fluorescent protein-fused receptors and an Alexa-555-conjugated precursor of methoctramine, and simulation of molecular dynamics of methoctramine association with the receptor. We confirm the hypothesis that methoctramine high-affinity binding to the M(2) receptors involves simultaneous interaction with both the orthosteric binding site and the allosteric binding site located between the second and third extracellular loops. Methoctramine can bind solely with low affinity to the allosteric binding site on the extracellular domain of NMS-occupied M(2) receptors by interacting primarily with glutamate 175 in the second extracellular loop. In this mode, methoctramine physically prevents dissociation of NMS from the orthosteric binding site. Our results also demonstrate that lysine 523 in the third extracellular loop of the M(3) receptors forms a hydrogen bond with glutamate 219 of the second extracellular loop that hinders methoctramine binding to the allosteric site at this receptor subtype. Impaired interaction with the allosteric binding site manifests as low-affinity binding of methoctramine at the M(3) receptor.

Agonist-selective activation of individual G-proteins by muscarinic receptors.

Selective activation of individual subtypes of muscarinic receptors is a promising way to safely alleviate a wide range of pathological conditions in the central nervous system and the periphery as well. The flexible G-protein interface of muscarinic receptors allows them to interact with several G-proteins with various efficacy, potency, and kinetics. Agonists biased to the particular G-protein mediated pathway may result in selectivity among muscarinic subtypes and, due to the non-uniform expression of individual G-protein alpha subunits, possibly achieve tissue specificity. Here, we demonstrate that novel tetrahydropyridine-based agonists exert specific signalling profiles in coupling with individual G-protein α subunits. These signalling profiles profoundly differ from the reference agonist carbachol. Moreover, coupling with individual Gα induced by these novel agonists varies among subtypes of muscarinic receptors which may lead to subtype selectivity. Thus, the novel tetrahydropyridine-based agonist can contribute to the elucidation of the mechanism of pathway-specific activation of muscarinic receptors and serve as a starting point for the development of desired selective muscarinic agonists.

On homology modeling of the M2 muscarinic acetylcholine receptor subtype.

Twelve homology models of the human M2 muscarinic receptor using different sets of templates have been designed using the Prime program or the modeller program and compared to crystallographic structure (PDB:3UON). The best models were obtained using single template of the closest published structure, the M3 muscarinic receptor (PDB:4DAJ). Adding more (structurally distant) templates led to worse models. Data document a key role of the template in homology modeling. The models differ substantially. The quality checks built into the programs do not correlate with the RMSDs to the crystallographic structure and cannot be used to select the best model. Re-docking of the antagonists present in crystallographic structure and relative binding energy estimation by calculating MM/GBSA in Prime and the binding energy function in YASARA suggested it could be possible to evaluate the quality of the orthosteric binding site based on the prediction of relative binding energies. Although estimation of relative binding energies distinguishes between relatively good and bad models it does not indicate the best one. On the other hand, visual inspection of the models for known features and knowledge-based analysis of the intramolecular interactions allows an experimenter to select overall best models manually.

NMR structure and action on nicotinic acetylcholine receptors of water-soluble domain of human LYNX1.

Discovery of proteins expressed in the central nervous system sharing the three-finger structure with snake α-neurotoxins provoked much interest to their role in brain functions. Prototoxin LYNX1, having homology both to Ly6 proteins and three-finger neurotoxins, is the first identified member of this family membrane-tethered by a GPI anchor, which considerably complicates in vitro studies. We report for the first time the NMR spatial structure for the water-soluble domain of human LYNX1 lacking a GPI anchor (ws-LYNX1) and its concentration-dependent activity on nicotinic acetylcholine receptors (nAChRs). At 5-30 µM, ws-LYNX1 competed with (125)I-α-bungarotoxin for binding to the acetylcholine-binding proteins (AChBPs) and to Torpedo nAChR. Exposure of Xenopus oocytes expressing α7 nAChRs to 1 µM ws-LYNX1 enhanced the response to acetylcholine, but no effect was detected on α4ß2 and α3ß2 nAChRs. Increasing ws-LYNX1 concentration to 10 µM caused a modest inhibition of these three nAChR subtypes. A common feature for ws-LYNX1 and LYNX1 is a decrease of nAChR sensitivity to high concentrations of acetylcholine. NMR and functional analysis both demonstrate that ws-LYNX1 is an appropriate model to shed light on the mechanism of LYNX1 action. Computer modeling, based on ws-LYNX1 NMR structure and AChBP x-ray structure, revealed a possible mode of ws-LYNX1 binding.

Characterization of the Drosophila adenosine receptor: the effect of adenosine analogs on cAMP signaling in Drosophila cells and their utility for in vivo experiments.

Adenosine receptors (AR) belonging to the G protein-coupled receptor family influence a wide range of physiological processes. Recent elucidation of the structure of human A2AR revealed the conserved amino acids necessary for contact with the Ado moiety. However, the selectivity of Ado analogs for AR subtypes is still not well understood. We have shown previously that the Drosophila adenosine receptor (DmAdoR) evokes an increase in cAMP and calcium concentration in heterologous cells. In this study, we have characterized the second-messenger stimulation by endogenous DmAdoR in a Drosophila neuroblast cell line and examined a number of Ado analogs for their ability to interact with DmAdoR. We show that Ado can stimulate cAMP but not calcium levels in Drosophila cells. We found one full and four partial DmAdoR agonists, as well as four antagonists. The employment of the full agonist, 2-chloroadenosine, in flies mimicked in vivo the phenotype of DmAdoR over-expression, whereas the antagonist, SCH58261, rescued the flies from the lethality caused by DmAdoR over-expression. Differences in pharmacological effect of the tested analogs between DmAdoR and human A2AR can be partially explained by the dissimilarity of specific key amino acid residues disclosed by the alignment of these receptors.

A specific multi-nutrient formulation enhances M1 muscarinic acetylcholine receptor responses in vitro.

Recent evidence indicates that supplementation with a specific combination of nutrients may affect cell membrane synthesis and composition. To investigate whether such nutrients may also modify the physical properties of membranes, and affect membrane-bound processes involved in signal transduction pathways, we studied the effects of nutrient supplementation on G protein-coupled receptor activation in vitro. In particular, we investigated muscarinic receptors, which are important for the progression of memory deterioration and pathology of Alzheimer's disease. Nerve growth factor differentiated pheochromocytoma cells that were supplemented with specific combinations of nutrients showed enhanced responses to muscarinic receptor agonists in a membrane potential assay. The largest effects were obtained with a combination of nutrients known as Fortasyn™ Connect, comprising docosahexaenoic acid, eicosapentaenoic acid, uridine monophosphate as a uridine source, choline, vitamin B6, vitamin B12, folic acid, phospholipids, vitamin C, vitamin E, and selenium. In subsequent experiments, it was shown that the effects of supplementation could not be attributed to single nutrients. In addition, it was shown that the agonist-induced response and the supplement-induced enhancement of the response were blocked with the muscarinic receptor antagonists atropine, telenzepine, and AF-DX 384. In order to determine whether the effects of Fortasyn™ Connect supplementation were receptor subtype specific, we investigated binding properties and activation of human muscarinic M1, M2 and M4 receptors in stably transfected Chinese hamster ovary cells after supplementation. Multi-nutrient supplementation did not change M1 receptor density in plasma membranes. However, M1 receptor-mediated G protein activation was significantly enhanced. In contrast, supplementation of M2- or M4-expressing cells did not affect receptor signaling. Taken together, these results indicate that a specific combination of nutrients acts synergistically in enhancing muscarinic M1 receptor responses, probably by facilitating receptor-mediated G protein activation.

Neurosteroids and steroid hormones are allosteric modulators of muscarinic receptors.

The membrane cholesterol was found to bind and modulate the function of several G-protein coupled receptors including muscarinic acetylcholine receptors. We investigated the binding of 20 steroidal compounds including neurosteroids and steroid hormones to muscarinic receptors. Corticosterone, progesterone and some neurosteroids bound to muscarinic receptors with the affinity of 100 nM or greater. We established a structure-activity relationship for steroid-based allosteric modulators of muscarinic receptors. Further, we show that corticosterone and progesterone allosterically modulate the functional response of muscarinic receptors to acetylcholine at physiologically relevant concentrations. It can play a role in stress control or in pregnancy, conditions where levels of these hormones dramatically oscillate. Allosteric modulation of muscarinic receptors via the cholesterol-binding site represents a new pharmacological approach at diseases associated with altered cholinergic signalling.

Neuroactive steroids, WIN-compounds and cholesterol share a common binding site on muscarinic acetylcholine receptors.

Endogenous neurosteroids and their synthetic analogues-neuroactive steroids-have been found to bind to muscarinic acetylcholine receptors and allosterically modulate acetylcholine binding and function. Using radioligand binding experiments we investigated their binding mode. We show that neuroactive steroids bind to two binding sites on muscarinic receptors. Their affinity for the high-affinity binding site is about 100 nM. Their affinity for the low-affinity binding site is about 10 µM. The high-affinity binding occurs at the same site as binding of steroid-based WIN-compounds that is different from the common allosteric binding site for alcuronium or gallamine that is located between the second and third extracellular loop of the receptor. This binding site is also different from the allosteric binding site for the structurally related aminosteroid-based myorelaxants pancuronium and rapacuronium. Membrane cholesterol competes with neurosteroids/neuroactive steroids binding to both high- and low-affinity binding site, indicating that both sites are oriented towards the cell membrane..

Functional cholinergic damage develops with amyloid accumulation in young adult APPswe/PS1dE9 transgenic mice.

We investigated the functional characteristics of pre- and postsynaptic cholinergic transmission in APPswe/PS1dE9 double transgenic mice at a young age (7-10 weeks) before the onset of amyloid plaque formation and at adult age (5-6 months) at its onset. We compared brain slices from cerebral cortex and hippocampus with amyloid deposits to slices from striatum with no amyloid plaques by 6 months of age. In young transgenic mice we found no impairments of preformed and newly synthesized [(3)H]-ACh release, indicating intact releasing machinery and release turnover, respectively. Adult transgenic mice displayed a significant increase in preformed [(3)H]-ACh release in cortex but a decrease in hippocampus and striatum. The extent of presynaptic muscarinic autoregulation was unchanged. Evoked release of newly synthesized [(3)H]-ACh was significantly reduced in the cortex and hippocampus but unchanged in the striatum. Carbachol-induced G-protein activation in cortical membranes displayed decreased potency but normal efficacy in adult animals and no changes in young animals. These results indicate that functional pre- and postsynaptic cholinergic deficits are not present in APPswe/PS1dE9 transgenic mice before 10 weeks of age, but develop along with beta-amyloid accumulation in the brain.

The operational model of allosteric modulation of pharmacological agonism.

Proper determination of agonist efficacy is indispensable in the evaluation of agonist selectivity and bias to activation of specific signalling pathways. The operational model (OM) of pharmacological agonism is a useful means for achieving this goal. Allosteric ligands bind to receptors at sites that are distinct from those of endogenous agonists that interact with the orthosteric domain on the receptor. An allosteric modulator and an orthosteric agonist bind simultaneously to the receptor to form a ternary complex, where the allosteric modulator affects the binding affinity and operational efficacy of the agonist. Allosteric modulators are an intensively studied group of receptor ligands because of their selectivity and preservation of physiological space-time pattern of the signals they modulate. We analysed the operational model of allosterically-modulated agonism (OMAM) including modulation by allosteric agonists. Similar to OM, several parameters of OMAM are inter-dependent. We derived equations describing mutual relationships among parameters of the functional response and OMAM. We present a workflow for the robust fitting of OMAM to experimental data using derived equations.

Novel M2 -selective, Gi -biased agonists of muscarinic acetylcholine receptors.

BACKGROUND AND PURPOSE: More than 30% of currently marketed medications act via GPCRs. Thus, GPCRs represent one of the most important pharmacotherapeutic targets. In contrast to traditional agonists activating multiple signalling pathways, agonists activating a single signalling pathway represent a new generation of drugs with increased specificity and fewer adverse effects. EXPERIMENTAL APPROACH: We have synthesized novel agonists of muscarinic ACh receptors and tested their binding and function (on levels of cAMP and inositol phosphates) in CHO cells expressing individual subtypes of muscarinic receptors, primary cultures of rat aortic smooth muscle cells and suspensions of digested native tissues from rats. Binding of the novel compounds to M2 receptors was modelled in silico. KEY RESULTS: Two of the tested new compounds (1-(thiophen-2-ylmethyl)-3,6-dihydro-2H-pyridinium and 1-methyl-1-(thiophen-2-ylmethyl)-3,6-dihydro-2H-pyridinium) only inhibited cAMP synthesis in CHO cells, primary cultures, and native tissues, with selectivity for M2 muscarinic receptors and displaying bias towards the Gi signalling pathway at all subtypes of muscarinic receptors. Molecular modelling revealed interactions with the orthosteric binding site in a way specific for a given agonist followed by agonist-specific changes in the conformation of the receptor. CONCLUSIONS AND IMPLICATIONS: The identified compounds may serve as lead structures in the search for novel non-steroidal and non-opioid analgesics acting via M2 and M4 muscarinic receptors with reduced side effects associated with activation of the phospholipase C signalling pathway.

Detection of choline transporter-like 1 protein CTL1 in neuroblastoma x glioma cells and in the CNS, and its role in choline uptake.

Choline is an essential nutrient necessary for synthesis of membrane phospholipids, cell signalling molecules and acetylcholine. The aim of this study was to detect and characterize the choline transporter-like 1 (CTL1/SLC44A1) protein in CNS tissues and the hybrid neuroblastoma x glioma cell line NG108-15, which synthesizes acetylcholine and has high affinity choline transport but does not express the cholinergic high affinity choline transporter 1. The presence of CTL1 protein in NG108-15 cells was confirmed using our antibody G103 which recognizes the C-terminal domain of human CTL1. Three different cognate small interfering RNAs were used to decrease CTL1 mRNA in NG108-15 cells, causing lowered CTL1 protein expression, choline uptake and cell growth. None of the small interfering RNAs influenced carnitine transport, demonstrating the absence of major non-specific effects. In parental C6 cells knockdown of CTL1 also reduced high affinity choline transport. Our results support the concept that CTL1 protein is necessary for the high affinity choline transport which supplies choline for cell growth. The presence of CTL1 protein in rat and human CNS regions, where it is found in neuronal, glial and endothelial cells, suggests that malfunction of this transporter could have important implications in nervous system development and repair following injury, and in neurodegenerative diseases.

Divergence of allosteric effects of rapacuronium on binding and function of muscarinic receptors.

BACKGROUND: Many neuromuscular blockers act as negative allosteric modulators of muscarinic acetylcholine receptors by decreasing affinity and potency of acetylcholine. The neuromuscular blocker rapacuronium has been shown to have facilitatory effects at muscarinic receptors leading to bronchospasm. We examined the influence of rapacuronium on acetylcholine (ACh) binding to and activation of individual subtypes of muscarinic receptors expressed in Chinese hamster ovary cells to determine its receptor selectivity. RESULTS: At equilibrium rapacuronium bound to all subtypes of muscarinic receptors with micromolar affinity (2.7-17 microM) and displayed negative cooperativity with both high- and low-affinity ACh binding states. Rapacuronium accelerated [3H]ACh association with and dissociation from odd-numbered receptor subtypes. With respect to [35S]GTPgammaS binding rapacuronium alone behaved as an inverse agonist at all subtypes. Rapacuronium concentration-dependently decreased the potency of ACh-induced [35S]GTPgammaS binding at M2 and M4 receptors. In contrast, 0.1 microM rapacuronium significantly increased ACh potency at M1, M3, and M5 receptors. Kinetic measurements at M3 receptors showed acceleration of the rate of ACh-induced [35S]GTPgammaS binding by rapacuronium. CONCLUSIONS: Our data demonstrate a novel dichotomy in rapacuronium effects at odd-numbered muscarinic receptors. Rapacuronium accelerates the rate of ACh binding but decreases its affinity under equilibrium conditions. This results in potentiation of receptor activation at low concentrations of rapacuronium (1 microM) but not at high concentrations (10 microM). These observations highlight the relevance and necessity of performing physiological tests under non-equilibrium conditions in evaluating the functional effects of allosteric modulators at muscarinic receptors. They also provide molecular basis for potentiating M3 receptor-mediated bronchoconstriction.

Analysis of equilibrium binding of an orthosteric tracer and two allosteric modulators.

Allosteric ligands bind to receptors at sites that are distinct from those endogenous agonists and orthosteric pharmacological agents interact with. Both an allosteric and orthosteric ligand bind simultaneously to the receptor to form a ternary complex, where each ligand influences binding affinity of the other to the receptor, either positively or negatively. Allosteric modulators are an intensively studied group of receptor ligands because of their potentially greater selectivity over orthosteric ligands, with the possibility of fine tuning of the effects of endogenous neurotransmitters and hormones. The affinity of an unlabelled allosteric ligand is commonly estimated by measuring its effects on binding of a radio-labelled orthosteric tracer. This scenario is complicated by many folds when one studies the kinetics of interactions of two allosteric agents, added simultaneously, on binding of an orthosteric tracer. In this paper, we provide, for the first time, theoretical basis for analysis of such complex interactions. We have expanded our analysis to include the possibility of having two allosteric modulators interact with the same or different sites on the receptor. An added value of our analysis is to provide a tool to distinguish between the two situations. Finally, we also modelled binding of two molecules of one allosteric modulator to one receptor.

Applications and limitations of fitting of the operational model to determine relative efficacies of agonists.

Proper determination of agonist efficacy is essential in the assessment of agonist selectivity and signalling bias. Agonist efficacy is a relative term that is dependent on the system in which it is measured, especially being dependent on receptor expression level. The operational model (OM) of functional receptor agonism is a useful means for the determination of agonist functional efficacy using the maximal response to agonist and ratio of agonist functional potency to its equilibrium dissociation constant (KA) at the active state of the receptor. However, the functional efficacy parameter τ is inter-dependent on two other parameters of OM; agonist's KA and the highest response that could be evoked in the system by any stimulus (EMAX). Thus, fitting of OM to functional response data is a tricky process. In this work we analyse pitfalls of fitting OM to experimental data and propose a rigorous fitting procedure where KA and EMAX are derived from half-efficient concentration of agonist and apparent maximal responses obtained from a series of functional response curves. Subsequently, OM with fixed KA and EMAX is fitted to functional response data to obtain τ. The procedure was verified at M2 and M4 muscarinic receptors fused with the G15 G-protein α-subunit. The procedure, however, is applicable to any receptor-effector system.