Bitopic Ligands and Metastable Binding Sites: Opportunities for G

Jan 31, 2017 - However, it remains a challenge to develop subtype selective GPCR ligands due to the high conservation of orthosteric binding sites. Bi...
4 downloads 30 Views 3MB Size
Perspective pubs.acs.org/jmc

Bitopic Ligands and Metastable Binding Sites: Opportunities for G Protein-Coupled Receptor (GPCR) Medicinal Chemistry Philipp Fronik,† Birgit I. Gaiser,† and Daniel Sejer Pedersen* Department of Drug Design and Pharmacology, Faculty of Health and Medical Sciences, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark ABSTRACT: G protein-coupled receptors (GPCRs) belong to a large superfamily of membrane receptors mediating a variety of physiological functions. As such they are attractive targets for drug therapy. However, it remains a challenge to develop subtype selective GPCR ligands due to the high conservation of orthosteric binding sites. Bitopic ligands have been employed to address the selectivity problem by combining (linking) an orthosteric ligand with an allosteric modulator, theoretically leading to high-affinity subtype selective ligands. However, it remains a challenge to identify suitable allosteric binding sites. Computational studies on ligand binding to GPCRs have revealed transient, low-affinity binding sites, termed metastable binding sites. Metastable binding sites may provide a new source of allosteric binding sites that could be exploited in the design of bitopic ligands. Unlike the bitopic ligands that have been reported to date, this type of bitopic ligands would be composed of two identical pharmacophores. Herein, we outline the concept of bitopic ligands, review metastable binding sites, and discuss their potential as a new source of allosteric binding sites.



INTRODUCTION G protein-coupled receptors (GPCRs), also referred to as seven-transmembrane (7TM) receptors, belong to a large superfamily of receptors that are divided into five families, termed: rhodopsin, secretin, glutamate, adhesion, and frizzled/ TAS2.1 GPCRs respond to a wide variety of ligands ranging from photons, ions, and small molecules to peptides.2 The physiological functions of GPCRs are wide-ranging and include smell, taste, vision, and the reaction to neurotransmitters and hormones that affect cardiovascular, neurological, and reproductive functions among others.3,4 The abundance of physiological functions of GPCRs has made them an attractive target for drug discovery, and approximately 30% of all FDAapproved drugs act on GPCRs.1,5 However, upon closer inspection of marketed GPCR drugs, it becomes apparent that the vast majority of these drugs target the aminergic receptors belonging to the rhodopsin family. Thus, to date medicinal chemists have only scratched the surface with respect to exploiting the full therapeutic potential of GPCRs. The biggest challenge facing medicinal chemists is to develop subtype selective GPCR ligands. Orthosteric binding sites are highly conserved and hence do not offer many options in terms of selectivity. Allosteric binding sites in contrast are generally distant from the orthosteric site and have been altered significantly throughout evolutionary processes. Targeting allosteric sites to achieve subtype selectivity has been a successful strategy, but it is hampered by the lack of knowledge regarding their location and structure.6,7 Thus, any advances in the prediction and identification of new allosteric binding sites are highly desirable. Bitopic ligands are yet another class of © 2017 American Chemical Society

ligands introduced to address the subtype selectivity problem. They contain two distinct pharmacophores which are connected by a linker, thus allowing concomitant binding to the orthosteric and an allosteric binding site of the same receptor monomer. In addition to the benefit of inducing selectivity into the orthosteric ligand, bitopic ligands can provide advantages such as higher affinity, improved off-rates, etc.,8−13 and can induce signaling bias (i.e., functional selectivity).14,15 With increasing computational power, computational chemists have undertaken extensive molecular dynamics simulations to probe the mechanism of action for an orthosteric ligand’s entry and exit from the binding pocket. Several independent reports have shown that various orthosteric ligands often engage in stable interactions at the entrance of the receptor and adopt energetically stable conformations before either continuing toward the binding pocket or dissociating. These transient binding sites have been termed metastable binding sites (or meta-binding sites).16−21 Metastable binding sites may provide a new source of allosteric binding sites that could be exploited in the design of allosteric modulators and bitopic ligands. Herein, we provide an overview of bitopic ligands targeting different GPCR subfamilies to showcase the opportunities the approach provides within medicinal chemistry. Moreover, the current knowledge on metastable binding sites and how these potential allosteric binding sites could be used in the design of bitopic ligands is discussed. Received: October 28, 2016 Published: January 31, 2017 4126

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

Figure 1. Bitopic/dualsteric ligands target the orthosteric and an allosteric binding site in a single receptor and commonly are heterobivalent. Dimeric ligands target two separate receptors/protomers simultaneously, here illustrated for a GPCR homodimer. Metastable binding sites may provide a new source of allosteric binding sites that could be targeted by homobivalent bitopic ligands.



BITOPIC LIGANDS Bivalent ligands are defined as consisting of two pharmacophores that simultaneously engage a combination of two allosteric and/or orthosteric binding sites.22−24 Dimeric ligands (ligands that bind to two separate receptors) and bitopic ligands (also termed dualsteric) are subclasses of bivalent ligands (Figure 1). The term bitopic ligand describes a molecule able to simultaneously bind to both the orthosteric and an allosteric binding site of a single receptor monomer. The two pharmacophores are connected by a linker adapted to fit the target. This idea dates back to Porthogese’s messageaddress concept that was applied to the development of opioid receptor dimer antagonists.25 All bivalent ligands can be further categorized as homo- or heterobivalent, depending on whether the same pharmacophore is engaging both binding sites or not.26−29 Bitopic ligands in theory combine the advantages of orthosteric and allosteric ligands by combining high affinity and subtype selective pharmacophores, respectively (Figure 2).

can be designed as orthosteric (inverse) agonists or antagonists connected to either positive or negative allosteric modulators (PAMs or NAMs), in theory providing access to many regulatory options. However, it is difficult to predict the outcome for such ligands. Chen et al. combined a PAM of the M1-muscarinic acetylcholine receptor (mAChR) with a full agonist, but rather than increasing potency a partial agonist was created.31 Mismatches such as a NAM-agonist hybrid can result in partial loss of affinity, presumably because each pharmacophore stabilizes a different receptor conformation.32 The M1-Muscarinic Acetylcholine Receptor (mAChR). The history of bitopic ligands dates back to 1961, when 1 (McN-A-343, Figure 3) was first described.34 At the time, the

Figure 3. Examples of bitopic mAChR ligands 1,34 2,37 3,38 4,39 5,40 6,41 and 7.42

Figure 2. Generally observed trends of orthosteric, allosteric, and bitopic ligands. Inspired by Kruse et al.33

ligand was believed to be a selective agonist of the mAChR M1 subtype. In the 1980s, evidence suggested that the compound was an allosteric ligand,35 however, recently it was shown that 1 in fact is bitopic. Progressively truncating 1 led to the identification of both an orthosteric and an allosteric fragment, and the pharmacological results were supported by molecular modeling studies. The allosteric fragment was found to positively modulate orthosteric antagonists, while negatively modulating orthosteric agonists, and was considered respon-

Moreover, bitopic ligands bear the potential to induce biased signaling, which for example has been shown for an adenosine A1 receptor ligand, an orthosteric-allosteric hybrid derived from the unbiased orthosteric ligand adenosine (vide infra).30 When designing bitopic ligands, both the orthosteric and the allosteric fragments need to be selected carefully and the linker must be optimized for the target with respect to its attachment point, length, composition, and flexibility.12,13 Bitopic ligands 4127

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

