Data-driven analysis in drug discovery

In the process of drug discovery for new chemical entities, application of appropriate pharmacological models often is not possible because the molecular mechanism of the compound is not yet elucidated. Therefore, a data-driven approach using generic tools designed to quantify characteristic patterns of concentration-response curves is required. This article outlines the options available for quantifying agonist and antagonist activity. Specifically, for agonists, the use of the Operational model for the determination of functional effects (equimolar potency ratios for full agonists, calculation of relative efficacy) is described. For antagonists, the measurement of pKB (-log of the equilibrium dissociation constant of the antagonist-receptor complex) for orthosteric antagonists that do not alter basal response (simple competitive antagonists), increase basal response (partial agonists), and decrease basal response (in constitutively active systems; inverse agonists) is discussed. In addition, this article considers methods to discern orthosteric receptor antagonism from allosteric antagonism whereby the agonist and antagonist bind to separate sites and interact through a conformational change in the receptor. Methods for the measurement of the pKB for allosteric modulators as well as co-operativity constants for these modulators is described.     

Allosteric and orthosteric sites in CC chemokine receptor (CCR5), a chimeric receptor approach

Chemokine receptors play a major role in immune system regulation and have consequently been targets for drug development leading to the discovery of several small molecule antagonists. Given the large size and predominantly extracellular receptor interaction of endogenous chemokines, small molecules often act more deeply in an allosteric mode. However, opposed to the well described molecular interaction of allosteric modulators in class C 7-transmembrane helix (7TM) receptors, the interaction in class A, to which the chemokine receptors belong, is more sparsely described. Using the CCR5 chemokine receptor as a model system, we studied the molecular interaction and conformational interchange required for proper action of various orthosteric chemokines and allosteric small molecules, including the well known CCR5 antagonists TAK-779, SCH-C, and aplaviroc, and four novel CCR5 ago-allosteric molecules. A chimera was successfully constructed between CCR5 and the closely related CCR2 by transferring all extracellular regions of CCR2 to CCR5, i.e. a Trojan horse that resembles CCR2 extracellularly but signals through a CCR5 transmembrane unit. The chimera bound CCR2 (CCL2 and CCL7), but not CCR5 chemokines (CCL3 and CCL5), with CCR2-like high affinities and potencies throughout the CCR5 signaling unit. Concomitantly, high affinity binding of small molecule CCR5 agonists and antagonists was retained in the transmembrane region. Importantly, whereas the agonistic and antagonistic properties were preserved, the allosteric enhancement of chemokine binding was disrupted. In summary, the Trojan horse chimera revealed that orthosteric and allosteric sites could be structurally separated and still act together with transmission of agonism and antagonism across the different receptor units.     

Identification and pharmacological characterization of domains involved in binding of CGRP receptor antagonists to the calcitonin-like receptor

The calcitonin-like receptor (CLR) and the calcitonin receptor (CTR) interact with receptor activity-modifying protein 1 (RAMP1) at the cell surface to form heterodimeric receptor complexes. CLR and CTR are members of the class II (family B) G-protein-coupled receptors (GPCR) and bind calcitonin gene-related peptide (CGRP) with similar affinities when coexpressed with RAMP1. The observation that various nonpeptide CGRP receptor antagonists display a higher affinity for the CLR/RAMP1 complex than for CTR/RAMP1 provided an opportunity to investigate the molecular determinants of the differential receptor affinities of these antagonists. A chimeric receptor approach was utilized to identify key domains within CLR responsible for conferring high-affinity antagonist binding. Initial chimera experiments implicated distinct regions within CLR as responsible for the affinities of structurally diverse CGRP receptor antagonists. Dissection of these key regions implicated amino acids 37-63 located in the amino terminus of CLR as responsible for the high-affinity interaction of one structural class, while transmembrane domain (TM) 7 was responsible for the interaction of a second class of antagonist. A unique binding interaction in the amino terminus of CLR is consistent with the observation that these compounds also interact with the extracellular region of RAMP1 and could suggest the formation of a binding pocket between the two proteins. Conversely, a compound which interacted with TM7 did not display a similar RAMP1 dependence, suggesting an allosteric mechanism of antagonism. Collectively, these data provide insight into two alternative mechanisms of antagonism for this unique heterodimeric receptor complex.     

Data-Driven Analysis in Drug Discovery