A1AR ligands are rarely used as drugs because the same receptor mediates bradycardia as a significant side effect.43 Thus, it is desirable to develop biased ligands that will induce the beneficial effect of cytoprotection while eliminating bradycardia side effects. The first bitopic ligand for A1AR was reported in 2010 in an effort to locate the allosteric binding site of the receptor.57 The agonist adenosine was linked to a PAM via the exocyclic amino group of the nucleobase, and the linker length was increased systematically. The bitopic ligands were evaluated pharmacologically in orthosteric radioligand displacement studies and in functional assays in the presence and absence of the PAM. The ligand 8 (LUF6258) with a linker length of nine atoms displayed the best overall properties (Figure 4).57 Later, Lane et al. reanalyzed the original data on 8

sible for both the partial agonism and the functional selectivity of 1. It was further hypothesized that 1 in addition to its bitopic mode of action acts as an allosteric modulator by only binding to the allosteric binding site when the orthosteric pocket is already occupied by another ligand. On the basis of these findings, 1 was classified as a bitopic, functionally selective mAChR partial agonist.36 In recent years, several new and selective mAChR agonists have been discovered, particularly for the M1 subtype. Just as in the case of 1, the selectivity for these ligands were ascribed to targeting allosteric binding sites, but later a bitopic mode of action was revealed. When 2 (AC-42)37 and 3 (77-LH-28-1)38 were discovered, they were believed to have a bitopic mode of action because their pharmacological profiles were similar to that of 1.36,43 However, this was only confirmed several years later based on further pharmacological evidence and molecular modeling.44 As in the case of 1, these ligands were found to adopt different orientations depending on whether an orthosteric ligand was present or not. The bitopic mode of action of 2 was probed and substantiated further by FRET experiments with fluorescent 2 derivatives and a fluorescent protein-fused mAChR.45,46 Further examples of such ligands are the selective M1 mAChR agonists 4 (VU0364572),39 5 (VU0357017),40 and several analogues, which were found to be bitopic during an optimization study.47,48 Likewise, 6 (TBPB) was first described as a novel, highly selective allosteric M1 mAChR agonist41 but later shown to be a bitopic ligand by a truncating study similar to that reported for 1.15,49 Since the discovery of their bitopic mode of action, the above-mentioned ligands have been used as tool compounds to investigate mAChR pharmacology. Site-directed mutagenesis, analytical pharmacology, and molecular modeling studies with 1−3 and 6 led to the identification of key receptor areas responsible for selective binding, signaling efficacy, and pathway-specific receptor conformations, i.e., functional selectivity.49−51 Taken together, these studies highlight the potential of targeting allosteric pockets, which can promote unique receptor states and provide biased ligands.51 The knowledge gained from the discovery of bitopic ligands and their pharmacological characterization has later been applied in the rational design of bitopic ligands. Disingrini et al. reported the first de novo synthesis of ligands that were designed to target the orthosteric and an allosteric pocket simultaneously.52 They combined nonselective, high affinity agonists with subtype selective allosteric ligands and obtained a new class of agonists with a pharmacological profile similar to that of the allosteric parent compound and described them as hybrid compounds. Several examples of rationally designed bitopic ligands have followed since then.31,42,53 One of these, 7 (THRX-160209), nicely illustrates that careful design, i.e., the appropriate combination of orthosteric and allosteric fragments, and choice of linker, can give rise to ligands with notable improvements in affinity and selectivity.42 To study receptor conformational transformations Bock et al. designed a series of bitopic agonists as chemical probes and found that the allosteric vestibule of GPCRs controls their functional selectivity.54 In 2014, the principle was extended to include inactive receptors when the first antagonist hybrid derived from an orthosteric antagonist and an allosteric antagonist was reported.55 The Adenosine A1 Receptor (A1AR). The adenosine A1 receptor (A1AR) is an important target within tissue protection, in particular in the context of cardiac diseases.56 However,

Figure 4. Bitopic A1AR ligands 857 and 9,30 both composed by the orthosteric agonist adenosine (blue) linked to a PAM (green).

and its analogues with shorter linkers and found a clear signaling bias toward cytoprotection for 8 in comparison to its shorter counterparts.29 Hence this study demonstrates that by careful design of bitopic ligands it is possible to access unique receptor conformations that can result in functionally selective compounds. This study was followed up by Valant et al., who developed 9,43 which fused adenosine with the PAM (2-amino4-(2-(trifluoromethyl)phenyl)thiophen-3-yl)(phenyl)methanone (VCP171, structure not disclosed).30 The latter was a known ligand that had been shown to bias cytoprotection at A1AR, and it was hoped that the bradycardia signaling pathway could be reduced further in this manner. As anticipated, the bitopic ligand 9 was found to induce significant biased cellular signaling in favor of cytoprotection. The finding was confirmed in a more recent study in which 9 was profiled pharmacologically against a series of known agonists.14 The Dopamine D2-Like Receptor Family. The dopamine D2-like family is associated with many neuropsychiatric disorders as well as drug addiction, and in particular the D3 receptor (D3R) has been investigated extensively as a drug target.58−61 The first selective bitopic ligand 10 (SB269652, Figure 5) for D3R was published by Stemp et al. in 2000, but its bitopic mode of action was only determined recently.62,63 In the same study, 10 was shown to exert allostery at D2R as well as across a D2R dimer. This led to the hypothesis that bitopic 4128

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

Figure 6. 1376 and 1474 are bitopic ligands for the S1P3 receptor.



METASTABLE BINDING SITES During the past decade, it has become possible to perform wide-spanning molecular dynamics (MD) simulations that include entire receptors, cell membranes, and the surrounding solvent. Consequently, many computational studies that characterize ligand recognition as well as entry and exit pathways from a protein’s binding site have been reported. On the basis of MD simulations, it has on multiple occasions been reported that ligands and proteins form loose, transient complexes before the ligand enters the binding site.16,18−21,77−80 These transient binding sites have been referred to as “metastable binding sites” or “meta-binding sites”. Moreover, the resulting complexes have been termed “metastable state”, “loose binding”, “persistent intermediate conformations”, etc.19−21,77,80,81 The current hypothesis regarding metastable binding sites is that they function as selectivity filters. Thus, upon binding of a suitable ligand to the metastable binding site, the receptor will undergo a conformational change that allows the ligand to move further into the receptor to engage the orthosteric binding pocket. Other molecules that do not have the required features are denied access to the receptor and diffuse into the surroundings.82,83 Metastable Binding Sites in GPCRs. The existence of metastable binding sites was first suggested in 1999.84 However, it took more than a decade before experimental evidence of metastable binding sites was provided by Dror et al. based on MD simulations on the β1- and β2-adrenergic receptors (β1AR and β2AR).16 These two receptors represent some of the best characterized GPCRs with >15 crystal structures of the receptors alone or in complex with an agonist or antagonist reported.85,86 Dror et al. found that the structurally related βblockers alprenolol, dihydroalprenolol (15), and propranolol, and the β-agonist isoproterenol upon entering the β2AR occupy a similar metastable binding site in the extracellular vestibule of the receptor (Figure 7).16 The ammonium groups of the ligands were found to make a key interaction with Asp300ECL3 and resided in the vestibule for several hundred nanoseconds (pose 2) before moving on to the orthosteric binding pocket (pose 5). When the same MD simulation was performed with 15 on the β1AR, a similar result to that observed for β2AR was obtained. By a similar method, González et al. studied ligand exit from the β2AR by the antagonist/partial inverse agonist carazolol.17 The existence of a metastable binding site in the extracellular vestibule and the involvement of Asp300ECL3 among other key residues was recapitulated. Two channels for ligand entry for the endogenous agonist histamine at the histamine H4 receptor