In the process of drug discovery for new chemical entities, application of appropriate pharmacological models often is not possible because the molecular mechanism of the compound is not yet elucidated. Therefore, a data-driven approach using generic tools designed to quantify characteristic patterns of concentration-response curves is required. This article outlines the options available for quantifying agonist and antagonist activity. Specifically, for agonists, the use of the Operational model for the determination of functional effects (equimolar potency ratios for full agonists, calculation of relative efficacy) is described. For antagonists, the measurement of pK_B (-log of the equilibrium dissociation constant of the antagonist-receptor complex) for orthosteric antagonists that do not alter basal response (simple competitive antagonists), increase basal response (partial agonists), and decrease basal response (in constitutively active systems; inverse agonists) is discussed. In addition, this article considers methods to discern orthosteric receptor antagonism from allosteric antagonism whereby the agonist and antagonist bind to separate sites and interact through a conformational change in the receptor. Methods for the measurement of the pK_B for allosteric modulators as well as co-operativity constants for these modulators is described.     

Regulation of signal transduction at M2 muscarinic receptor

Muscarinic acetylcholine receptors mediate transmission of an extracellular signal represented by released acetylcholine to neuronal or effector cells. There are five subtypes of closely homologous muscarinic receptors which are coupled by means of heterotrimeric G-proteins to a variety of signaling pathways resulting in a multitude of target cell effects. Endogenous agonist acetylcholine does not discriminate among individual subtypes and due to the close homology of the orthosteric binding site the same holds true for most of exogenous agonists. In addition to the classical binding site muscarinic receptors have one or more allosteric binding sites at extracellular domains. Binding of allosteric modulators induces conformational changes in the receptor that result in subtype-specific changes in orthosteric binding site affinity for both muscarinic agonists and antagonists. This overview summarizes our recent experimental effort in investigating certain aspects of M2 muscarinic receptor functioning concerning i) the molecular determinants that contribute to the binding of allosteric modulators, ii) G-protein coupling specificity and subsequent cellular responses and iii) possible functional assays that exploit the unique properties of allosteric modulators for characterization of muscarinic receptor subtypes in intact tissue. A detailed knowledge of allosteric properties of muscarinic receptors is required to permit drug design that will modulate signal transmission strength of specific muscarinic receptor subtypes. Furthermore, allosteric modulation of signal transmission strength is determined by cooperativity rather than concentration of allosteric modulator and thus reduces the danger of overdose.     

Inverse agonism and neutral antagonism at cannabinoid CB1 receptors

There are at least two types of cannabinoid receptor, CB1 and CB2, both G protein coupled. CB1 receptors are expressed predominantly at nerve terminals and mediate inhibition of transmitter release whereas CB2 receptors are found mainly on immune cells, one of their roles being to modulate cytokine release. Endogenous cannabinoid receptor agonists also exist and these "endocannabinoids" together with their receptors constitute the "endocannabinoid system". These discoveries were followed by the development of a number of CB1- and CB2-selective antagonists that in some CB1 or CB2 receptor-containing systems also produce "inverse cannabimimetic effects", effects opposite in direction from those produced by cannabinoid receptor agonists. This review focuses on the CB1-selective antagonists, SR141716A, AM251, AM281 and LY320135, and discusses possible mechanisms by which these ligands produce their inverse effects: (1) competitive surmountable antagonism at CB1 receptors of endogenously released endocannabinoids, (2) inverse agonism resulting from negative, possibly allosteric, modulation of the constitutive activity of CB1 receptors in which CB1 receptors are shifted from a constitutively active "on" state to one or more constitutively inactive "off" states and (3) CB1 receptor-independent mechanisms, for example antagonism of endogenously released adenosine at A1 receptors. Recently developed neutral competitive CB1 receptor antagonists, which are expected to produce inverse effects through antagonism of endogenously released endocannabinoids but not by modulating CB1 receptor constitutive activity, are also discussed. So too are possible clinical consequences of the production of inverse cannabimimetic effects, there being convincing evidence that released endocannabinoids can have "autoprotective" roles.     

Kinetic versus allosteric mechanisms to explain insurmountable antagonism and delayed ligand dissociation