Figure 5. Structures of bitopic dopamine receptor ligands 10,63 11,66 and 12.67 In the dopamine area, some research groups have defined PAMs and NAMs as allosteric ligands, whereas other binding sites (not including the orthosteric binding site) have been classified as “secondary binding pockets (SBP)”.68−70 Herein, we describe all other binding sites than the orthosteric as allosteric binding sites. Moreover, a bitopic ligand is defined as a ligand that concomitantly binds to the orthosteric and an allosteric binding site in a single receptor monomer.

binding of 10 at one protomer of a D2R dimer allosterically modulates the binding of dopamine at the other protomer. Through fragmentation of 10 closely related and more potent bitopic antagonists were discovered, but no data on selectivity was reported.64,65 In 2009, Newman et al. reported 11 (R-22) to be a 400-fold selective bitopic ligand for D3R (over D2R).66 In a molecular modeling study, 11 was docked into the D3R crystal structure and was found to bind in a bitopic mode.68 Several studies were performed to determine what induced the marked selectivity for D3R, ranging from computational methods to mutagenesis studies and functional assays with 11 fragments.71−73 The size and shape of extracellular loop (ECL) 1 was found to be an important determinant for subtype selectivity. Unlike in D3R, the ECL1 in D2R is relatively short and thus cannot accommodate the allosteric pharmacophore of 11.71 Moreover, the short ECL1 for D2R influences the orientation of transmembrane helix 2 unfavorably.72 Vass et al. reported a fragment based docking approach on D3R, which resulted in the synthesis of a new type of bitopic ligands (e.g., 12).67 Subnanomolar affinities were achieved and some of the ligands showed moderate selectivity for D2R over D3R. Other GPCR Targets. The majority of bitopic GPCR ligands have been reported for muscarinic acetylcholine, adenosine, and dopamine receptors. However, one successful example of a bitopic ligand developed for another target, the sphingosine-1-phosphate receptor (S1PR), illustrates that the principle is applicable across the GPCR landscape.74 The S1PRs are implicated in vasculogenesis and immune cell regulation and development, and the nonselective S1PR modulator fingolimod is used in the treatment of multiple sclerosis.75 Two closely related bitopic ligands, 13 (SPM-242) and 14 (SPM-354) (Figure 6), have been reported for the S1PR subtype 3 (S1P3), for which the latter has improved potency and in vivo efficacy.74,76 Both ligands occupy the allosteric pocket while competing with sphingosine-1-phosphate for the orthosteric binding site. 4129

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

studies have also been reported for nonrhodopsin-like GPCRs, with comparable results. A metastable binding site was discovered for the agonist strychnine at the TAS2R46 bitter taste receptor in the vicinity of ECL1.89 The data is not directly comparable to that outlined above for the rhodopsinlike GPCRs, however, it does fortify the hypothesis that metastable binding sites act as selectivity filters for ligand recognition and binding across GPCR families. Metastable Binding Sites in Other Systems. Recognition pathways similar to those described for GPCRs have also been discovered in other protein ligand systems, thus suggesting that metastable binding sites are prevalent across large parts of the human proteome. Other receptors, transporters, and enzymes alike display a variety of low affinity metastable complexes with their respective ligands. The entry pathways of the inhibitor DADMe-immucillin-H at the enzyme purine nucleoside phosphorylase were simulated by MD.77 It was discovered that the ligand can adopt two separate metastable poses before reaching its final binding site. In a similar manner, trypsin binds its ligand benzamidine via two metastable states.18 At the enzyme Src kinase, the situation appears to be similar, as the ligands PP1 and dasatinib enter the active center via multiple “persistent intermediate conformations”.80 However, in this study, the binding positions do not seem to be well-defined. Other reports do not explicitly mention transient low affinity binding sites but do acknowledge the presence of certain points of interaction and recognition on or near the surface of proteins. Hasenhuetl et al. proposed a “selectivity filter” in the entry pathways of the dopamine and serotonin transporters based on binding kinetic experiments.82 Another computational study dealt with two separate enzyme ligand systems: the haloalkane dehalogenase LinB and its ligand 1,2-dichloroethane, as well as levansucrase and its ligands glucose and sucrose.90 These ligands appear to have defined entry pathways along their respective enzyme’s surfaces, thus suggesting that key interactions are made along the ligands’ trajectory.

Figure 7. Process of 15 approaching the orthosteric binding pocket of β2AR via several metastable binding sites.16 The timeline for the simulation is indicated in μs for each binding pose. (1 (red)) The ligand moves from bulk solvent and attaches to the receptor. (2 and 3 (green)) The ligand pauses at two metastable binding sites in the extracellular vestibule. (4 and 5 (blue)) The ligand at the orthosteric binding site. Pose 4 is a long-lived metastable binding site overlapping with the final orthosteric binding pose 5. Pose 5 matches that seen in the X-ray crystal structure of β2AR in complex with alprenolol (PDB: 3NYA). The figure was adapted from Dror et al.16

were recently described by Wittmann and Strasser.87 Similarly to that described above for the β2AR they reported that Glu160ECL2 served as a point of first ligand attachment for both possible pathways. Metastable binding sites for the adenosine A2A receptor have also been proposed by Sabbadin et al. for several antagonists as well as the orthosteric agonist adenosine.20,21 In every case, Glu169ECL2 was seen to form a persistent interaction with the ligand and Asn2536.55 was also suggested as a ligand anchoring point upon receptor entry. Computational studies on the mAChR family showed comparable results to those for the A2A receptor. Ligands paused near known allosteric binding sites in the extracellular vestibule when entering and exiting the orthosteric binding pocket.54,78,79,88 For the mAChR M3 subtype, the large inverse agonist tiotropium bound at a metastable site interacting with Phe221ECL2 and Leu225ECL2, which coincidentally correlates with a known allosteric site in the extracellular vestibule reported by Dror et al.88 The small M3 (partial) agonists acetylcholine and arecoline paused in two separate “clusters”. One cluster was located just above the orthosteric binding site maintaining most of the key orthosteric interactions. The second cluster was found to be in the extracellular vestibule with important interactions to Phe1242.60 and Tyr1272.63. Tiotropium binds to the M2 subtype of mAChR in a similar manner as seen for M3 with key interactions to Tyr177ECL2 and Phe181ECL2 (corresponding to Phe221ECL2 and Leu225ECL2 in M3). Supervised MD simulations on the P2Y12 receptor with antagonist ticagrelor also indicate the presence of metastable binding sites.19 For the ligand to enter the orthosteric binding pocket of the P2Y12 receptor, a stable salt bridge must dissociate. The ligand interaction with the metastable binding site may be the driving force necessary for this energetically unfavorable rearrangement to occur. Some computational



CONCLUSION AND FUTURE DIRECTIONS In recent years, the development of bitopic ligands has transcended the subfamily of mAChRs and moved on to new targets such as the adenosine and dopamine receptors. The ability of such ligands to induce signaling bias, as illustrated with the A1AR ligand 9, showcases one of many unique properties that ligands of this class can possess. Unprecedented subtype selectivity can also be achieved with bitopic ligands and recent advances with dopamine receptors show promising results for the future. With the discovery of new allosteric binding sites in spatial proximity to the orthosteric binding site, the potential for bitopic ligands will increase and the first therapeutically promising bivalent ligands may see the light of day. When the presence of metastable binding sites in the extracellular vestibule was first computationally determined by Dror et al., they suggested that these transient binding sites could serve as targets for allosteric modulators.16 The presence of allosteric binding sites in the extracellular vestibule is now well established, and the PAM 3-amino-5-chloro-N-cyclopropyl-4-methyl-6-(2-(4-methylpiperazin-1-yl)-2-oxoethoxy)thieno[2,3-b]pyridine-2-carboxamide (LY2119620, structure not disclosed) serves as a prominent example of such a ligand targeting the mAChR M2 subtype.91 Besides in silico studies, no experimental evidence for the existence of metastable binding 4130

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

human β2AR in complex with the neutral antagonist alprenolol and docked 15 into both the orthosteric binding pocket and the metastable binding site in the extracellular vestibule (Figure 8) to obtain poses for both ligands in agreement with that published by others. The two docked ligands (15) are positioned in such a way that they may in principle be connected by a linker either via the channel that the ligand normally traverses in MD simulations (left)16 or alternatively through the wider channel (right). The idea of developing homobivalent bitopic ligands targeting the orthosteric and a metastable (i.e., allosteric) binding site simultaneously is in principle easily probed by the design of ligands based on known orthosteric ligands (e.g., 15, Figure 8). Work in our laboratories and those of others are underway to explore this strategy on the adrenergic and adenosine receptors and will be reported in due course.



AUTHOR INFORMATION

Corresponding Author

Figure 8. 15 (dark blue) docked into: (i) the orthosteric binding site of β2AR (PDB: 3NYA) overlaid with the crystallized ligand alprenolol (cyan), and (ii) the metastable binding site (dark blue).16 Dotted yellow lines highlight the two channels via which one could imagine connecting the ligands to create a potential homobivalent bitopic ligand. In MD simulations, ligands normally traverse the narrow channel to the left.16 The ligands were docked using the Schrödinger software suite with the OPLS3 force field.

*Phone: +45 35 33 60 00. Fax: +45 35 33 60 01. E-mail: daniel. [email protected]. ORCID

Daniel Sejer Pedersen: 0000-0003-3926-7047 Author Contributions †

These authors contributed equally.

Notes

The authors declare no competing financial interest.

sites has been provided to date. However, the presence of metastable binding sites appears likely based on the independent observation in silico by many different research groups on a diverse range of receptors. Interestingly, if the orthosteric ligand for a receptor is known it follows that the same ligand should bind to the metastable binding sites, thus providing a template for medicinal chemists. One could imagine elaborating on an orthosteric ligand to obtain an allosteric ligand binding in the extracellular vestibule. Alternatively, metastable binding sites may provide an entry to homobivalent bitopic ligands that simultaneously target the orthosteric and a metastable (i.e., allosteric) binding site (Figure 1). In theory, homobivalent bitopic ligands combine the benefits of high affinity and selectivity by targeting both the highly conserved orthosteric binding site and an evolutionarily less conserved metastable binding site. To the best of our knowledge, no one has reported this ligand design strategy. It should be mentioned that Gmeiner and co-workers published the synthesis and pharmacological evaluation of homodimeric ligands targeting D2-like receptor dimers.92 They discovered three ligands with increased affinity and enhanced selectivity on the D3 subtype, displacing only 1 equiv of radioligand. Thus, they reasoned that they were targeting the orthosteric binding pocket and an undefined allosteric binding site within one receptor rather than bridging the two orthosteric pockets on the GPCR homodimer. However, when considering the very long linker employed in this study, it is unlikely a metastable binding site is engaged. More likely the second pharmacophore is interacting with the extracellular surface of either one of the receptors. Jörg et al. reached a similar conclusion when targeting the same receptor with another class of bivalent ligands.93 Thus, the bivalent D2R ligands reported in these two instances are probably not binding to a metastable binding site in the extracellular vestibule. On the basis of the original reports by Dror et al. and González et al., we utilized the X-ray crystal structure of the

Biographies Philipp Fronik received his B.Sc. degree in Molecular Biotechnology from the University of Applied Sciences “FH Campus Wien” in Vienna, Austria (2013). In 2016, he completed a double M.Sc. degree in Medicinal Chemistry at the University of Copenhagen, Denmark/ VU University Amsterdam, The Netherlands. Birgit I. Gaiser obtained her B.Sc. in Chemistry (2011) from the University of Stuttgart, Germany, and her M.Sc. in Medicinal Chemistry (2014) from the University of Copenhagen, Denmark, where she currently is undertaking her graduate studies under the supervision of Associate Professor Daniel Sejer Pedersen. Daniel Sejer Pedersen received his M.Sc. degree under the supervision of Associate Professor Otto Dahl in 1999 (University of Copenhagen), and his Ph.D. degree under the supervision of Professor Stuart Warren in 2009 (University of Cambridge). He worked for several years in the biotech industry and as a Postdoctoral research fellow both in Australia and Denmark. In 2010, he was appointed Associate Professor in medicinal chemistry at the Department of Drug Design and Pharmacology, University of Copenhagen.



ACKNOWLEDGMENTS The Lundbeck Foundation, The Carlsberg Foundation, and the Faculty of Health and Medical Sciences, University of Copenhagen, are thanked for financial support.



ABBREVIATIONS USED A1AR, adenosine A1 receptor; ABS, allosteric binding site; AR, adrenergic receptor; D2R, dopamine D2 receptor; D3R, dopamine D3 receptor; ECL, extracellular loop; GPCR, G protein-coupled receptor; mAChR, muscarinic acetylcholine receptor; MBS, metastable binding site; MD, molecular dynamics; NAM, negative allosteric modulator; OBS, orthosteric binding site; OPLS, Optimized Potentials for Liquid Simulations; PAM, positive allosteric modulator; S1PR, 4131

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