The present review addresses the theories which have been advanced to explain experimental observations dealing with insurmountable antagonism and accelerated radioligand dissociation in the presence of an excess unlabelled ligand. We came to the perception that, for each of these phenomena, the theories can be placed into two distinctive categories. The "kinetic" interpretations attribute these phenomena to, respectively, the ability of antagonists to form long-lasting complexes with their cognate receptor and the ability of dissociated ligands to bind again to the same or neighbouring receptors rather than to diffuse away from the cell surface. On the other hand, these observations can also be explained by negative allosteric interactions among topographically distinct ligand binding sites at the same receptor or di/multimeric receptor complex.     

Allosteric modulation of kinases and GPCRs: design principles and structural diversity

Allosteric binding sites, as opposed to traditional orthosteric binding sites, offer unparalleled opportunities for drug discovery by providing high levels of selectivity, mimicking physiological conditions, affording fewer side effects because of desensitization/downregulation, and engendering ligands with chemotypes divergent from orthosteric ligands. For kinases, allosteric mechanisms described to date include alteration of protein kinase conformation blocking productive ATP binding which appear 'ATP competitive' or blocking kinase activation by conformational changes that are 'ATP non-competitive'. For GPCRs, allosteric mechanisms impart multiple modes of target modulation (positive allosteric modulation (PAM), negative allosteric modulation (NAM), neutral cooperativity, partial antagonism (PA), allosteric agonism and allosteric antagonism). Here, we review recent developments in the design principles and structural diversity of allosteric ligands for kinases and GPCRs.     

Allosteric alpha 1-adrenoreceptor antagonism by the conopeptide rho-TIA

A peptide contained in the venom of the predatory marine snail Conus tulipa, rho-TIA, has previously been shown to possess alpha1-adrenoreceptor antagonist activity. Here, we further characterize its pharmacological activity as well as its structure-activity relationships. In the isolated rat vas deferens, rho-TIA inhibited alpha1-adrenoreceptor-mediated increases in cytosolic Ca2+ concentration that were triggered by norepinephrine, but did not affect presynaptic alpha2-adrenoreceptor-mediated responses. In radioligand binding assays using [125I]HEAT, rho-TIA displayed slightly greater potency at the alpha 1B than at the alpha 1A or alpha 1D subtypes. Moreover, although it did not affect the rate of association for [3H]prazosin binding to the alpha 1B-adrenoreceptor, the dissociation rate was increased, indicating non-competitive antagonism by rho-TIA. N-terminally truncated analogs of rho-TIA were less active than the full-length peptide, with a large decline in activity observed upon removal of the fourth residue of rho-TIA (Arg4). An alanine walk of rho-TIA confirmed the importance of Arg4 for activity and revealed a number of other residues clustered around Arg4 that contribute to the potency of rho-TIA. The unique allosteric antagonism of rho-TIA resulting from its interaction with receptor residues that constitute a binding site that is distinct from that of the classical competitive alpha1-adrenoreceptor antagonists may allow the development of inhibitors that are highly subtype selective.     

Structure-activity relationship studies of carboxamido-biaryl ethers as opioid receptor antagonists (OpRAs). Part 1

A structurally unique and new class of opioid receptor antagonists (OpRAs) that bear no structural resemblance with morphine or endogenous opioid peptides has been discovered. A series of carboxamido-biaryl ethers were identified as potent receptor antagonists against mu, kappa and delta opioid receptors. The structure-activity relationship indicated para-substituted aryloxyaryl primary carboxamide bearing an amine tether on the distal phenyl ring was optimal for potent in vitro functional antagonism against three opioid receptor subtypes.     

Synthesis and SAR of a mGluR5 allosteric partial antagonist lead: unexpected modulation of pharmacology with slight structural modifications to a 5-(phenylethynyl)pyrimidine scaffold

This Letter describes the synthesis and SAR, developed through an iterative analogue library approach, of a mGluR5 allosteric partial antagonist lead based on a 5-(phenylethynyl)pyrimidine scaffold. With slight structural modifications to the distal phenyl ring, analogues demonstrated a range of pharmacological activities from mGluR5 partial antagonism to full antagonism/negative allosteric modulation to positive allosteric modulation.     

Further exploration of M1 allosteric agonists: Subtle structural changes abolish M1 allosteric agonism and result in pan-mAChR orthosteric antagonism