ligand recognition at the P2Y12 receptor in light of X-ray structural information. J. Comput.-Aided Mol. Des. 2015, 29, 737−756. (20) Sabbadin, D.; Ciancetta, A.; Moro, S. Bridging molecular docking to membrane molecular dynamics to investigate GPCR-ligand recognition: The human A2A adenosine receptor as a key study. J. Chem. Inf. Model. 2014, 54, 169−183. (21) Sabbadin, D.; Ciancetta, A.; Deganutti, G.; Cuzzolin, A.; Moro, S. Exploring the recognition pathway at the human A2A adenosine receptor of the endogenous agonist adenosine using supervised molecular dynamics simulations. MedChemComm 2015, 6, 1081− 1085. (22) Adam, D.; Kessler, U.; Kronke, M. Cross-linking of the p55 tumor necrosis factor receptor cytoplasmic domain by a dimeric ligand induces nuclear factor-kappa B and mediates cell death. J. Biol. Chem. 1995, 270, 17482−17487. (23) Ho, S. N.; Biggar, S. R.; Spencer, D. M.; Schreiber, S. L.; Crabtree, G. R. Dimeric ligands define a role for transcriptional activation domains in reinitiation. Nature 1996, 382, 822−826. (24) Gangloff, M.; Murali, A.; Xiong, J.; Arnot, C. J.; Weber, A. N.; Sandercock, A. M.; Robinson, C. V.; Sarisky, R.; Holzenburg, A.; Kao, C.; Gay, N. J. Structural insight into the mechanism of activation of the toll receptor by the dimeric ligand Spätzle. J. Biol. Chem. 2008, 283, 14629−14635. (25) Portoghese, P. S. Bivalent ligands and the message-address concept in the design of selective opioid receptor antagonists. Trends Pharmacol. Sci. 1989, 10, 230−235. (26) Hiller, C.; Kühhorn, J.; Gmeiner, P. Class A G-protein-coupled receptor (GPCR) dimers and bivalent ligands. J. Med. Chem. 2013, 56, 6542−6559. (27) Huang, G.; Nimczick, M.; Decker, M. Rational modification of the biological profile of GPCR ligands through combination with other biologically active moieties. Arch. Pharm. 2015, 348, 531−540. (28) Karellas, P.; McNaughton, M.; Baker, S. P.; Scammells, P. J. Synthesis of bivalent β2-adrenergic and adenosine A1 receptor ligands. J. Med. Chem. 2008, 51, 6128−6137. (29) Lane, J. R.; Sexton, P. M.; Christopoulos, A. Bridging the gap: bitopic ligands of G-protein-coupled receptors. Trends Pharmacol. Sci. 2013, 34, 59−66. (30) Aurelio, L.; Valant, C.; Flynn, B. L.; Sexton, P. M.; Christopoulos, A.; Scammells, P. J. Allosteric modulators of the adenosine A1 receptor: synthesis and pharmacological evaluation of 4substituted 2-amino-3-benzoylthiophenes. J. Med. Chem. 2009, 52, 4543−4547. (31) Chen, X.; Klöckner, J.; Holze, J.; Zimmermann, C.; Seemann, W. K.; Schrage, R.; Bock, A.; Mohr, K.; Tränkle, C.; Holzgrabe, U.; Decker, M. Rational design of partial agonists for the muscarinic M1 acetylcholine receptor. J. Med. Chem. 2015, 58, 560−576. (32) Antony, J.; Kellershohn, K.; Mohr-Andrä, M.; Kebig, A.; Prilla, S.; Muth, M.; Heller, E.; Disingrini, T.; Dallanoce, C.; Bertoni, S.; Schrobang, J.; Tränkle, C.; Kostenis, E.; Christopoulos, A.; Höltje, H. D.; Barocelli, E.; De Amici, M.; Holzgrabe, U.; Mohr, K. Dualsteric GPCR targeting: a novel route to binding and signaling pathway selectivity. FASEB J. 2009, 23, 442−450. (33) Kruse, A. C.; Kobilka, B. K.; Gautam, D.; Sexton, P. M.; Christopoulos, A.; Wess, J. Muscarinic acetylcholine receptors: novel opportunities for drug development. Nat. Rev. Drug Discovery 2014, 13, 549−560. (34) Roszkowski, A. P. An unusual type of sympathetic ganglionic stimulant. J. Pharmacol. Exp. Ther. 1961, 132, 156−170. (35) Birdsall, N. J.; Burgen, A. S.; Hulme, E. C.; Stockton, J. M.; Zigmond, M. J. The effect of McN-A-343 on muscarinic receptors in the cerebral cortex and heart. Br. J. Pharmacol. 1983, 78, 257−259. (36) Valant, C.; Gregory, K. J.; Hall, N. E.; Scammells, P. J.; Lew, M. J.; Sexton, P. M.; Christopoulos, A. A novel mechanism of G proteincoupled receptor functional selectivity. J. Biol. Chem. 2008, 283, 29312−29321. (37) Spalding, T. A.; Trotter, C.; Skjærbæk, N.; Messier, T. L.; Currier, E. A.; Burstein, E. S.; Li, D.; Hacksell, U.; Brann, M. R.

sphingosine-1-phosphate receptor; S1P3, sphingosine-1-phosphate receptor 3; SAR, structure−activity relationship; TM, transmembrane



REFERENCES

(1) Stevens, R. C.; Cherezov, V.; Katritch, V.; Abagyan, R.; Kuhn, P.; Rosen, H.; Wü thrich, K. The GPCR network: a large-scale collaboration to determine human GPCR structure and function. Nat. Rev. Drug Discovery 2013, 12, 25−34. (2) Fredriksson, R.; Lagerström, M. C.; Lundin, L. G.; Schiöth, H. B. The G-protein-coupled receptors in the human genome form five main families. Phylogenetic analysis, paralogon groups, and fingerprints. Mol. Pharmacol. 2003, 63, 1256−1272. (3) Katritch, V.; Cherezov, V.; Stevens, R. C. Structure-function of the G protein-coupled receptor superfamily. Annu. Rev. Pharmacol. Toxicol. 2013, 53, 531−556. (4) Lagerström, M. C.; Schiöth, H. B. Structural diversity of G protein-coupled receptors and significance for drug discovery. Nat. Rev. Drug Discovery 2008, 7, 339−357. (5) Overington, J. P.; Al-Lazikani, B.; Hopkins, A. L. How many drug targets are there? Nat. Rev. Drug Discovery 2006, 5, 993−996. (6) Wootten, D.; Christopoulos, A.; Sexton, P. M. Emerging paradigms in GPCR allostery: implications for drug discovery. Nat. Rev. Drug Discovery 2013, 12, 630−644. (7) Lane, J. R.; Abdul-Ridha, A.; Canals, M. Regulation of G proteincoupled receptors by allosteric ligands. ACS Chem. Neurosci. 2013, 4, 527−534. (8) Bradley, S. J.; Tobin, A. B. Design of next-generation G protein− coupled receptor drugs: linking novel pharmacology and in vivo animal models. Annu. Rev. Pharmacol. Toxicol. 2016, 56, 535−559. (9) Christopoulos, A. Advances in GPCR allostery: from function to structure. Mol. Pharmacol. 2014, 86, 463−478. (10) Feng, Z.; Hu, G.; Ma, S.; Xie, X. Q. Computational advances for the development of allosteric modulators and bitopic ligands in G protein-coupled receptors. AAPS J. 2015, 17, 1080−1095. (11) Gentry, P. R.; Sexton, P. M.; Christopoulos, A. Novel allosteric modulators of G protein-coupled receptors. J. Biol. Chem. 2015, 290, 19478−19488. (12) Mohr, K.; Schmitz, J.; Schrage, R.; Tränkle, C.; Holzgrabe, U. Molecular alliance-from orthosteric and allosteric ligands to dualsteric/ bitopic agonists at G−protein coupled receptors. Angew. Chem., Int. Ed. 2013, 52, 508−516. (13) Soriano-Ursúa, M. A.; Trujillo-Ferrara, J. G.; Arias-Montaño, J. A.; Villalobos-Molina, R. Insights into a defined secondary binding region on β-adrenoceptors and putative roles in ligand binding and drug design. MedChemComm 2015, 6, 991−1002. (14) Baltos, J. A.; Gregory, K. J.; White, P. J.; Sexton, P. M.; Christopoulos, A.; May, L. T. Quantification of adenosine A1 receptor biased agonism: Implications for drug discovery. Biochem. Pharmacol. 2016, 99, 101−112. (15) Keov, P.; Valant, C.; Devine, S. M.; Lane, J. R.; Scammells, P. J.; Sexton, P. M.; Christopoulos, A. Reverse engineering of the selective agonist TBPB unveils both orthosteric and allosteric modes of action at the M1 muscarinic acetylcholine receptor. Mol. Pharmacol. 2013, 84, 425−437. (16) Dror, R. O.; Pan, A. C.; Arlow, D. H.; Borhani, D. W.; Maragakis, P.; Shan, Y.; Xu, H.; Shaw, D. E. Pathway and mechanism of drug binding to G-protein-coupled receptors. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 13118−13123. (17) González, A.; Perez-Acle, T.; Pardo, L.; Deupi, X. Molecular basis of ligand dissociation in β-adrenergic receptors. PLoS One 2011, 6, e23815. (18) Buch, I.; Giorgino, T.; De Fabritiis, G. Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 10184−10189. (19) Paoletta, S.; Sabbadin, D.; von Kügelgen, I.; Hinz, S.; Katritch, V.; Hoffmann, K.; Abdelrahman, A.; Straßburger, J.; Baqi, Y.; Zhao, Q.; Stevens, R. C.; Moro, S.; Müller, C. E.; Jacobson, K. A. Modeling 4132

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

Discovery of an ectopic activation site on the M1muscarinic receptor. Mol. Pharmacol. 2002, 61, 1297−1302. (38) Langmead, C. J.; Austin, N. E.; Branch, C. L.; Brown, J. T.; Buchanan, K. A.; Davies, C. H.; Forbes, I. T.; Fry, V. A. H.; Hagan, J. J.; Herdon, H. J.; Jones, G. A.; Jeggo, R.; Kew, J. N. C.; Mazzali, A.; Melarange, R.; Patel, N.; Pardoe, J.; Randall, A. D.; Roberts, C.; Roopun, A.; Starr, K. R.; Teriakidis, A.; Wood, M. D.; Whittington, M.; Wu, Z.; Watson, J. Characterization of a CNS penetrant, selective M1 muscarinic receptor agonist, 77-LH-28−1. Br. J. Pharmacol. 2008, 154, 1104−1115. (39) Lebois, E. P.; Digby, G. J.; Sheffler, D. J.; Melancon, B. J.; Tarr, J. C.; Cho, H. P.; Miller, N. R.; Morrison, R.; Bridges, T. M.; Xiang, Z.; Scott Daniels, J.; Wood, M. R.; Conn, P. J.; Lindsley, C. W. Development of a highly selective, orally bioavailable and CNS penetrant M1 agonist derived from the MLPCN probe ML071. Bioorg. Med. Chem. Lett. 2011, 21, 6451−6455. (40) Lebois, E. P.; Bridges, T. M.; Lewis, L. M.; Dawson, E. S.; Kane, A. S.; Xiang, Z.; Jadhav, S. B.; Yin, H.; Kennedy, J. P.; Meiler, J.; Niswender, C. M.; Jones, C. K.; Conn, P. J.; Weaver, C. D.; Lindsley, C. W. Discovery and characterization of novel subtype-selective allosteric agonists for the investigation of M1 receptor function in the central nervous system. ACS Chem. Neurosci. 2010, 1, 104−121. (41) Jones, C. K.; Brady, A. E.; Davis, A. A.; Xiang, Z.; Bubser, M.; Tantawy, M. N.; Kane, A. S.; Bridges, T. M.; Kennedy, J. P.; Bradley, S. R.; Peterson, T. E.; Ansari, M. S.; Baldwin, R. M.; Kessler, R. M.; Deutch, A. Y.; Lah, J. J.; Levey, A. I.; Lindsley, C. W.; Conn, P. J. Novel selective allosteric activator of the M1 muscarinic acetylcholine receptor regulates amyloid processing and produces antipsychoticlike activity in rats. J. Neurosci. 2008, 28, 10422−10433. (42) Steinfeld, T.; Mammen, M.; Smith, J.; Wilson, R. D.; Jasper, J. R. A novel multivalent ligand that bridges the allosteric and orthosteric binding sites of the M2 muscarinic receptor. Mol. Pharmacol. 2007, 72, 291−302. (43) Valant, C.; May, L. T.; Aurelio, L.; Chuo, C. H.; White, P. J.; Baltos, J. A.; Sexton, P. M.; Scammells, P. J.; Christopoulos, A. Separation of on-target efficacy from adverse effects through rational design of a bitopic adenosine receptor agonist. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 4614−4619. (44) Avlani, V. A.; Langmead, C. J.; Guida, E.; Wood, M. D.; Tehan, B. G.; Herdon, H. J.; Watson, J. M.; Sexton, P. M.; Christopoulos, A. Orthosteric and allosteric modes of interaction of novel selective agonists of the M1 muscarinic acetylcholine receptor. Mol. Pharmacol. 2010, 78, 94−104. (45) Daval, S. B.; Valant, C.; Bonnet, D.; Kellenberger, E.; Hibert, M.; Galzi, J. L.; Ilien, B. Fluorescent derivatives of AC-42 to probe bitopic orthosteric/allosteric binding mechanisms on muscarinic M1 receptors. J. Med. Chem. 2012, 55, 2125−2143. (46) Daval, S. B.; Kellenberger, E.; Bonnet, D.; Utard, V.; Galzi, J. L.; Ilien, B. Exploration of the orthosteric/allosteric interface in human M1 muscarinic receptors by bitopic fluorescent ligands. Mol. Pharmacol. 2013, 84, 71−85. (47) Digby, G. J.; Utley, T. J.; Lamsal, A.; Sevel, C.; Sheffler, D. J.; Lebois, E. P.; Bridges, T. M.; Wood, M. R.; Niswender, C. M.; Lindsley, C. W.; Conn, P. J. Chemical modification of the M1 agonist VU0364572 reveals molecular switches in pharmacology and a bitopic binding mode. ACS Chem. Neurosci. 2012, 3, 1025−1036. (48) Melancon, B. J.; Gogliotti, R. D.; Tarr, J. C.; Saleh, S.; Chauder, B.; Lebois, E. P.; Cho, H. P.; Utley, T. J.; Sheffler, D. J.; Bridges, T. M.; Morrison, R. D.; Daniels, J. S.; Niswender, C. M.; Conn, P. J.; Lindsley, C. W.; Wood, M. R. Continued optimization of the MLPCN probe ML071 into highly potent agonists of the hM1 muscarinic acetylcholine receptor. Bioorg. Med. Chem. Lett. 2012, 22, 3467−3472. (49) Keov, P.; López, L.; Devine, S. M.; Valant, C.; Lane, J.; Scammells, P. J.; Sexton, P. M.; Christopoulos, A. Molecular mechanisms of bitopic ligand engagement with the M1 muscarinic acetylcholine receptor. J. Biol. Chem. 2014, 289, 23817−23837. (50) Gregory, K. J.; Hall, N. E.; Tobin, A. B.; Sexton, P. M.; Christopoulos, A. Identification of orthosteric and allosteric site