This letter describes the further exploration of two series of M₁ allosteric agonists, TBPB and VU0357017, previously reported from our lab. Within the TPBP scaffold, either electronic or steric perturbations to the central piperidine ring led to a loss of selective M₁ allosteric agonism and afforded pan-mAChR antagonism, which was demonstrated to be mediated via the orthosteric site. Additional SAR around a related M₁ allosteric agonist family (VU0357017) identified similar, subtle ‘molecular switches’ that modulated modes of pharmacology from allosteric agonism to pan-mAChR orthosteric antagonism. Therefore, all of these ligands are best classified as bi-topic ligands that possess high affinity binding at an allosteric site to engender selective M₁ activation, but all bind, at higher concentrations, to the orthosteric ACh site, leading to non-selective orthosteric site binding and mAChR antagonism.     

Opiate receptors and their interactions with agnoists and antagonists.

Independent evidence, derived from studies of the kinetics of receptor inactivation by sulfhydryl reagents, is presented in support of a previously postulated allosteric model of the opiate receptor, in which sodium ions act as effectors. Scatchard analysis provides evidence for cooperative binding, a characteristic of allosteric systems. Both new and old evidence is cited which is consistent with the hypothesis that agonists and antagonists bind to the same receptor sites. These sites can exist in alternative conformations that differ in their relative affinities for agonists and antagonists.     

Self-peptide/MHC and TCR antagonism: physiological role and therapeutic potential

TCR antagonists are peptides that bind MHC molecules and can specifically inhibit T cell activation induced by antigens. Studying TCR antagonism has taken an important place in immunology for both theoretical and practical reasons. Deciphering the mechanism(s) of action of TCR antagonists can yield important information about interactions of the TCR with ligands, T cell development, and TCR signaling. Moreover, microorganisms may employ TCR antagonism to elude the attention of the immune system. Finally, specificity of inhibition makes TCR antagonists an ideal tool to seek antigen-specific immunomodulation. Present state of knowledge on these topics is reviewed.     

The position of the double bond in monounsaturated free fatty acids is essential for the inhibition of the nicotinic acetylcholine receptor

Free fatty acids (FFAs) are non-competitive antagonists of the nicotinic acetylcholine receptor (AChR). Their site of action is supposedly located at the lipid-AChR interface. To elucidate the mechanism involved in this antagonism, we studied the effect that FFAs with a single double-bond at different positions (ω6, ω9, ω11 and ω13 cis-18:1) have on different AChR properties. Electrophysiological studies showed that only two FFAs (ω6 and ω9) reduced the duration of the channel open-state. The briefest component of the closed-time distribution remained unaltered, suggesting that ω6 and ω9 behave as allosteric blockers. Fluorescence resonance energy transfer studies indicated that all FFAs locate at the lipid-AChR interface, ω6 being restricted to annular sites and all others occupying non-annular sites. The perturbation of the native membrane order by FFAs was evaluated by DPH (1,6-diphenyl-1,3,5-hexatriene) and Laurdan fluorescence polarization studies, with the greatest decrease observed for ω9 and ω11. AChR conformational changes produced by FFAs present at the lipid bilayer were evaluated by fluorescence quenching studies of pyrene-labeled AChR and also using the AChR conformational-sensitive probe crystal violet. All cis-FFAs produced AChR conformational changes at the transmembrane level, but only ω9, ω11 and ω13 perturbed the resting state. Thus, the position and isomerism of the torsion angle of unsaturated FFAs are probably a key factor in terms of AChR blockage, suggesting that FFAs with a unique cis double bond at a superficial position inside the membrane directly inhibit AChR function by perturbing a potential conserved core structure for AChR gating at that level.     

Low affinity channel blocking (uncompetitive) NMDA receptor antagonists as therapeutic agents – toward an understanding of their favorable tolerability

Studies in experimental models have suggested that NMDA receptor antagonists may have utility in the treatment of a wide variety of neurological and psychiatric disorders. However, clinical trials have not been encouraging largely because the antagonists evaluated to date have exhibited unacceptable neurobehavioral side effects. In animals, therapeutic doses of some low-affinity channel blocking (uncompetitive) NMDA receptor antagonists are associated with less gross neurological impairment and behavioral toxicity than other types of NMDA receptor antagonists. Favorable clinical experiences with several such agents has bolstered confidence in the neurotherapeutic potential of low affinity NMDA antagonists. This article reviews current research attempting to explain the improved tolerability of such antagonists. While no single mechanism appears to account for the reduced toxicity of such agents, kinetic properties, particularly rapid blocking rate, seem to be of key importance. Other factors include partial trapping, reduced agonist-independent (closed channel) block, subunit selectivity (particularly for receptors that do not contain the NR2A subunit), combined block at allosteric (voltage-independent) sites, and synergistic therapeutic effects produced by additional actions at receptor targets apart from NMDA receptors (e.g., weak positive allosteric modulation of GABA(A) receptors or state-dependent Na+ channel block).     