mutations in M2 muscarinic acetylcholine receptors that contribute to ligand-selective signaling bias. J. Biol. Chem. 2010, 285, 7459−7474. (51) Nawaratne, V.; Leach, K.; Felder, C. C.; Sexton, P. M.; Christopoulos, A. Structural determinants of allosteric agonism and modulation at the M4 muscarinic acetylcholine receptor. J. Biol. Chem. 2010, 285, 19012−19021. (52) Disingrini, T.; Muth, M.; Dallanoce, C.; Barocelli, E.; Bertoni, S.; Kellershohn, K.; Mohr, K.; De Amici, M.; Holzgrabe, U. Design, synthesis, and action of oxotremorine-related hybrid-type allosteric modulators of muscarinic acetylcholine receptors. J. Med. Chem. 2006, 49, 366−372. (53) Bonifazi, A.; Yano, H.; Del Bello, F.; Farande, A.; Quaglia, W.; Petrelli, R.; Matucci, R.; Nesi, M.; Vistoli, G.; Ferré, S.; Piergentili, A. Synthesis and biological evaluation of a novel series of heterobivalent muscarinic ligands based on xanomeline and 1-[3-(4-butylpiperidin-1yl)propyl]-1,2,3,4-tetrahydroquinolin-2-one (77-LH-28−1). J. Med. Chem. 2014, 57, 9065−9077. (54) Bock, A.; Merten, N.; Schrage, R.; Dallanoce, C.; Bätz, J.; Klöckner, J.; Schmitz, J.; Matera, C.; Simon, K.; Kebig, A.; Peters, L.; Müller, A.; Schrobang-Ley, J.; Tränkle, C.; Hoffmann, C.; De Amici, M.; Holzgrabe, U.; Kostenis, E.; Mohr, K. The allosteric vestibule of a seven transmembrane helical receptor controls G-protein coupling. Nat. Commun. 2012, 3, 1044. (55) Schmitz, J.; van der Mey, D.; Bermudez, M.; Klöckner, J.; Schrage, R.; Kostenis, E.; Tränkle, C.; Wolber, G.; Mohr, K.; Holzgrabe, U. Dualsteric muscarinic antagonists–orthosteric binding pose controls allosteric subtype selectivity. J. Med. Chem. 2014, 57, 6739−6750. (56) Romagnoli, R.; Baraldi, P. G.; IJzerman, A. P.; Massink, A.; Cruz-Lopez, O.; Lopez-Cara, L. C.; Saponaro, G.; Preti, D.; Aghazadeh Tabrizi, M.; Baraldi, S.; Moorman, A. R.; Vincenzi, F.; Borea, P. A.; Varani, K. Synthesis and biological evaluation of novel allosteric enhancers of the A1 adenosine receptor based on 2-amino-3-(4′chlorobenzoyl)-4-substituted-5-arylethynyl thiophene. J. Med. Chem. 2014, 57, 7673−7686. (57) Narlawar, R.; Lane, J. R.; Doddareddy, M.; Lin, J.; Brussee, J.; IJzerman, A. P. Hybrid ortho/allosteric ligands for the adenosine A1 receptor. J. Med. Chem. 2010, 53, 3028−3037. (58) Beaulieu, J. M.; Gainetdinov, R. R. The physiology, signaling, and pharmacology of dopamine receptors. Pharmacol. Rev. 2011, 63, 182−217. (59) Heidbreder, C. A.; Gardner, E. L.; Xi, Z. X.; Thanos, P. K.; Mugnaini, M.; Hagan, J. J.; Ashby, C. R. The role of central dopamine D3 receptors in drug addiction: a review of pharmacological evidence. Brain Res. Rev. 2005, 49, 77−105. (60) Joyce, J. N.; Millan, M. J. Dopamine D3 receptor antagonists as therapeutic agents. Drug Discovery Today 2005, 10, 917−925. (61) Newman, A. H.; Grundt, P.; Nader, M. A. Dopamine D3 receptor partial agonist and antagonist as potential drug abuse therapeutic agents. J. Med. Chem. 2005, 48, 3663−3679. (62) Lane, J. R.; Donthamsetti, P.; Shonberg, J.; Draper-Joyce, C. J.; Dentry, S.; Michino, M.; Shi, L.; López, L.; Scammells, P. J.; Capuano, B.; Sexton, P. M.; Javitch, J. A.; Christopoulos, A. A new mechanism of allostery in a G protein−coupled receptor dimer. Nat. Chem. Biol. 2014, 10, 745−752. (63) Stemp, G.; Ashmeade, T.; Branch, C. L.; Hadley, M. S.; Hunter, A. J.; Johnson, C. N.; Nash, D. J.; Thewlis, K. M.; Vong, A. K.; Austin, N. E.; Jeffrey, P.; Avenell, K. Y.; Boyfield, I.; Hagan, J. J.; Middlemiss, D. N.; Reavill, C.; Riley, G. J.; Routledge, C.; Wood, M. Design and synthesis of trans-N-[4-[2-(6-cyano-1,2,3, 4-tetrahydroisoquinolin-2yl)ethyl]cyclohexyl]-4-quinolinecarboxamide (SB-277011): A potent and selective dopamine D3 receptor antagonist with high oral bioavailability and CNS penetration in the rat. J. Med. Chem. 2000, 43, 1878−1885. (64) Mistry, S. N.; Shonberg, J.; Draper-Joyce, C. J.; Klein Herenbrink, C.; Michino, M.; Shi, L.; Christopoulos, A.; Capuano, B.; Scammells, P. J.; Lane, J. R. Discovery of a novel class of negative allosteric modulator of the dopamine D 2 receptor through fragmentation of a bitopic ligand. J. Med. Chem. 2015, 58, 6819−6843. 4133

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134

Journal of Medicinal Chemistry

Perspective

(80) Shan, Y.; Kim, E. T.; Eastwood, M. P.; Dror, R. O.; Seeliger, M. A.; Shaw, D. E. How does a drug molecule find its target binding site ? J. Am. Chem. Soc. 2011, 133, 9181−9183. (81) Spiliotopoulos, D.; Caflisch, A. Molecular dynamics simulations of bromodomains reveal binding-site flexibility and multiple binding modes of the natural ligand acetyl-lysine. Isr. J. Chem. 2014, 54, 1084− 1092. (82) Hasenhuetl, P. S.; Schicker, K.; Koenig, X.; Li, Y.; Sarker, S.; Stockner, T.; Sucic, S.; Sitte, H. H.; Freissmuth, M.; Sandtner, W. Ligand selectivity among the dopamine and serotonin transporters specified by the forward binding reaction. Mol. Pharmacol. 2015, 88, 12−18. (83) Kass, I.; Buckle, A. M.; Borg, N. Understanding the structural dynamics of TCR-pMHC complex interactions. Trends Immunol. 2014, 35, 604−612. (84) Moro, S.; Hoffmann, C.; Jacobson, K. Role of the extracellular loops of G protein-coupled receptors in ligand recognition: A molecular modeling study of the human P2Y1 receptor. Biochemistry 1999, 38, 3498−3507. (85) Cooke, R. M.; Brown, A. J. H.; Marshall, F. H.; Mason, J. S. Structures of G protein-coupled receptors reveal new opportunities for drug discovery. Drug Discovery Today 2015, 20, 1355−1364. (86) Rasmussen, S. G. F.; Choi, H. J.; Rosenbaum, D. M.; Kobilka, T. S.; Thian, F. S.; Edwards, P. C.; Burghammer, M.; Ratnala, V. R. P.; Sanishvili, R.; Fischetti, R. F.; Schertler, G. F. X.; Weis, W. I.; Kobilka, B. K. Crystal structure of the human β2 adrenergic G-protein-coupled receptor. Nature 2007, 450, 383−387. (87) Wittmann, H. J.; Strasser, A. Binding pathway of histamine to the hH4R, observed by unconstrained molecular dynamics. Bioorg. Med. Chem. Lett. 2015, 25, 1259−1268. (88) Dror, R. O.; Green, H. F.; Valant, C.; Borhani, D. W.; Valcourt, J. R.; Pan, A. C.; Arlow, D. H.; Canals, M.; Lane, J. R.; Rahmani, R.; Baell, J. B.; Sexton, P. M.; Christopoulos, A.; Shaw, D. E. Structural basis for modulation of a G-protein-coupled receptor by allosteric drugs. Nature 2013, 503, 295−299. (89) Sandal, M.; Behrens, M.; Brockhoff, A.; Musiani, F.; Giorgetti, A.; Carloni, P.; Meyerhof, W. Evidence for a transient additional ligand binding site in the TAS2R46 bitter taste receptor. J. Chem. Theory Comput. 2015, 11, 4439−4449. (90) Negami, T.; Shimizu, K.; Terada, T. Coarse-grained molecular dynamics simulations of protein-ligand binding. J. Comput. Chem. 2014, 35, 1835−1845. (91) Kruse, A. C.; Ring, A. M.; Manglik, A.; Hu, J.; Hu, K.; Eitel, K.; Hubner, H.; Pardon, E.; Valant, C.; Sexton, P. M.; Christopoulos, A.; Felder, C. C.; Gmeiner, P.; Steyaert, J.; Weis, W. I.; Garcia, K. C.; Wess, J.; Kobilka, B. K. Activation and allosteric modulation of a muscarinic acetylcholine receptor. Nature 2013, 504, 101−106. (92) Huber, D.; Löber, S.; Hübner, H.; Gmeiner, P. Bivalent molecular probes for dopamine D2-like receptors. Bioorg. Med. Chem. 2012, 20, 455−466. (93) Jörg, M.; Kaczor, A. A.; Mak, F. S.; Lee, K. C.; Poso, A.; Miller, N. D.; Scammells, P. J.; Capuano, B. Investigation of novel ropinirole analogues: synthesis, pharmacological evaluation and computational analysis of dopamine D2 receptor functionalized congeners and homobivalent ligands. MedChemComm 2014, 5, 891.

(65) Shonberg, J.; Draper-Joyce, C.; Mistry, S. N.; Christopoulos, A.; Scammells, P. J.; Lane, J. R.; Capuano, B. A structure-activity study of N-((trans)-4-(2-(7-cyano-3,4-dihydroisoquinolin-2(1 H)-yl)ethyl)cyclohexyl)-1H-indole-2-carboxamide (SB269652); A bitopic ligand that acts as a negative allosteric modulator of the dopamine D2 receptor. J. Med. Chem. 2015, 58, 5287−5307. (66) Newman, A. H.; Grundt, P.; Cyriac, G.; Deschamps, J. R.; Taylor, M.; Kumar, R.; Ho, D.; Luedtke, R. R. N-(4-(4-(2,3-dichloroor 2-methoxyphenyl)piperazin-1-yl)butyl) heterobiarylcarboxamides with functionalized linking chains as high affinity and enantioselective D3 receptor antagonists. J. Med. Chem. 2009, 52, 2559−2570. (67) Vass, M.; Agai-Csongor, E.; Horti, F.; Keserü, G. M. Multiple fragment docking and linking in primary and secondary pockets of dopamine receptors. ACS Med. Chem. Lett. 2014, 5, 1010−1014. (68) Chien, E. Y. T.; Liu, W.; Zhao, Q.; Katritch, V.; Won Han, G.; Hanson, M. A.; Shi, L.; Newman, A. H.; Javitch, J. A.; Cherezov, V.; Stevens, R. C. Structure of the human dopamine D3 receptor in complex with a D2/D3 selective antagonist. Science 2010, 330, 1091. (69) Michino, M.; Donthamsetti, P.; Beuming, T.; Banala, A.; Duan, L.; Roux, T.; Han, Y.; Trinquet, E.; Newman, A. H.; Javitch, J. A.; Shi, L. A single glycine in extracellular loop 1 is the critical determinant for pharmacological specificity of dopamine D2 and D3 receptors. Mol. Pharmacol. 2013, 84, 854. (70) Furman, C. A.; Roof, R. A.; Moritz, A. E.; Miller, B. N.; Doyle, T. B.; Free, R. B.; Banala, A. K.; Paul, N. M.; Kumar, V.; Sibley, C. D.; Newman, A. H.; Sibley, D. R. Investigation of the binding and functional properties of extended length D3 dopamine receptorselective antagonists. Eur. Neuropsychopharmacol. 2015, 25, 1448− 1461. (71) Cao, Y.; Min, C.; Acharya, S.; Kim, K. M.; Cheon, S. H. Design, synthesis and evaluation of bitopic arylpiperazinephenyl-1,2,4oxadiazoles as preferential dopamine D3 receptor ligands. Bioorg. Med. Chem. 2016, 24, 191−200. (72) Michino, M.; Donthamsetti, P.; Beuming, T.; Banala, A.; Duan, L.; Roux, T.; Han, Y.; Trinquet, E.; Newman, A. H.; Javitch, J.; Shi, L. A single glycine in extracellular loop 1 is the critical determinant for pharmacological specificity of dopamine D2 and D3 receptors. Mol. Pharmacol. 2013, 84, 854−864. (73) Newman, A. H.; Beuming, T.; Banala, A. K.; Donthamsetti, P.; Pongetti, K.; Labounty, A.; Levy, B.; Cao, J.; Michino, M.; Luedtke, R. R.; Javitch, J.; Shi, L. Molecular determinants of selectivity and efficacy at the dopamine D3 receptor. J. Med. Chem. 2012, 55, 6689−6699. (74) Sanna, M. G.; Vincent, K. P.; Repetto, E.; Nguyen, N.; Brown, S. J.; Abgaryan, L.; Riley, S. W.; Leaf, N. B.; Cahalan, S. M.; Kiosses, W. B.; Kohno, Y.; Brown, J. H.; McCulloch, A. D.; Rosen, H.; GonzalezCabrera, P. J. Bitopic S1P3 antagonist rescue from complete heart block: pharmacological and genetic evidence for direct S1P3 regulation of mouse cardiac conduction. Mol. Pharmacol. 2016, 89, 176−186. (75) Brinkmann, V.; Billich, A.; Baumruker, T.; Heining, P.; Schmouder, R.; Francis, G.; Aradhye, S.; Burtin, P. Fingolimod (FTY720): discovery and development of an oral drug to treat multiple sclerosis. Nat. Rev. Drug Discovery 2010, 9, 883−897. (76) Jo, E.; Bhhatarai, B.; Repetto, E.; Guerrero, M.; Riley, S.; Brown, S. J.; Kohno, Y.; Roberts, E.; Schürer, S. C.; Rosen, H. Novel selective allosteric and bitopic ligands for the S1P3 receptor. ACS Chem. Biol. 2012, 7, 1975−1983. (77) Decherchi, S.; Berteotti, A.; Bottegoni, G.; Rocchia, W.; Cavalli, A. The ligand binding mechanism to purine nucleoside phosphorylase elucidated via molecular dynamics and machine learning. Nat. Commun. 2015, 6, 6155. (78) Kappel, K.; Miao, Y.; McCammon, J. A. Accelerated molecular dynamics simulations of ligand binding to a muscarinic G-proteincoupled receptor. Q. Rev. Biophys. 2015, 48, 479−487. (79) Kruse, A. C.; Hu, J.; Pan, A. C.; Arlow, D. H.; Rosenbaum, D. M.; Rosemond, E.; Green, H. F.; Liu, T.; Chae, P. S.; Dror, R. O.; Shaw, D. E.; Weis, W. I.; Wess, J.; Kobilka, B. K. Structure and dynamics of the M3 muscarinic acetylcholine receptor. Nature 2012, 482, 552−556. 4134

DOI: 10.1021/acs.jmedchem.6b01601 J. Med. Chem. 2017, 60, 4126−4134