G protein coupled receptors as allosteric proteins and the role of allosteric modulators

Seven transmembrane receptors (7TMRs) are proteins that convey signals through changes in conformation. These conformations are stabilized by external molecules (i.e. agonists, antagonists, modulators) and act upon other bodies (termed ‘guests’) which can be other molecules in the extracellular space, or proteins along the plane of the membrane (receptor oligomerization) or signaling proteins in the cytosol (i.e. G protein, β-arrestin). These elements comprise allosteric systems and a great deal of 7TMR pharmacology can be considered in terms of allosteric behavior. Allosteric ligands acting on 7TMRs possess four unique behaviors that can be valuable therapeutically; () the ability to alter the interaction of very large proteins, () probe dependence, () saturable effect, and () induction of separate changes in affinity and efficacy of other ligands. Two of these behaviors (namely probe dependence for CCR5-based HIV-1 entry inhibitors and functional selectivity for biased agonism) will be highlighted with examples.     

Dynamic control of allosteric antagonism of leukocyte function antigen-1 and intercellular adhesion molecule-1 interaction

Leukocyte function associated antigen-1 (LFA-1) plays a critical role in T cell migration and has been recognized as a therapeutic target for immune disorders. Several classes of small molecule antagonists have been developed to block LFA-1 interaction with intercellular adhesion molecule-1 (ICAM-1). Recent structural studies show that the antagonists bind to an allosteric site in the I-domain of LFA-1. However, it is not yet clear how these small molecules work as antagonists since no significant conformational change is observed in the I-domain-antagonist complex structures. Here we present a computational study suggesting how these allosteric antagonists affect the dynamics of the I-domain. The lowest frequency vibrational mode calculated from an LFA-1 I-domain structure shows large scale "coil-down" motion of the C-terminal alpha7 helix, which may lead to the open form of the I-domain. The presence of an allosteric antagonist greatly reduces this motion of the alpha7 helix as well as other parts of the I-domain. Thus, our study suggests that allosteric antagonists work by eliminating breathing motion that leads to the open conformation of the I-domain.     

Allosteric disulfide bonds

Disulfide bonds have been generally considered to be either structural or catalytic. Structural bonds stabilize a protein, while catalytic bonds mediate thiol-disulfide interchange reactions in substrate proteins. There is emerging evidence for a third type of disulfide bond that can control protein function by triggering a conformational change when it breaks and/or forms. These bonds can be thought of as allosteric disulfides. To better define the properties of allosteric disulfides, we have analyzed the geometry and dihedral strain of 6874 unique disulfide bonds in 2776 X-ray structures. A total of 20 types of disulfide bonds were identified in the dataset based on the sign of the five chi angles that make up the bond. The known allosteric disulfides were all contained in 1 of the 20 groups, the -RHStaple bonds. This bond group has a high mean potential energy and narrow energy distribution, which is consistent with a functional role. We suggest that the -RHStaple configuration is a hallmark of allosteric disulfides. About 1 in 15 of all structurally determined disulfides is a potential allosteric bond.     

3- and 6-Substituted 2-amino-4,5,6,7-tetrahydrothieno[2,3-c]pyridines as A1 adenosine receptor allosteric modulators and antagonists

A series of 2-amino-4,5,6,7-tetrahydrothieno[2,3-c]pyridines were prepared and evaluated as potential allosteric modulators at the A₁ adenosine receptor. The structure–activity relationships of the 3- and 6-positions of a series of 2-amino-4,5,6,7-tetrahydrothieno[2,3-c]pyridines were explored. Despite finding that 3- and 6-substituted 2-amino-4,5,6,7-tetrahydrothieno[2,3-c]pyridines possess the ability to recognize an allosteric site on the agonist-occupied A₁AR at relatively high concentrations, the structural modifications we have performed on this scaffold favor the expression of orthosteric antagonist properties over allosteric properties. This research has identified 2-amino-4,5,6,7-tetrahydrothieno[2,3-c]pyridines as novel class of orthosteric antagonist of the A₁AR and highlighted the close relationship between structural elements governing allosteric modulation and orthosteric antagonism of agonist function at the A₁AR.