The Impact of the Secondary Binding Pocket on the Pharmacology of Class A GPCRs

Abstract

G-protein coupled receptors (GPCRs) are considered important therapeutic targets due to their pathophysiological significance and pharmacological relevance. Class A receptors represent the largest group of GPCRs that gives the highest number of validated drug targets. Endogenous ligands bind to the orthosteric binding pocket (OBP) embedded in the intrahelical space of the receptor. During the last 10 years, however, it has been turned out that in many receptors there is secondary binding pocket (SBP) located in the extracellular vestibule that is much less conserved. In some cases, it serves as a stable allosteric site harbouring allosteric ligands that modulate the pharmacology of orthosteric binders. In other cases it is used by bitopic compounds occupying both the OBP and SBP. In these terms, SBP binding moieties might influence the pharmacology of the bitopic ligands. Together with others, our research group showed that SBP binders contribute significantly to the affinity, selectivity, functional activity, functional selectivity and binding kinetics of bitopic ligands. Based on these observations we developed a structure-based protocol for designing bitopic compounds with desired pharmacological profile.

Introduction

G-protein coupled receptors (Figure 1) are among the most popular targets for drug discovery and the development of novel therapeutic and pharmacological tools. One third of the drugs currently approved by the Food and Drug Administration affects one of the GPCRs (Sriram and Insel, 2018). They are critical in signal transduction of hormones and neurotransmitters, and consequently are pharmacological targets for many diseases (Overington et al., 2006). Furthermore, studying these receptors may help to elucidate the signaling mechanisms in cells, as they play a crucial role in the regulation of both central and peripherial neurological and physiological processes. Detailed understanding of these processes facilitates the development of more targeted therapies (Christopoulos, 2014).

FIGURE 1

GPCRs have multiple ligand binding sites, the orthosteric binding pocket and a generally separated less conserved allosteric secondary binding pocket (Christopoulos, 2014). Basically, the endogenous ligand binds to the OBP. SBPs are found in both the extracellular and intracellular parts of the receptor (Figure 2), some of these binding sites are well separated from the OBP while others may have extended binding pocket-like features such as the 5-HT_2A aripirazole structure (PDB: 7VOE) (Chen et al., 2021).

FIGURE 2

These secondary binding sites have become key to achieve the right subtype selectivity and functionality. Therefore, a lot of effort was given to the research of allosteric binding sites and allosteric modulators. A large number of allosteric modulators of GPCRs that bind to the extracellular or intracellular domains were identified. The combination of a primary pharmacophore (PP) binding to the OBP and a secondary pharmacophore (SP) binding to the SBP resulted in bitopic compounds (Figure 2C) that combine the pharmacological properties of both types of ligands defining a new unique pharmacological profile. One of the first published bitopic molecules of this type is methoctramine that acts as an antagonist at the muscarinic receptor M₂R (Melchiorre et al., 1987).

In this review we would like to give only a brief insight into class A GPCR structures and the world of allosteric modulators as several reviews have been published in the field. Mainly, we discuss in detail the recent advances in bitopic ligands, while we close the review with an outlook towards the design approaches in the field.

Ligand Binding Pocket Revealed by Experimental Structures

Recent advances in X-ray crystallography and cryo electron microscopy provided many new structures of GPCRs complexed with allosteric ligands. As of early December 2021, 57 GPCR structures containing allosteric ligands have been found in GPCRdb (Kooistra et al., 2021), these structures cover 20 receptor types and three different states; active, inactive and intermediate. Among allosteric ligands, examples of positive (PAM) and negative allosteric modulators (NAM) can be found. The collection of the published GPCR structures with allosteric ligands is available in the supporting information (Kooistra et al., 2021) (Supplementary Table S1). In addition, a significant number of active structures have become accessible, which may provide more information on the mechanism of receptor activation and offer considerable support for drug design, although few of these are allosteric ligands. Among them, 35 active aminergic GPCR structures have been published in the last 2 years (Supplementary Table S2) (Kooistra et al., 2021). These include 7 serotonin (Kim et al., 2020; Peiyu Xu et al., 2021a; Huang et al., 2021) (5-HTR), 15 dopamine (Zhuang et al., 2021a; Xiao et al., 2021; Zhuang et al., 2021b; Yin et al., 2020; Peiyu Xu et al., 2021b) (DR), 1 histamine (Xia et al., 2021) (HR), 1 muscarinic (Staus et al., 2020) (MR) and 11 adrenergic (Lee et al., 2020; Fan Yang et al., 2021; Yuan et al., 2020; Su et al., 2020; Xinyu Xu et al., 2021; Zhang et al., 2020; Nagiri et al., 2021) (AR) receptor structures. Out of these complexes, 20 structures contain allosteric modulators but not obviously in the SBP, while 10 were co-crystallized with bitopic ligands bound both the OBP and the SBP. The discussion of the structures in detail is out of scope of this review, however we highlight here the new cariprazine and aripiprazole bound 5-HT_2A structures (Figure 3A). (Chen et al., 2021) Interestingly, both compounds display an unexpected binding mode with their secondary binding motif exploring a binding pocket deep in the receptor instead of engaging with the extracellular secondary binding pocket. In the dopamine D₂ and D₃ receptors (D₂R, D₃R) the docking positions of aripiprazole so far have shown that 4-(2,3-dichlorophenyl)piperazine PP is located roughly parallel to the membrane plane and close to S5.42 and F6.51. The dihydroquinoline secondary pharmacophore is located at the junction of transmembrane helices (TM) 1, TM2, TM7 or TM3, TM5 and extracellular loop (ECL) 2. However, in the 5-HT_2A crystal structures of aripiprazole and cariprazine the ligands are located in an “upside-down” binding mode. The 2,3-dichlorophenyl PP occupies the orthosteric site and faces the extracellular region, but the dihydroquinoline SP vertically penetrates the hydrophobic pocket formed between TM5 and TM6 and interacts with residues L247^5.51, V333^6.45 and C337^6.49 and forms π-π interactions with residues F332^6.44 and W338^6.48. Upon binding of aripiprazole, a conformational rearrangement occurs resulting in an increase in the size of the binding pocket. Induced docking with D₂R was used to reproduce the “upside-down” binding pose of aripiprazole and cariprazine. Compared with the rigid docking, a much lower binding energy was calculated in the induced-fit docking, indicating that the upside-down binding mode represents a more stable conformation of D₂R (Chen et al., 2021).

FIGURE 3

Egyed et al., 2020). (C) The binding mode of salbutamol (cyan) and salmeterol (green) (Masureel et al., 2018) in the β₂R highlighting the important role of ECL2 as discussed in more detail in the binding kinetic section of this review.

Allosteric Modulators in the Class A GPCR Field

Allosteric binding sites (Figures 2A,B) have attracted increasing interest in order to develop more selective agents with fewer side effects (Congreve et al., 2017a; Chan et al., 2019). Allosteric sites are typically less conserved than orthosteric pockets and therefore they could provide greater selectivity and better control over the dynamical equilibrium of the receptor. Following the classic structural architecture of a class A GPCR, the orthosteric binding pocket is formed by the transmembrane helixes while the extracellular loops and the N-terminus of the peptide chain define the secondary binding domain. It should be mentioned, however, that there are other allosteric sites (e.g., extrahelical sites at the protein-membrane interface, intracellular sites at the signalling domain or intrahelical sodium site) available. Allosteric ligands can modify the biological response, they can stabilise the active or inactive conformation that is potentially linked to biased signalling or partial agonism (Wakefield et al., 2019). Based on spectroscopic and structural studies, conformational changes in the receptor govern the activation of signalling pathways. Characterization of interactions with intracellular partners guiding the allosteric process is a major challenge and can only be fully understood by using a combination of different methodologies (Liu et al., 2012; Masureel et al., 2018; Frei et al., 2020). Most allosteric modulators have been discovered serendipitously by high throughput screening (HTS) campaigns (Bian et al., 2020). Due to the vastness of the topic and the number of reviews published in the last years, we will only provide a brief insight into the world of allosteric modulators.

The tissue distribution and relative expression of the four adenosine receptor (AR) subtypes A₁R, A_2AR, A_2BR and A₃R regulate the physiological effects of endogenous adenosine. Adenosine receptors are expressed in most tissues and major organs, including brain, heart, kidney, skin, adipose tissue, immune cells, lung and liver. The four adenosine receptor subtypes can be broadly classified into two classes. Baressi et al. described a type of A_2BR allosteric modulators with good selectivity over the other subtypes, these compounds contain a 1,3-substituted indole unit (Barresi et al., 2021a; Barresi et al., 2021b). Lu et al. established a fragment screening method using mass spectrometry to screen GPCR ligands, identifying an A_2AR NAM. Fg754 (Figure 4) contains a specific acetidine moiety that forms bonds in the sodium ion pocket. Based on molecular dynamics (MD) simulations, it may overlap with the orthosteric binding site, probably acting in a mixed mode. The compound could thus be a new starting point for the development of allosteric modulators or bitopic compounds (Yan Lu et al., 2021). The A₁R and A₃R preferentially bind to G_i/o proteins to inhibit adenylate cyclase activity, while the A_2AR and A_2BR preferentially bind to G_s proteins to stimulate adenylate cyclase activity. Like other GPCRs, adenosine receptors can interact with different G-protein subtypes. In addition, A_2BR has been suggested to couple to both G_i/o and G_q proteins (Linden et al., 1999; Gao et al., 2018), while A₁R has been shown to couple to both G_s and G_q proteins (Cordeaux et al., 2004). In addition to G_α signalling, G_βγ dimers released following G protein activation can interact with effector proteins to modulate intracellular signalling. Beside G-protein-dependent signalling, adenosine receptors can also signal through G-protein-independent effectors. One of the best described G-protein-independent pathway is initiated following recruitment of arrestin adaptor proteins (β-arrestin1 and β-arrestin2). This process is typically preceded by G-protein coupled receptor kinase mediated phosphorylation, but recent studies have shown the possibility of phosphorylation independent β-arrestin recruitment for several GPCRs, including A₃R. Arrestin recruitment has been investigated primarily in A_2BR and A₃R and there is limited evidence that A₁R or A_2AR can recruit β-arrestin. McNeill et al. have discussed in detail the effects of allosteric modulators belonging to different subtypes on distorted signalling, which will not be discussed in detail below (McNeill et al., 2021).

FIGURE 4

Free fatty acids may act as signalling molecules at FFA receptors (FFARs). Free fatty acids of different chain lengths and saturation states activate FFARs as endogenous agonists by binding at the orthosteric receptor site. Following FFAR deorphanisation, a number of ligands targeting allosteric sites on FFARs have been identified with the aim of developing drugs for metabolic, (auto)inflammatory, infectious, endocrine, cardiovascular and renal diseases. In 2021, Grundmann et al. published a detailed review (Grundmann et al., 2021) on free fatty acid receptors, describing in detail the subtypes (FFAR1, FFAR2, FFAR3, FFAR4, GPR84), their function, structures and outlined the importance and challenges of allosteric modulators. FFAR1 is the most studied subtype. Although the biology of the receptors is still largely elusive, a large body of research evidence has accumulated around ligand-receptor interactions and their associated signalling capabilities. At least three distinct groups of FFAR1-activating ligands can be distinguished: 1) endogenous/orthosteric agonists (long-chain fatty acids), partial allosteric agonists (fasiglifam, MK-8666, AM 837), and full allosteric agonists (AM 1638, AP8) (Figure 2B, Figure 4). These groups differ not only in their apparent binding sites (Figure 2B) on the receptor, but also in their ability to induce different downstream signalling pathways of FFAR1, ultimately leading to different results in the phenotype of the FFA1 receptor in vivo. New results on allosteric FFAR2 ligands (AMG 7703, AZ1729, Compound 58) (Figure 4), show promising pharmacological properties and have generated new interest in this target, considering new allosteric modalities. GLPG1205 (Figure 4), an antagonist and negative allosteric modulator of GPR84, showed promising preclinical results in models of idiopathic pulmonary fibrosis, but was later discontinued from development. Allosteric targeting of small-, medium-, and long-chain fatty acid receptors is a promising approach to address a variety of therapeutic areas, demonstrating the biological diversity and drug target attractiveness of members of this receptor family (Grundmann et al., 2021).

The cannabinoid receptor type 1 (CB₁) was first discovered as the main target for Δ9-tetrahydrocannabinol (THC), the psychoactive compound in Cannabis. CB₁ was first identified in rat and later cloned from a human brain cDNA library. Widely known CB₁ agonists are synthetic cannabinoids and THC analogues, such as HU-210 (Howlett et al., 1990), CP55940 (Kapur et al., 2009), and WIN55212 (Felder et al., 1995). The CB₁ receptor preferentially binds a G_i protein and its activation leads to a decrease in cyclic adenosine monophosphate (cAMP) levels in cells. Other signalling pathways have also been investigated, focusing primarily on ERK1/2 phosphorylation. ERK signalling is hypothesised to play a role in cocaine addiction, and together with cAMP, to be an important regulator of synaptic plasticity, memory and learning. Inhibition of CB₁ proved effective in the treatment of obesity with antagonists or inverse agonists, but they were later withdrawn from the market due to adverse psychiatric side effects (anxiety, suicidal ideation). Several new strategies to avoid potential side effects have been analysed, one of them being the development of allosteric modulators. Leo and Abood reviewed the physiological and pathophysiological roles of CB₁, described the signalling mechanisms, and investigated CB₁ biased signaling (Leo and Abood, 2021). Based on agonist-bound solvated molecular structures and biased allosteric modulators they look at possible molecular mechanisms of CB₁ signalling. Mielnk et al. present the in vitro and in vivo profiles of several NAMs (Org27569, PSNCBAM-1, ABM300, Pepcan-12, Pregnenolone, and cannabidiol) and PAMs (ZCZ011, GAT211, Lipoxin A4, LDK1258) in detail (Figure 2B, Figure 4). They concluded that CB₁ PAMs in anxiety and depression while CB₁ NAMs—in combination with cannabidiol—in psychosis could be promising (Mielnik et al., 2021).

Che and Roth have provided a detailed summary of the pharmacology, ligands (orthosteric, allosteric), and structures of opioid receptors (OR) (Che and Roth, 2021). Activating µ-opioid receptor (MOR) causes serious side effects, which are the root of the current opioid crisis. In their review, potential strategies and targets for developing opioid alternatives were discussed. Separately, they list OR biased agonists, allosteric modulators, multitarget ligands and peripherally restricted ligands. The complexity of signalling pathways should be considered in the therapeutic potential of biased agonists, and allosteric modulators are alternative means to modulate more precisely the action of endogenous or exogenous ligands. As opioid receptors are widely expressed in the peripheral system, the use of ligands restricted to this system would avoid central nervous system induced side effects. Simultaneous targeting of multiple opioid and non-opioid receptors may result in safer analgesics (Che and Roth, 2021).

The family of aminergic GPCRs includes adrenergic, dopamine, serotonin, histamine, muscarinic and trace amine receptors. These receptors have several similarities, they bind monoamine neurotransmitters, acetylcholine, or trace amines. They share common features in sequence, structure and function. Ergotamine (Figure 3B) can bind to 22 aminergic receptors with K_i values less than 1 µM (Peng et al., 2018). Other examples can be found in the literature, such as chlorpromazine, clozapine, thioridazine, olanzapine which have good affinity for several aminergic GPCRs (Roth et al., 2004). On the other hand, it would be important to produce drugs that have subtype and functional selectivity to avoid side effects.

In the field of adrenergic receptors, Wu and co-workers have discussed in detail the binding of endogenous ligands to different receptors, the mechanism of β-adrenergic and α₂ receptor attenuation, distorted signal transduction, subtype selectivity, and selectivity between the main types. Insights into the allosteric modulation of β_2AR were provided. They also reported on the results obtained with different modalities. The cholesterol binding site was recently described in detail by Sarkar and Chattopadhyay (2020) The arrangement of the 7 TMs in each class of GPCRs results in a groove at the lipid interface formed by TM3/4/5, and in β_2AR, to this site the binding of PAMs and NAMs were identified. GPCRs use the cytoplasmic surface to interact with intracellular partners with small molecules binding at this site discovered primarily in chemokine receptors. Only Cmpd15PA (Figure 2B, Figure 4) in β_2AR targets this site outside the chemokine subfamily. These small molecules are all NAMs. Cmpd15PA has little interaction with the G protein, but stabilizes the receptor inactive state through extensive interactions with TM1, TM2, TM6, TM7, H8 and intracellular loop 1 (Wu et al., 2021).

The five dopamine receptor subtypes (D1–5) are activated by the endogenous catecholamine dopamine. The D1-like family comprises dopamine D₁ and D₅ receptors that mainly couple to the G_s G-protein and thereby stimulate cAMP production. The D₂-like family includes D₂, D₃, and D₄ receptors, that couple to G_i/o G-proteins and attenuate cAMP production (British Pharmacological Society, 2021). Fasciani et al. have presented allosteric modulators of the DR, the bitopic compound SB269652 has been analysed in detail. Mao et al. describe the role of different dopamine receptor allosteric modulators in the treatment of Parkinson’s disease. DR allosteric modulators represent an alternative and promising strategy for drug discovery of GPCRs with high selectivity and low side effects (Mao et al., 2020). Like many other receptors, the classical approach to D₁R is the development of orthosteric ligands, but this has several drawbacks from a therapeutic point of view. D₁R agonists have narrow therapeutic window, can induce seizures and hypotensive side effect. PAMs are a more useful approach because they potentiate the effect of endogenic dopamine, the available dopmaine level provides a natural ceiling effect for PAM activity, and endogenous spatial and temporal regulation of dopamine-mediated stimulation is maintained. To date, seven D₁R PAM structural classes have been discovered. Two of these (MLS1082 and MLS6585) were discovered in 2018 by Luderman and colleagues using HTS (Luderman et al., 2018) (Figure 4). Subsequently, MLS1082 was investigated in a SAR study and they identified several analogues that enhanced dopamine-induced D₁R activation (Luderman et al., 2021).

There are five subtypes of the muscarinic acetylcholine receptor. The different subtypes show high degree of homology in the transmembrane domains. In recent years, the structures of all five have been resolved by X-ray crystallography (Vuckovic et al., 2019; Thal et al., 2016; Kruse et al., 2013; Kruse et al., 2012; Haga et al., 2012). In a review, Jakunik and El-Fakahany provide a detailed analysis of allosteric adhesion, the molecular mechanisms of action, and present specific modulators. The diversity of the effects of allosteric modulators and the studies on them will greatly influence the development of new therapies. Selective PAMs (LY2119620) (Figure 2B, Figure 4), which have therapeutic potential in the treatment of Alzheimer’s disease or schizophrenia, show encouraging results (Conn et al., 2009; Bock et al., 2018; Jakubik and El-Fakahany, 2020).

Biochemically, 5-hydroxytryptamine (5-HT) is derived from the amino acid tryptophan, undergoing hydroxylation and decarboxylation processes that are catalyzed by tryptophan hydroxylase and aromatic L-amino acid decarboxylase, respectively. As a biogenic amine, 5-HT plays important roles in cardiovascular function, bowel motility, platelet aggregation, hormone release and psychiatric disorders. 5-HT achieves its physiological functions by targeting various 5-HT receptors (5-HTRs), which are composed of six classes (5-HT₁, 5-HT₂, 5-HT₄, 5-HT₅, 5-HT₆, and 5-HT₇ receptors, a total of 13 subtypes) and a class of cation-selective ligand-gated ion channels, the 5-HT₃ receptor. Barnes et al. have published a review (Barnes et al., 2021) detailing each subtype, describing their functions and pharmacology one by one and discuss known allosteric ligands. They find that 5-HT receptors are less involved in allosteric modulation than other GPCRs (e.g., muscarinic, GABA), with the possible exception of 5-HT₃R. However, from some structures with ergoline, it becomes clear that, in addition to the classical OBP, some 5-HT receptors have an extended binding site very similar to that described for muscarinic allosteric ligands. Such molecular targets may offer attractive strategies for new therapies (Barnes et al., 2021).

Bitopic Ligands to Study Selectivity and Functional Selectivity of Class A GPCRs

As outlined in the introduction, our primary focus is on bitopic compounds in this review. These compounds combine the efficiency of orthosteric ligands and the diversity of allosteric SPs by interacting with both binding sites simultaneously. This gives bitopic ligands an advantage over allosteric modulators, as the latter need an orthosteric ligand to exert their effect. This may be important in cases where endogenous substrate depletion contributes to the pathogenesis of disease, such as in Parkinson’s and Alzheimer’s diseases, but there are further examples in metabolic disorders. The key strucutural moieties of bitopic compounds (PP, SP and linker, depicted on Figure 3C) have different roles. PP is classically considered to be responsible for functionality while SP can modulate binding affinity, selectivity as well as functional character and efficacy. The linker connects the two pharmacophores and may be responsible for the optimal binding poses by positioning the pharmacophores and affecting the pharmacology profile (Bethany et al., 2019).

In the design of bitopic compounds, the desired orthosteric binding motif should have high affinity for the selected receptor and ideally, the SP should provide high subtype selectivity while maintaining or even increasing affinity. In the case of a linker, the choice of attachment points and length must be appropriate, and the linker must be moderately flexible to allow the pharmacophores to bind properly. For agonists, it is important that the linker does not interfere with conformational changes induced by receptor activation (Valant et al., 2012; Lane et al., 2013; Fronik et al., 2017; Bethany et al., 2019). Reinecke et al. published a review on bitopic compounds in 2019, summarizing the new bitopic compounds that have been published in the last 5 years (Bethany et al., 2019). Here we therefore focus on compounds published in 2020–21, with a contextual analysis of previously published compounds where appropriate. In the following subsections, we discuss subtype selectivity and functional selectivity results separately.

Receptor and Subtype Selectivity

Receptor and subtype selectivity is an important criterion for minimizing side effects, therefore tremendous efforts go into the development of compounds with designed binding profile.

Keserű et al. have developed a fragment based docking protocol to design specific receptor ligands. Based on the docking results, they have synthesized several compounds and demonstrated the usefulness of the method for the designing D₂/D₃, 5-HT_1B/5-HT_2B and H₁/M₁ receptor ligands with improved selectivity (Figure 5). In the first two cases, the selectivity of the PP was reversed using the SP moiety, while in the third case, a selective compound was designed and synthesized for a receptor pair with very similar PP (Egyed et al., 2021).

FIGURE 5

Egyed et al., 2021).

The importance of bitopic compounds in the inhibition of dopamine receptors is demonstrated by second and third generation antipsychotics, including aripiprazole (Burris et al., 2002) and cariprazine (Ágai-Csongor et al., 2012). 2,3-dichlorophenyl-piperazine, that serves as PP in these compounds was changed to 2-methoxyphenylpiperazine (1) PP. Although this PP exhibits weak D₂R selectivity, combined with a suitable SP group (2,3) its profile has been changed to mild D₃R selectivity. The efficacy of this methodology was further tested on serotonin receptors. The LSD-like PP of ergotamine (Figure 3B) did not show subtype selectivity at the two selected serotonin receptors, but the designed compounds (4,5) with the modified SP already had significantly higher affinity at 5-HT_2BR over 5-HT_1BR. Although ergotamine was more potent, compounds 4 and 5 had much greater selectivity over it. Among the first-generation antihistamines, muscarinic acetylcholine M₁ activity was a major problem due to side effects. Therefore, huge efforts were dedicated to the development of compounds with significant H₁R receptor selectivity. Starting from amitriptyline having only 7-fold selectivity, bitopic compounds (6,7) were designed that demonstrated 50–80 fold selectivity over M₁R (Egyed et al., 2021). The proposed protocol detailed in the design section of this review may be applied to other targets to achieve designed selectivity with bitopic compounds.

Tan et al. have exploited the basic 2-phenylcyclopropylmethylamine (PCPMA) scaffold (8, 9), whose analogues are known 5-HT_2CR agonists (Cheng et al., 2015; Cheng et al., 2016a; Cheng et al., 2016b; Zhang et al., 2017), to design new bitopic compounds (Tan et al., 2020) (Supplementary Table S3). Here we discuss only a subset of these compounds. As secondary pharmacophore, 1,2,4-triazolylthiol ethers were used and a propyl chain was employed as a linker. The introduction of SP alone improved D₃R activity 3-fold. A major leap forward was the realization that the alkyl side chain introduced on the amino group of PCPMA significantly improves subtype selectivity and D₃R affinity. Next, they investigated the substituents of the aromatic ring of PCPMA. First, the ortho positioned 2-fluoroethoxy group was changed, whereby methoxy was found to be the optimal one, thus significantly improving the D₃R affinity. The replacement of the fluorine atom by chlorine resulted in a moderate selectivity towards D₂R, D₄R, 5-HT_2CR and a strong selectivity towards D₁R and D₅R (10–12). As these results could only approximate the values of the reference compound 13 (BP-897) (Supplementary Table S3) the strategy was changed and a buthylene linker was used instead of the propylene group, the SP was replaced by other aromatic rings (naphthyl, indolyl, and 4-pyridylphenyl) and an amide bond between the linker and the SP was introduced instead of thioether (14–20). For these compounds, only N-alkyl substituted variants have been prepared and the effects of several PPs have been investigated. When examining the racemic compounds, the compound containing 4-pyridylphenyl SP and dichlorophenyl PP (20) has more than 1000-fold selectivity towards the other DRs, with milder but still significant selectivity in the range of 17, 18, and 19 (Tan et al., 2020) (Table 1).

TABLE 1

Cmpd	Structure	Ki (nM)
		D1R	D2R	D3R	D4R	D5R	5-HT2C
(1S,2S)-17a		1,071	1,230	3.8	851	>5,000	50.1
(1R,2R)-17b		4,898	1,349	4.1	575	>5,000	1,122
(1S,2S)-18a		1,047	1,148	20.8	776	>5,000	138
(1R,2R)-18b		1,288	676	4.4	813	>5,000	513
(1S,2S)-19a		1,122	992	12.8	676	>5,000	61.7
(1R,2R)-19b		1,380	537	2.2	1,047	>5,000	513
(1S,2S)-20a		2344	1,023	5.3	912	>5,000	44.7
(1R,2R)-20b		1,349	550	1.5	676	>5,000	417
Cmpd	Structure	D2R Ki(nM)	D3R Ki(nM)	D4R Ki(nM)	D2R/D3R	D4R/D3R
24		2600	24200	ND	0.110	ND
25		34.6	31.2	ND	1.1	ND
27		134	5.96	357	22.5	59.9
28a		87.8	1.85	286	47.5	155
28b		831	282	2930	2.95	10.4
39		648	1.4	-	467	-

Selected compounds from DR related selectivity studies (Battiti et al., 2019; Tan et al., 2020; Lee et al., 2021).

Battiti and co-workers performed a SAR analysis combining two PPs for the synthesis of bitopic compounds; one is a selective dopamine agonist PF-592379 (Allerton et al., 2005; Ackley, 2008) and the other is PD-128907, which is a D₂R/D₃R agonist. (Supplementary Table S4) They concluded that the structural features of PD-128907 avoided the construction of bitopic compound. Therefore, they focused to PF-592379 to synthesize D₂R/D₃R active bitopic compounds. Here we discuss a representative example for different SPs (Supplementary Table S4). D₂R and D₃R binding data clearly show that the (S,S) enantiomer of the PP is more favourable for receptor binding. The (S,S) enantiomer already plays a prominent role in PP (22, 23), with a 3-fold activity difference between the enantiomers. The same effect can be observed when using tetrahydroisoquinoline (24,25) or indole (26,27) SP, although here the difference in activity at the D₃ receptor is about 100-fold (Table 1, Supplementary Table S4). Compound 27 show a 22.5-fold subtype selectivity towards D₃R that is due to the SP moiety. Next the authors investigated the effect of the linker. Changing the original cyclopropylethyl linker (Supplementary Table S4) for the racemic derivative (rac-trans-28) resulted in 37.3-fold selectivity towards D₃R. Separating the enantiomers, (1S, 2R)-trans-cylopropyl stereochemistry (28a) showed D₃R K_i of 1.85 nM and an unprecedented 47.5-fold selectivity for D₃R over D₂R (D₂R Ki = 87.8 nM), while the other enantiomer (28b) has much weaker activity coupled with poor selectivity (Battiti et al., 2019). Finally, two additional linkers were used (31, 32) that are widely used among D₃R bitopic compounds including several high selectivity partial agonists or antagonists (Kumar et al., 2016; Michino et al., 2017; Verma et al., 2018). Compound 31 showed reduced affinity compared to 28a, inferring that the hydroxyl group on the linker is optimal for antagonism but cannot be directly transferred to the agonist binding mode due to different receptor conformations in the active and inactive states. Compound 32 shows good affinity but neither affinity nor selectivity reaches that of 28a. Compounds 31, 32, 28a were tested at MOR. For 32 there is a decrease in affinity at the dopamine receptor but the weak subtype selectivity is retained, however there is a 22.9-fold increase in activity at the MOR receptor (Supplementary Table S4) (Battiti et al., 2020). The same group synthesized a number of eticlopride analogues using different SPs in the 2-N or 4-C position of pyrrolidine via lycerol (Battiti et al., 2020). They found that O-alkylated analogues had better affinity for D₂ and D₃ receptors than the N-substituted derivatives. In BRET assays, these compounds exhibited antagonist or very weak partial agonist behaviour. Docking studies revealed that the SPs of the O-alkylated analogues form aromatic stacking interactions with conserved residues His6.55 and Tyr7.35 both in the D₂ and D₃ receptors, while the SPs of the N-alkylated derivatives extend towards the extracellular site that is less conserved (Shaik et al., 2021).

N-phenylpiperazine analogues were used extensively for constructing bitopic ligands against dopamine receptors. Lee et al. synthesized and evaluated a series of N-phenylpiperazine analogues substituted with 3-thiophen and 4-thiazolylphenylfluoride (Supplementary Table S5). They identified several ligands that bind with high affinity to D₃R and exhibit considerable selectivity towards D₂R. Comparison of the binding results of compounds 33–38 and 39–44 suggests that 39–44 binds to D₃R but not to D₂R. The replacement of the thiophene ring by a thiazole ring (45–50) led to a decrease in receptor binding selectivity. Compound 39 (Table 1) possessed the highest D₃R affinity (K_i = 1.4 nM) and 450-fold selectivity that nominated this compound for in vivo testing. Intraperitoneal administration of 39 led to a significant reduction in DOI-dependent head twitch response in mice and a reduction in AIM scores in dyskinetic hemiparkinsonian rats. These data suggest that compound 39 is able to cross the blood-brain barrier and achieves therapeutic concentrations (Lee et al., 2021).

Starting from the 5-HT_2A receptor-bound structure of aripiprazole and cariprazine Chen et al. designed D₂/D₃ receptor ligands with no significant 5-HT_2A affinity (Chen et al., 2021). The authors suggested that the unusal “upside-down” binding mode (Figure 3A) might affect the observed selectivity. According to the structural rearrangements, the location of the SP of aripiprazole in the exosite is important for its signal transduction efficiency. In the interest of identifying residues critical for efficacy, the exosite sequence of the 5-HT_2A and D₂ receptors was aligned, with an important difference between the two found at position 5.51, which is Leu in 5-HT_2AR and Phe in D₂R. Mutations in D₂R demonstrated that substitution of F202^5.51 with Leu or Ala reduces the G-protein activity and β-arrestin2 recruitment of aripiprazole. In addition, a derivative of aripiprazole substituted with benzothiazole for the dihydroquinoline ring of D₂R had reduced efficiencies of both G protein activity and β-arrestin2 recruitment. Substitution of L247^5.51F in 5-HT_2AR did not increase the efficacy of aripiprazole. The results suggest that aripiprazole may stabilize different conformations of TM5 and TM6 between the two receptors. Alignment of 5-HT_2AR and D₂R structures (active and inactive) shows that activation of 5-HT_2AR requires a larger downstream swing of W6.48 from the CWxP motif than that observed for D₂R activation. In the 5-HT_2AR, dihydroquinoline is located deeper in the binding pocket interacting with W336^6.48, restricting its movement, whereas it can move gently upon D₂R receptor activation. Similar observations were made for cariprazine. Here, the dynamic coupling between F/L5.51, W6.48 and the PIF motif by the exosite may partly explain why the compounds tested have different efficacies at 5-HT_2AR and D₂R receptors. Compared with inactive and active D₂R constructs, the 5-HT_2AR-aripiprazole complex in the extracellular compartment shows inward movement of TM6, TM7, and ECL2 toward the seven transmembrane cores. These rearrangements suggest that the 5HT_2AR affinity of the bitopic compounds can be reduced by increasing the size of the PP. Changing the arylpiperazine PP to aza-ergoline the authors identified IHCH7009 (D₂R K_i = 33.65 nM, 5-HT_2AR K_i = 3639.15 nM), IHCH7010 (D₂R K_i = 9.03 nM, 5-HT_2AR K_i = 906.78 nM) and IHCH7041 (D₂R K_i = 50.64 nM, 5-HT_2AR K_i = 2371.37 nM) all with very weak 5-HT_2AR affinity. IHCH7041 retains partial agonism of D₂R while IHCH7009 and IHCH7010 are full D₂R agonists (Chen et al., 2021).

Kling et al. investigated the neurotensin receptor type (NTS) 1 receptor crystal structures (White et al., 2012; Egloff et al., 2014) and found that an allosteric binding site was saturated at the C-terminus of NT (8–13). Following sequence analysis, they confirmed that there is a difference between NTS₁R (Arg149^3.32) and NTS₂R (His115^3.32) that may allow for subtype selectivity. Several bitopic ligands of type NT (8–13) were synthesized (Table 2) and compounds (51–56) showed a promising trend in the NTS₁R selectivity. The best compound (54) has K_i value of 1.3 nM associated with 26-fold selectivity towards NTS₂R. Homology modelling and MD simulations confirmed that the compounds bind in a bitopic mode, with NT (8–13) occupying the orthosteric binding site and the amino acid extension occupying the secondary binding site. These results provide a promising starting point for the design of NTS₁R selective agonists (Kling et al., 2019).

Functional Selectivity

Advances in GPCR structural biology and pharmacology have opened up new opportunities for functional drug design. Modulation of GPCRs through allosteric binding sites can alter receptor structure, dynamics and function, resulting in increased spatial and temporal variation. One important aspect of these changes is functional selectivity or otherwise termed biased signalling. Biased signalling can contribute to the enhancement of the intended effect, but can also cause side effects, so one of the most intriguing areas of current research is investigating the functional character of the ligands in different signalling pathways (Hauser et al., 2017).

Egyed at al. reported a systematic study exploring the extracellular SBP to fine-tune the functional profile of D₂R and D₃R ligands. Introduction of the SP increased affinity at both D₂ and D₃ receptors for each ligand. The study demonstrated that the G_i/o and β-arrestin pathways can be specifically modulated from the extracellular vestibule incorporating different SPs to the ligands. Molecular dynamics simulations revealed that G-protein signalling could be linked to the orientation of the PP that is influenced by the SBP binding part of the bitopic compounds (Figure 6). Three PPs and two SPs (Figure 6) were tested using an ethylcyclohexyl linker in analogy to cariprazine. In the G_i/o-mediated signalling pathway, dichlorophenylpiperazine (57) (PP 1) was a partial agonist on both D₂R and D₃R (Table 3). Application of N,N-dimethylurea (SP 1) (cariprazine) also resulted in a partial agonist with significantly increased potency (D₂R pEC₅₀ = 8.85 nM, E_max = 77.4%, D₃R pEC₅₀ = 8.58 nM E_max = 27%). The use of the OtBu motif (SP 2) (61) led to a full agonist, the potency on D₂R was superior to that on D₃R. For 2-methoxyphenylpiperazine (2, 58, 62) (PP 2), no prominent change was observed, all were partial agonists. The 3-(piperazin-1-yl)-5-(trifluoromethyl)benzonitrile (59) (PP 3) with the N,N-dimethylurea SP (60), showed antagonist effects on the G protein coupled signalling pathway of D₂R and D₃R, with an increase in potency. Interestingly, incorporating SP 2 (63) turned the function of PP to a weak partial agonist at both receptors. These results suggest that PP and SP affect functionality together. In the β-arrestin signalling pathway, compounds with SP 2 achieve the largest increase in E_max values, while this was lower for cariprazine. (Table 3). This suggests that cariprazine shows a significant bias towards the G-protein controlled pathway on D₂R. In all cases, the bitopic compounds with 2-methoxyphenylpiperazine PP (2, 58, 62) exhibited antagonist behaviour in contrast to the partial agonism observed in the G-protein coupled signalling pathway. The antagonistic behaviour of 59 was also preserved in the β-arrestin signalling pathway; following the previous trends introduction of any SP led to an increase in pIC₅₀ values here as well. In general, the efficacy data measured at both receptors followed similar trends in both modalities as the receptor affinities (Egyed et al., 2020).

FIGURE 6

Egyed et al., 2020). The binding mode of compound 60 and 63 was extracted from the MD simulations. The simulations revealed that the SP motif influence the position of the PP and that might be linked to the observed different functional profile. The figure representing the binding mode is reproduced from the TOC Figure of our original article Egyed, A et al. Controlling Receptor Function from the Extracellular Vestibule of G-Protein Coupled Receptors. Chem. Commun. 2020, 56 (91), 14167–14170.

High affinity binders, such as 39, 40, 41, 42, 49 (Supplementary Table S5) were also tested for their efficacy on D₃R, both by examining forskolin-dependent inhibition of adenylyl cyclase and by measuring β-arrestin binding. Compounds 42 and 49 were found antagonists in both assays. Compound 41 display functional selectivity, being a weak partial agonist in the adenylyl cyclase assay and a very weak partial agonist/antagonist in the β-arrestin binding assay. Compounds 39 and 40 exhibit weak partial agonism in both the adenylyl cyclase inhibition and β-arrestin binding assays (Lee et al., 2021).

Investigating pure enantiomeric forms of compounds 17–20 (Supplementary Table S3) Tan et al. showed that the (R,R) enantiomers (17b-20b) have a better affinity for D₃R than (S,S) (17a-20a), with the exception of compound 17, which had an identical affinity for both of the enantiomers (17a, 17b) (Tan et al., 2020). The (R,R) isomers (17b-20b) showed weaker affinity (3–20-fold) towards 5-HT_2CR than their (S,S) counterparts (17a-20a). The data suggest that D3R is less sensitive to conformational changes than the 5-HT_2C receptor. Functional studies were also performed with the 17a,b-20a,b (Table 4). Compounds 18–20 were all full or partial agonist on D₃ receptors, whereas for 5-HT_2CR the (S,S) enantiomers (18a-20a) are weak partial agonists, whereas the (R,R) enantiomers (18b-20b) are weak antagonists. Compared to the binding assay, functional results indicate greater selectivity towards D₃R. Furthermore, these compounds showed only very weak partial agonism at 5-HT_2AR and no affinity at 5-HT_2BR. The two enantiomers of compound 17 exhibit opposite behaviour, while (1R,2R)-17b was a potent agonist (EC₅₀ = 3.6 nM, E_max = 77.9%), (1S,2S)-17a was an antagonist on D₃R with a K_i of 16.7 nM, and both derivatives were weak antagonists with micromolar activity on 5-HT_2C receptor. Docking studies suggested a difference between the two compounds (17a,17b) in the orientation of PP. In the case of the agonist (1R,2R)-17b, the 2-methoxy group is deep in the OBP and forms hydrophobic interactions with residues C114^3.36, S196^5.46, and F346^6.52. In the case of the antagonist (1S,2S)-17a, the 2-methoxy group flips out to the extracellular side and the cyclopropane linker between the benzene ring and the protonated N overlays perfectly with the amide linker of eticlopride, which is not present in the agonist. Compounds (1S,2S)-17a, (1R,2R)-18b, (1R,2R)- 19b, and (1R,2R)-20b were inactive in the Tango assay on D₃R, indicating their preference for the G-protein signalling pathway. For further profiling (1R,2R)-17b and (1R,2R)-19b were tested on 29 other aminergic GPCRs that confirmed their good selectivity for D₃R (Tan et al., 2020).

Yan et al. also used PCPMA analogues as PP, with propyl, butyl, pentyl, or cyclohexylethyl linkers, and SP groups taken from aripiprazole, brexipirazole, and cariprazine, respectively. The synthesized library was measured in D₂R binding, D₂R G_i and D₂R β-arrestin BRET assays (Table 5). The starting compound (64) exhibits good affinity (K_i = 61.9 nM) and partial agonist activity in both G_i (EC₅₀ = 49.0 nM, E_max = 25%) and β-arrestin (EC₅₀ = 67.6 nM, E_max = 30%) BRET assays. In comparison, replacement of SP with quinolone (65) increased the potency two-fold with unchanged binding. Changing the linker to propyl (66,67) led to a small decrease in binding affinity but an increase in efficacy (∼10 nM EC₅₀ values and E_max values higher than 50%). Lengthening the linker to 5C units (68,69) led to a decrease in binding affinity and functional activity. The cariprazine-like SP (dimethylamine) and linker (cyclohexyl) with this PP did not show significant activity. The best compound from this series (70) has very potent partial agonist character in both G_i BRET (EC₅₀ = 8.45 nM, E_max = 68%) and β-arrestin2 recruitment assays (EC₅₀ = 9.49 nM, E_max = 16%), with a much lower E_max in the latter. The significant difference between binding affinity and potency for many of these compounds likely reflects the use of an antagonist radioligand [(³H)-N-methylspiperone] in the competitive binding assay, from which an agonist ligand tends to show much lower apparent binding affinity. Attempts have been made to use several PPs but these have been shown to give significantly worse results than the methoxy derivative. In the case of isoquinoline and tetrahydroisoquinoline SP, it was not practical to use the dichlorophenyl motif in the PP (71,72). The best results were obtained with derivatives containing halogen in the meta position on the phenyl group of PP and methoxy in the ortho position (73a,b-76a,b). Pure forms of the enantiomers were also investigated. The majority of the fluorinated derivatives ((1S,2S)-42a, (1R,2R)-73b, (1S,2S)-74a) showed K_i values below 50 nM on binding assay and EC₅₀ values below 20 nM in both G_i and β-arrestin2 BRET assays. The same trend was observed for the chlorinated derivatives [(1S,2S)-75a,(1S,2S)-76a]. Higher E_max was observed for the halogenated derivatives in the G_i signal transduction than in the β-arrestin. After separation of the enantiomers, it was confirmed that the (S,S)-isomers were more efficient in D₂R binding and functional assay. The (R,R) compounds exhibit partial agonist behaviour and the E_max values are higher for G_i signaling. The selectivity of the compounds [(1S,2S)-73a, (1S,2S)-74a, (1S,2S)-75a, (1S,2S)-76a] was investigated on D₁R, D₂R, D₄R, D₅R, 5-HT_1AR, 5-HT_2AR, and 5-HT_2CR, with low selectivity observed towards the D₃ receptor and potent activity on the 5-HT_1A receptor, and good or acceptable selectivity on the other receptors (Table 6). In the case of D₃R, these compounds showed weak partial agonist activity in both G_o and β-arrestin2 BRET assays, albeit with different efficacies. For the 5-HT_1A receptor, all four compounds ((1S,2S)-73a, (1S,2S)-74a, (1S,2S)-75a, (1S,2S)-76a) were similar partial agonists in G_i BRET assays. The lack of selectivity over D₃R and 5-HT_1AR should not be a concern for these compounds, as both D₃R and 5-HT_1AR have been shown to be involved in the therapeutic effects of some antipsychotics. Overall, these four compounds have shown an interesting pharmacological profile (Yan et al., 2021).

TABLE 5

Cmpd	Structure	D2R binding Ki nM (pKi±SEM)	D2R Gαi1 BRET EC50 nM (Emax%) (pEC50 ± SEM)	D2R β-arrestin2 BRET EC50 nM (Emax%) (pEC50 ± SEM)
64		61.9 (7.21 ± 0.04)	49.0 (25 ± 2%) (7.31 ± 0.09)	67.6 (30 ± 1%) (7.17 ± 0.07)
65		59.9 (7.22 ± 0.13)	26.3 (52 ± 1%) (7.58 ± 0.08)	32.4 (53 ± 2%) (7.49 ± 0.14)
66		125.7 (6.90 ± 0.08)	9.30 (58 ± 3%) (8.03 ± 0.01)	10.0 (52 ± 1) (8.00 ± 0.11)
67		155.7 (6.81 ± 0.03)	11.2 (65 ± 3%) (7.95 ± 0.04)	7.08 (60 ± 1%) (8.15 ± 0.12)
68		259.2 (6.59 ± 0.05)	891.2 (12 ± 1%) (6.05 ± 0.42)	416.9 (14 ± 4%) (6.38 ± 0.64)
69		217.8 (6.66 ± 0.08)	77.6 (18 ± 1%) (7.11 ± 0.12)	190.6 (19 ± 1%) (6.72 ± 0.49)
70		977.2 (6.01 ± 0.11)	8.45 (68 ± 1%) (8.07 ± 0.11)	9.49 (16 ± 1%) (8.02 ± 0.06)
71		244.3 (6.61 ± 0.07)	34.8 (51 ± 5%) (7.46 ± 0.10)	94.0 (39 ± 4%) (7.03 ± 0.20)
72		128.1 (6.89 ± 0.112)	14.73 (66 ± 3%) (7.83 ± 0.12)	27.6 (33 ± 1%) (7.56 ± 0.09)
(1S,2S)-73a		20.8 (7.68 ± 0.06)	9.43 (29 ± 3%) (8.03 ± 0.05)	3.63 (18 ± 1%) (8.44 ± 0.17)
(1R,2R)-73b		43.8 (7.36 ± 0.07)	12.9 (13 ± 3%) (7.89 ± 0.14)	1.86 (10 ± 2%) (8.71 ± 0.15)
(1S,2S)-74a		6.58 (8.18 ± 0.04)	4.12 (55 ± 2%) (8.39 ± 0.08)	4.66 (29 ± 1%) (8.33 ± 0.15)
(1R,2R)-74b		362.5 (6.44 ± 0.07)	62.0 (7 ± 1%) (7.21 ± 0.16)	14.7 (17 ± 1%) (7.83 ± 0.12)
(1S,2S)-75a		11.5 (7.94 ± 0.07)	8.9 (40 ± 2%) (8.05 ± 0.04)	2.50 (20 ± 1%) (8.60 ± 0.10)
(1R,2R)-75b		30.1 (7.52 ± 0.02)	NT	NT
(1S,2S)-76a		12.8 (7.89 ± 0.05)	3.41 (71 ± 3%) (8.47 ± 0.08)	8.30 (47 ± 2%) (8.08 ± 0.06)
(1R,2R)-76b		317.0 (6.50 ± 0.04)	197.2 (41 ± 5%) (6.71 ± 0.05)	70.1 (18 ± 3%) (7.15 ± 0.15)

Pharmacological profiling of compounds (D2R binding and functional activity) (Yan et al., 2021).

Schramm et al. investigated the effect of bitopic compounds on muscarinic acetylcholine receptors. Carbachol (CCh) PP was cross-linked to allosteric ligands by linkers of different lengths (1C, 3C, 5C, 8C). The benzoimidazole-piperidine moiety of TBPB [1-(1′-(2-tolyl)-1,4′-bipiperidin-4-yl)-1H-benzo(d)imidazol-2(3H)-one], a known selective bitopic M₁R agonist, and BQCA (benzyl quinolone carboxylic acid) derivatives, that are PAMs, were used as allosteric modulators (Table 7). It was found that BQCA-CCh bitopic compounds act as agonists. The highest potency and efficacy was observed for the compound containing BQCA moiety 81. Comparing with reference compound 86, which does not contain a CCh moiety but only the linker, revealed that the CCh moiety provides some of the agonist activity. In contrast, the TBPB-CCh bitopic ligand (78) showed partial agonism, while the reference 84 was a full agonist. The binding mode of 81 was investigated by docking to an active receptor model. The ammonium group of the CCh moiety forms a charge-assisted hydrogen bond with D105^3.32, while the carbamate carbonyl group serves as a hydrogen bond acceptor for the hydroxyl group of Y408^7.43. This is different from the carbachol binding mode, in which the carbamate structure has a different orientation. The BQCA moiety, located in the region of the extracellular loop, is stabilized by hydrophobic contacts with L174^ECL2 and Y179^ECL2 and a charge-assisted H-bond with K392^ECL3. They concluded that partial agonism through bitopic compounds can be achieved not only by quenching orthosteric receptor activation by an allosteric moiety as in 81 but also by quenching bitopic activation of the receptor by an orthosteric moiety such as CCh in 78 (Schramm et al., 2019).

TABLE 7

Cmpd	N	R	pEC50 nM ± SEM	% Emax±SEM
CCh			6.97 ± 0.03	99 ± 1
TBPB			7.32 ± 0.02	83 ± 1
BQCA			7.20 ± 0.03	90 ± 1
77 (TBPB)	1		n.d.	n.d.
78 (TBPB)	3		5.09 ± 0.24	12 ± 2
79 (TBPB)	6		n.d.	n.d.
80 (BQCA)	1		5.89 ± 0.01	66 ± 0.5
81 (BQCA)	3		6.67 ± 0.02	78 ± 1
82 (BQCA)	6		6.62 ± 0.03	28 ± 0.5
83 (TBPB)	1	H	6.05 ± 0.01	99 ± 1
84 (TBPB)	3	H	6.42 ± 0.01	97 ± 1
85 (TBPB)	6	H	7.38 ± 0.04	98 ± 2
86 (BQCA)	1	H	5.82 ± 0.02	35 ± 1
87 (BQCA)	3	H	n.d.	n.d.
88 (BQCA)	6	H	n.d.	n.d.

Potency and efficacy induced by muscarinic agonists bitopic compounds HEK293t cells overexpressing the M1 receptor (Schramm et al., 2019).

Holze et al. have shown that allosteric coupling of the M₁R can induce conformational changes that affect intracellular signalling. They investigated two groups of M₁R bitopic agonists and varied the length of the linker. Iperoxo, a known agonist, was selected as the PP motif, while two negative allosteric modulators, phtp (89–91) and naph (92–94), were incorporated as SP. (Figure 7) The latter differs from the phth derivative in two main respects: naph contains a larger and branched aliphatic linker. The two pharmacophores were linked by alkyl chains of different length (6–8C) (89–94). While the ligand affinities for the allosteric binding site were very similar within a ligand set, the ligand affinities for the orthosteric binding site depended on the length of the linker, where increasing linker length was correlated with increasing ligand affinity. From this information, it was concluded that the same binding mode was adopted by iperoxo in a series of bitopic compounds driven by its high affinity, and this was confirmed by MD simulations. Therefore, a series of bitopic ligands differing only in the length of the linker may be suitable to investigate the effect of allosteric coupling on signal transduction with subnanometer accuracy. Whereas the longest bitopic agonist, 91, was able to stimulate all three G-protein families, 90 activated G_q/11 and G_s proteins, 89 promoted signal transduction only via G_q/11. 93 and 94 only activated G_q/11 protein signalling, while 92 did not activate any signalling pathway, unlike 89. None of the naph-based ligands were able to activate G_s and G_i/o signalling. These data suggest that different G-proteins show different sensitivities to M₁R activation by these bitopic compounds. While G_q/11 coupling is conserved in almost all bitopic ligands, G_s signalling is promoted only by two members of the phth series. G_i/o activation is particularly sensitive to the bitopic ligand structure with only 91 showing weak M₁R/G_i/o coupling among the compounds tested. MD simulations show that binding of iperoxo results in a complete contraction of the extracellular parts of the ligand binding pocket. In contrast, the bitopic ligands of the phth series bind in such a way that they sterically inhibit the closure of the binding pocket. The extent of the conformational interference depends on the length of the linker and hence the position of the allosteric building block. Since the phth part of 89 is located close to the orthosteric binding site, it inhibits closure, resulting in a more open extracellular conformation. Elongation of the linker with additional methylene groups allowed for subnanometer regulation of the position of the allosteric building block, thereby progressively reducing the closure of the binding pocket, ultimately resulting in greater G-protein binding capacity. FRET measurements have demonstrated that the more closed ligand-binding pocket is associated with greater receptor conformational changes at the G-protein binding surface via an allosteric coupling mechanism. Consistent with this idea, 92, a bitopic ligand with a branched and larger allosteric motif, did not induce conformational changes in M₁R (Holze et al., 2020).

FIGURE 7

Holze et al. (2020).

Wang et al. investigated two naltrexone derivatives substituted with isoquinoline at MOR. The isoquinoline moiety of these bitopic compounds is the SP that interacts with the allosteric site of MOR, and the epoxymorphinan moiety is the PP (Table 8). NAQ has a high affinity for MOR (K_i = 0.55 nM) and high selectivity for κ-opioid receptor (KOR) (48-fold) and δ-opioid receptor (DOR) (241-fold). Compared to DAMGO, it acts as a MOR antagonist in the ³⁵S-GTP [γS]-binding assay with CHO cell lines expressing MOR. It showed less significant withdrawal effects compared to the well-known opioid antagonists naloxone and naltrexone. Similar properties were observed for the compound NCQ (K_i = 0.55, 40-fold KOR, 62-fold DOR selectivity), which shares the same PP part as NAQ and differs only in the SP. NCQ contains a methoxy at position 1 and a chloro functional group at position 4 of isoquinoline. However, in ³⁵S-GTP (γS)-binding assay, NCQ behaved as a partial agonist. MD simulations and free energy calculations proposed that the allosteric part of NAQ and NCQ bind differently in the inactive structure and in the active structure, respectively. Docking studies have shown that the SP parts of NAQ and NCQ may occupy two different subdomains of the allosteric site of MOR, named ABD1 and ABD2. MD simulations were performed with three poses (NAQ inactive, NCQ active and inactive) obtained from the docking calculations and showed that the SP part of NAQ was bound to ABD1 in the inactive MOR. Although the SP motif occupied an allosteric site, no significant modulatory effect was observed on the binding of the PP, similar to the function of a silent allosteric modulator. In the inactive and active MOR the SP of NCQ showed positive allosteric modulation through binding to ABD2. Molecular modelling combined with interaction energy and distance analyses unravelled the molecular mechanisms of allosteric modulation of NAQ and NCQ and emphasized the importance of the chlorine and methoxy substituents of the isoquinoline ring for the allosteric modulatory function of NCQ (Wang et al., 2020).

TABLE 8

Cmpd	Ki (nM±SEM)			MOR vs. KOR	MOR vs. DOR	MOR (35S) GTPγS binding
	MOR	KOR	DOR			EC50 (nM ±SEM)	Emax of DAMGO % ± SEM
	0.55 ± 0.15	26.45 ± 5.22	132.50 ± 27.01	48	241	4.36 ± 0.72	15.83 ± 2.53
	0.55 ± 0.01	22.20 ± 2.10	33.90 ± 0.50	40	62	1.74 ± 0.13	51.00 ± 0.40

Binding affinities and functional efficacies of NAQ and NCQ (Wang et al., 2020).

Binding Kinetics

Although ligand-receptor binding kinetics might have a fundamental role in the development of drug candidates, it is still often overlooked in the early phase of drug discovery. In line with the increased interest in the field, more and more kinetics data (among others association and dissociation rate, residence time, etc.) have been published in the literature, however the magnitude still lags behind the amount of affinity and selectivity data available especially regarding only the allosteric and bitopic ligands. Furthermore, the interpretation of the kinetic data might be hindered by the probe dependence as observed in a prototypical competitive radioligand binding assay for H₁ receptor antagonists, although that aspect is often not considered (Bosma et al., 2019). In line with the relatively limited amount of recent papers, first we refer the readers to recent general review articles on binding kinetics (Sykes et al., 2019; Hoare et al., 2020; Rafael et al., 2020; van der Velden et al., 2020). Very recently a book chapter collecting available kinetic data of GPCR ligands together with experimental evidence for properties that influence the residence time were published (Potterton et al., 2022). The repository enables researchers to analyse the relationship between the structure and the kinetic parameters as well as provides data for the development of predictive algorithms. The authors also outline machine learning workflows to predict residence time. Sykes et al. reviewed recently the literature related to the binding kinetics of GPCR ligands (Sykes et al., 2019). They discussed the theoretical aspects, the experimental methods and their limitations, detailed several factors influencing binding kinetics among others they explored the role of allosteric modulators, that by definition act through the modulation of the binding kinetics of the endogenous or orthosteric ligands. The authors also discuss some molecular level features including shielding the hydrogen bonds from water that affects the binding kinetics.

Although shielding the hydrogen bonds was thought to decrease residence time, in a recent case study on CCR2 receptor, MD simulations of Magarkar et al. suggested that even shielding an intra protein hydrogen bond can enhance the residence time of ligands through the preservation of the binding site rigidity (Magarkar et al., 2019). The ECL2 loop, that is regularly engaged with bitopic compounds, was also proposed to modulate the binding kinetics (Sykes et al., 2019; van der Velden et al., 2020). Already one of the seminal works in the field of modelling the binding pathway to GPCRs, which investigated the binding of three antagonists and an agonists to the β2-adrenoreceptor and one agonist to the β1-adrenoreceptor with MD simulations, highlighted the role of the ECL2 loop and the extracellular vestibule. Interestingly, even the highest barrier of binding often corresponds to the association with the extracellular vestibule even though the binding requires conformational change of the receptor and the ligand has to enter through a narrow passage (Dror et al., 2011). In several receptors, ECL2 were proposed to function as a lid facilitating the entrance and exit of the ligands (Thomas et al., 2016; Wacker et al., 2017; Frank et al., 2018). One of these studies investigated the binding kinetics of cariprazine and aripiprazole. As a prototypical bitopic compounds we exemplify here the effect of the SBP on the binding kinetics through them (Frank et al., 2018). At the D₃ receptor, aripiprazole exhibits a slow monophasic dissociation, while cariprazine displays a rapid biphasic behaviour. Interestingly, in the D₂ receptor both compounds display a slow dissociation. These differences may influence the in vivo action of the drugs. Interactions with ECL2 residues influence the residence time in other receptors like in the β₂ and A_2A receptors, as well (Guo et al., 2016; Masureel et al., 2018). Gaussian accelerated molecular dynamics revealed the role of the ECL2 loop in the formation of allosteric sites for PAMs in the adenosine A₁ receptor (Miao et al., 2018) and unveiled an intermediate binding site between ECL2 and TM1 for caffein in the adenosine A_2A receptor. The authors analysed the effect of more general features like physicochemical properties of the ligand (e.g., lipophilicity) and close contact residue numbers on the drug-receptor dissociation.

Van der Velden et al. summarized structural considerations in relation to binding kinetics presenting the results through four case studies (van der Velden et al., 2020). They showcased the role of the ECL2 loop in the regulation of the ligand kinetics through tiotropium binding to the M₃R and M₂R receptors (Kruse et al., 2012; Tautermann et al., 2013). The more open, flexible ECL2 loop conformation was linked to the shorter residence time observed in the M₂R receptor. Through the example of ZM241385, an A_2A receptor antagonists they highlighted the role of molecular dynamics and mutation experiments in providing structural background for observed kinetics behaviour (Guo et al., 2016). Another example was focused on the β₂ adrenoreceptor. Salmeterol, a bitopic compound displays a 5–7 fold higher residence time compared to salbutamol and epinephrin, both binding only to the orthosteric site (Figure 3C). As salmoterol and salbutamol share the orthosteric binding motif, the interactions in the extracellular site are linked to the increased residence time (Masureel et al., 2018). They also discussed other aspects, like the effect of natural receptor variants, ligand variants and probe dependency.

Riddy et al. investigated the binding kinetics of H₃ receptor antagonists/inverse agonists (Riddy et al., 2019). Although the binding mode of the compounds were not investigated experimentally, they likely form interactions outside the orthosteric pocket, too therefore can be considered bitopic. The different pharmacological profile and the residence time of the compounds might be linked to their preclinical and clinical efficacy. Furthermore, H₃ and off-target sigma-1 receptor occupancy may contribute to paradoxical efficacy of some compounds. In the study of Pedersen et al. (2020) the differential binding kinetics profile of the agonists were not linked to the functional bias, as the bias profile of the selected agonists were not time-dependent and despite the difference in their binding kinetic properties they can display the same degree of bias.

Bitopic compounds and allosteric modulators may directly bind to the secondary binding pocket, however, during the association and dissociation process the secondary site plays a crucial role for the appropriate positioning of all compounds. While experiments rarely shed light on the structural details of binding, molecular dynamics simulations can explore the atomistic process and are useful to predict residence time (Potterton et al., 2019; Decherchi and Cavalli, 2020; Lamim Ribeiro et al., 2020; Salmaso and Jacobson, 2020; Bekker et al., 2021; Kokh and Wade, 2021). Ribeiro and co-workers recently used machine learning and infrequent metadynamics to efficiently predict kinetic rates, transient conformational states, and molecular determinants of drug dissociation on the MOR (Lamim Ribeiro et al., 2020). While both investigated compounds bind to the orthosteric pocket, the transient conformational state for the dissociation was identified around the secondary binding pocket suggesting a key role of the secondary site in the association/dissociation process. In dynamic docking simulations Bekker et al. investigated β₂-adrenoreceptor antagonists identifying several stable and metastable conformational states for the compounds along their association/dissociation path (Bekker et al., 2021). Based on these simulations they propose a way to develop allosteric modulators to inhibit the receptor by blocking the path of the endogenous ligand to the orthosteric site. Metastable binding sites play a crucial role in the study of Gaiser et al. as well (Gaiser et al., 2019). They developed homobivalent bitopic ligands for β₂AR to target the OBP and a previously identified metastable binding site as an allosteric site. Among others the residence time of orthosteric and bitopic A_2A receptor binders was predicted with ensemble based steered molecular dynamics (Potterton et al., 2019). Analysis of the pathways revealed dominant interactions, residues influencing the dissociation time and the calculations proposed that changes in water-ligand energy from the ligand in the binding pocket to the extracellular vestibule was the main factor in the determination of residence time. While hydrophilic ligands are expected to access the orthosteric binding site, that is deeply embedded in the center of the receptor, from the aqueous phase, hydrophobic compounds were proposed to entry through lipid pathways. The examples detailed in this part explore the traditional pathway, however cholesterol and other ligands might enter the receptor from the membrane. As an exciting study we refer to the work of Guixá-González et al. who investigated the cholesterol access to the A_2AR with combined computational and experimental methods. They showed that cholesterol’s impact on A_2AR-binding affinity goes beyond pure allosteric modulation and unveils a new interaction mode between cholesterol and the A2AR (Guixà-González et al., 2017). Similar findings were collected and analysed in a recent review dedicated to the role of the lipid bilayer in the binding of the ligands to the orthosteric and allosteric sites (Szlenk et al., 2019). Even though in this review we focused mainly on the secondary binding pocket in the extracellular vestibule that is accessible through the aqueous phase, some allosteric sites on the receptor surface can only be targeted through the membrane fortifying that investigation of the binding pathways through the membrane is also crucial.

Design Approaches for Allosteric and Bitopic Compounds

During the previous sections we often pointed out the value of computational approaches in the investigation of both allosteric and bitopic compounds. Due to the tremendous number of studies a comprehensive overview of the computational approaches to design allosteric (Wold et al., 2019; Chatzigoulas and Cournia, 2021) and bitopic ligands (Newman et al., 20162020; Fronik et al., 2017) for GPCRs warrant a separate review (Basith et al., 2018; Raschka and Kaufman, 2020; Ballante et al., 2021), we could only highlight here a few important studies to draw attention towards their usefulness in drug discovery settings (Dehua Yang et al., 2021).

Allosteric sites are less conserved and therefore they can be exploited to design ligands with high selectivity and modalities that could not be achieved from the orthosteric site. The increasing number of experimental GPCR structures urges the use of structure-based methods. However, the identification of the allosteric sites remains challenging as they often form fully only in the presence of an allosteric ligand following an induced fit mechanism. Nevertheless, several computational approach were developed to facilitate the spotting of new allosteric sites like Allosite (Huang et al., 2013), AlloFinder (Huang et al., 2018), ExProSE (Greener et al., 2017), Fpocket (Le Guilloux et al., 2009; Schmidtke et al., 2010), FTmap (Brenke et al., 2009; Kozakov et al., 2015), GRID (Goodford, 1985), LIGSITE^csc (Huang and Schroeder, 2006), SiteMap (Halgren, 2009) and MixMD (Ghanakota and Carlson, 2016). FTMap and FTSite was recently shown to perform well on identifying GPCR allosteric sites with limitations on those occurring on the protein-membrane interface that could be attributed to the development of the program originally for soluble globular proteins (Wakefield et al., 2019).

Even after the identification of the allosteric site, simple docking might not always be successful due to induced fit binding. Furthermore, allosteric modulators are prone to “steep” SAR, obscure relationship between the binding affinity and functional effect and slow kinetics (on and/or off rates) that hinders their discovery and design (Congreve et al., 2017b). Huang et al. developed a protocol combining homology modelling and docking to find novel allosteric modulators of the orphan GPR68 and GPR65 receptors (Huang et al., 2015). They generated over three thousand homology models and docked their experimentally validated active compound lorazepam and decoy compounds to identify putative binding sites. They optimized the binding site around the bound ligand and redocked the ligand and the decoys again until a stable docking mode emerged. That plausible binding site was utilized to dock over 3.1 million lead-like compounds. From the selected 17 hits four increased cAMP production. Docking close analogues of the hit compounds lead to another 25 compounds for testing among them 13 with higher activity than the reference compound lorazepam. Similar protocol was utilized for the GPR65 receptor as well showing that the protocol might be applied to a broader field. While this protocol might be applied to several—even orphan—GPCRs it requires at least one experimentally determined known active compound that might be hard to get for other orphan GPCRs and close enough homology to templates that warrant the homology modelling. Nevertheless, this is a great example how the combined experimental and computational approaches can lead to the identification of novel allosteric modulators even for orphan GPCR targets. Miao et al. focused on the identification of novel, chemically diverse allosteric modulators of the M₂ receptor (Miao et al., 2016). The authors used accelerated molecular dynamics to account for receptor flexibility and to generate an ensemble of structures for docking. After retrospective validation virtual screening coupled with induced fit docking (IFD) was applied to select compounds targeting the IXO-nanobody-bound active and the QNB-bound inactive M₂ mAChR for testing. The method successfully identified both positive and negative allosteric modulators and clearly demonstrate that accounting for receptor flexibility is a key in the discovery of allosteric modulators. Nevertheless, for less flexible binding sites even simple docking protocols might be plausible as demonstrated by Korczynska et al. identifying a positive allosteric modulator that potentiates antagonist binding leading to subtype selectivity at the M₂ muscarinic acetylcholine receptor (Korczynska et al., 2018). Since allosteric modulators are often small and rigid compounds, fragment based approaches (Keserű et al., 2016) emerge as a plausible choice for the design that is supported by several successful application (Christopher et al., 2015; Orgován et al., 2019). Furthermore, covalent approaches should not be overlooked either to aid structurally informed rational design (Lu and Zhang, 2017; Bian et al., 2020; Wenchao Lu et al., 2021).

Bitopic compounds are in the forefront of drug development for GPCRs as they can combine the advantages of targeting the orthosteric and a secondary site (Newman et al., 20162020; Fronik et al., 2017). Fragment based methods are often applied to design novel bitopic compounds (Vass et al., 2014; Egyed et al., 2021). Recently our group have developed a computational protocol to design specific, selective receptor ligands (Egyed et al., 2021). First fragments were docked to the orthosteric binding site of the receptors available in experimental structures (D₃: PDB ID: 3PBL (Chien et al., 2010), 5-HT_1B: PDB ID: 4IAQ (Wang et al., 2013), 4IAR (Wang et al., 2013); 5-HT_2B: PDB ID: 4IB4 (Wacker et al., 2013), 4MC3 (Liu et al., 2013); H₁: PDB ID 3RZE (Shimamura et al., 2011) and M₁: PDB ID: 5CXV (Thal et al., 2016)), or a homology model in case of the D₂ receptor. Then, virtual fragment screening was performed against the secondary binding site of the combined protein-ligand complex. The identified SBP fragment was then linked to the OBP core by a linker. As a control, the resulting bitopic compounds were docked back into the initial crystal structure. This protocol has been validated by designing selective D₂/D₃, 5-HT_1B/5-HT_2B and H₁/M₁ receptors. Docking-based fragment evolution approach utilizes the same methodology as exemplified on the design of β₁ and β₂ receptor bitopic compounds (Chevillard et al., 2021). The fragment evolution protocol merges fragment growing with a matrix-based strategy that was originally implemented for potency optimization (Chevillard et al., 2019). First, possible OBP fragments were docked and they were evaluated using the concept of ligand efficiency. Next, fragment growing surrogates suitable for reactive alkylation were defined and docked to the secondary binding pocket. Surrogates that overlap with the core OBP fragment or was marked favourably in both receptors were removed from the top ranked compounds, the remaining top surrogates were kept for further investigation. The OBP fragments were reacted in silico with the surrogates, the resulting compounds were docked into the receptors to ensure pose fidelity. Based on these calculations the best surrogates were selected as secondary binding motif for the β₁ and β₂ receptor, respectively. The approach was validated by the synthesis and experimental evaluation of the designed compounds. Classical docking and virtual screening approaches could be also utilized for the development of bitopic compounds (Cao et al., 2018) and even to develop fluorescent GPCR probes (Prokop et al., 2021). We highlight here a study that utilized structure guided design of GPCR polypharmacology (Kampen et al., 2021). Kampen et al. aimed to design dual A_2A/D₂ bitopic compounds that was very challenging due to the significantly different binding sites of the receptors. First, docking based structural analysis confirmed that dual-target ligands of the A_2AR and D₂R could be obtained by targeting the orthosteric and secondary pockets. Then, they designed potential dual targeting virtual chemical libraries that could be rapidly synthesized. The prepared libraries were screened virtually with docking on the A_2A and the D₂ receptor to select hits. From them one promising compound was selected and developed further with SAR investigations.

Discussing the recent advances in the allosteric and bitopic field we pointed out several times the usefulness of MD based methods. These simulations can explore the differences in the interaction patterns of congeneric molecules more sensitively compared to docking that could be important to understand the different functional outcome of these ligands (Egyed et al., 2020) and to design compounds with specific pharmacological profile (McCorvy et al., 2018) and they might reveal cryptic pockets opened by ligands (Ferruz et al., 2018) that might be overlooked in simple docking calculations. Mutation studies combined with extensive molecular dynamics modelling the dissociation of the ligands was utilized to clarify the structural basis of the long duration of action and kinetic selectivity of tiotropium for the M₃ receptor (Tautermann et al., 2013). A similar study aimed to clarify the molecular determinants of the bitopic binding mode of a negative allosteric modulator of the dopamine D₂ receptor (Draper-Joyce et al., 2018). MDs combined with docking linked the degree of closure of the extracellular loop region to the extent of ligand bias and highlighted the importance of the appropriate receptor conformation for virtual screening at the 5-HT_2B receptor (Denzinger et al., 2020). A similar concept was presented by Bermudez et al. proposing that agonists with extended binding modes selectively interfere with binding pocket closure and through divergent allosteric coupling that leads to ligand bias (Bermudez and Bock, 2019).

The structure-based methods clearly benefit from the increase of published GPCR structures, especially that more and more active structures are available, however the design still remains challenging. Nevertheless, with more template available for homology modelling and the publication of AlphaFold (Jumper et al., 2021) facilitate the structure-based methods for targets previously out of scope for these methods broadening the applicability spectrum. While we mainly highlighted structure based approaches classical ligand based methods and cheminformatics also contribute to the development of bitopic GPCR ligands (Basith et al., 2018; James and Heifetz, 2018; Raschka and Kaufman, 2020).

Conclusion

GPCRs are one of the largest families of receptors and are among the most targeted proteins for drug discovery. One of the major challenges in the field is the identification of subtype and functionally selective compounds with high potency, designed efficacy and appropriate binding kinetics profile, which are essential to avoid side effects. The secondary binding pocket plays a prominent role in achieving selectivity, while orthosteric ligands are mainly responsible for affinity and functional activity. Bitopic compounds combine the properties of orthosteric and allosteric pharmacophores. With the continuous expansion of available GPCR structures, the secondary binding sites of the receptors are becoming better understood, allowing the construction of complex ligands with designed pharmacological profile. In this review, we have provided an insight into allosteric modulators of class A GPCRs and a detailed review of bitopic compounds that have been released in the last years. We have highlighted the influence of the secondary site in affinity, selectivity, functional selectivity and binding kinetics. The increasing amount of pharmacological data and new structures together with appropriate modelling tools can contribute to the design of allosteric and bitopic drug candidates with an optimized pharmacology profile and thus accelerating the drug discovery against diseases with high unmet medical need.

References

• 1

  AckleyM. A. (2008). Morpholine Dopamine Agonists for the Treatment of Pain. WO 087512 A1.

  • Google Scholar
  • View reference in article

• 2

  Ágai-CsongorÉ.DományG.NógrádiK.GalambosJ.VágóI.KeserűG. M.et al (2012). Discovery of Cariprazine (RGH-188): A Novel Antipsychotic Acting on Dopamine D3/D2 Receptors. Bioorg. Med. Chem. Lett.22 (10), 3437–3440. 10.1016/j.bmcl.2012.03.104

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AllertonC. M. N.CookA. S.HepworthD.MillerD. C. (2005). Aminopyridine Derivatives as Selective Dopamine D3 Agonists. WO 115985 Al.

  • Google Scholar
  • View reference in article

• 4

  BallanteF.KooistraA. J.KampenS.de GraafC.CarlssonJ. (2021). Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You. Pharmacol. Rev.73 (4), 1698–1736. 10.1124/pharmrev.120.000246

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BarnesN. M.AhernG. P.BecamelC.BockaertJ.CamilleriM.Chaumont-DubelS.et al (2021). International Union of Basic and Clinical Pharmacology. CX. Classification of Receptors for 5-Hydroxytryptamine; Pharmacology and Function. Pharmacol. Rev.73 (1), 310–520. 10.1124/pr.118.015552

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BarresiE.MartiniC.Da SettimoF.GrecoG.TalianiS.GiacomelliC.et al (2021a). Allosterism vs. Orthosterism: Recent Findings and Future Perspectives on A2B AR Physio-Pathological Implications. Front. Pharmacol.12, 384. 10.3389/fphar.2021.652121

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BarresiE.GiacomelliC.MarchettiL.BagliniE.SalernoS.GrecoG.et al (2021b). Novel Positive Allosteric Modulators of A2B Adenosine Receptor Acting as Bone Mineralisation Promoters. J. Enzyme Inhib. Med. Chem.36 (1), 286–294. 10.1080/14756366.2020.1862103

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BasithS.CuiM.MacalinoS. J. Y.ParkJ.ClavioN. A. B.KangS.et al (2018). Exploring G Protein-Coupled Receptors (GPCRs) Ligand Space via Cheminformatics Approaches: Impact on Rational Drug Design. Front. Pharmacol.9, 128. 10.3389/fphar.2018.00128

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BattitiF. O.CemajS. L.GuerreroA. M.ShaikA. B.LamJ.RaisR.et al (2019). The Significance of Chirality in Drug Design and Synthesis of Bitopic Ligands as D3 Receptor (D3R) Selective Agonists. J. Med. Chem.62 (13), 6287–6314. 10.1021/acs.jmedchem.9b00702

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BattitiF. O.NewmanA. H.BonifaziA. (2020). Exception that Proves the Rule: Investigation of Privileged Stereochemistry in Designing Dopamine D3R Bitopic Agonists. ACS Med. Chem. Lett.11 (10), 1956–1964. 10.1021/acsmedchemlett.9b00660

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  BekkerG. J.ArakiM.OshimaK.OkunoY.KamiyaN. (2021). Accurate Binding Configuration Prediction of a G-Protein-Coupled Receptor to its Antagonist Using Multicanonical Molecular Dynamics-Based Dynamic Docking. J. Chem. Inf. Model.61 (10), 5161–5171. 10.1021/acs.jcim.1c00712

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  BermudezM.BockA. (2019). Does Divergent Binding Pocket Closure Drive Ligand Bias for Class A GPCRs. Trends Pharmacol. Sci.40 (4), 236–239. 10.1016/j.tips.2019.02.005

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  BethanyA. R.HuiqunW.YanZ. (2019). Recent Advances in the Drug Discovery and Development of Dualsteric/Bitopic Activators of G Protein-Coupled Receptors. Curr. Top. Med. Chem.19 (26), 2378–2392. 10.2174/1568026619666191009164609

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  BianY.JunJ. J.CuylerJ.XieX. Q. (2020). Covalent Allosteric Modulation: An Emerging Strategy for GPCRs Drug Discovery. Eur. J. Med. Chem.206, 112690. 10.1016/j.ejmech.2020.112690

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  BockA.SchrageR.MohrK. (2018). Allosteric Modulators Targeting CNS Muscarinic Receptors. Neuropharmacology136 (Pt C), 427–437. 10.1016/j.neuropharm.2017.09.024

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  BosmaR.StoddartL. A.GeorgiV.Bouzo-LorenzoM.BushbyN.InkoomL.et al (2019). Probe Dependency in the Determination of Ligand Binding Kinetics at a Prototypical G Protein-Coupled Receptor. Sci. Rep.9 (1), 7906. 10.1038/s41598-019-44025-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  BrenkeR.KozakovD.ChuangG. Y.BeglovD.HallD.LandonM. R.et al (2009). Fragment-based Identification of Druggable 'hot Spots' of Proteins Using Fourier Domain Correlation Techniques. BIOINFORMATICS25 (5), 621–627. 10.1093/bioinformatics/btp036

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  British Pharmacological Society (2021). Dopamine Receptors – IUPHAR Review 13 - Beaulieu - 2015 - British Journal of Pharmacology - Wiley Online Library. Availablr at: https://bpspubs.onlinelibrary.wiley.com/doi/10.1111/bph.12906 (accessed 12 2021, 19).

  • Google Scholar
  • View reference in article

• 19

  BurrisK. D.MolskiT. F.XuC.RyanE.TottoriK.KikuchiT.et al (2002). Aripiprazole, a Novel Antipsychotic, Is a High-Affinity Partial Agonist at Human Dopamine D2 Receptors. J. Pharmacol. Exp. Ther.302 (1), 381–389. 10.1124/jpet.102.033175

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  CaoY.SunN.ZhangJ.LiuZ.TangY. Z.WuZ.et al (2018). Design, Synthesis, and Evaluation of Bitopic Arylpiperazine-Phthalimides as Selective Dopamine D3 Receptor Agonists. MedChemComm9 (9), 1457–1465. 10.1039/C8MD00237A

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  ChanH. C. S.LiY.DahounT.VogelH.YuanS. (2019). New Binding Sites, New Opportunities for GPCR Drug Discovery. Trends Biochem. Sci.44 (4), 312–330. 10.1016/j.tibs.2018.11.011

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  ChatzigoulasA.CourniaZ. (2021). Rational Design of Allosteric Modulators: Challenges and Successes. Wires Comput. Mol. Sci.11 (6), e1529. 10.1002/wcms.1529

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  CheT.RothB. L. (2021). Structural Insights Accelerate the Discovery of Opioid Alternatives. Annu. Rev. Biochem.90, 739–761. 10.1146/annurev-biochem-061620-044044

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  ChenZ.FanL.WangH.YuJ.LuD.QiJ.et al (2021). Structure-Based Design of a Novel Third-Generation Antipsychotic Drug Lead with Potential Antidepressant Properties. Nat. Neurosci.25, 39–49. 10.1038/s41593-021-00971-w

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  ChengJ.GiguèreP. M.OnajoleO. K.LvW.GaisinA.GunosewoyoH.et al (2015). Optimization of 2-Phenylcyclopropylmethylamines as Selective Serotonin 2C Receptor Agonists and Their Evaluation as Potential Antipsychotic Agents. J. Med. Chem.58 (4), 1992–2002. 10.1021/jm5019274

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  ChengJ.GiguereP. M.SchmerbergC. M.PogorelovV. M.RodriguizR. M.HuangX. P.et al (2016a). Further Advances in Optimizing (2-Phenylcyclopropyl)Methylamines as Novel Serotonin 2C Agonists: Effects on Hyperlocomotion, Prepulse Inhibition, and Cognition Models. J. Med. Chem.59 (2), 578–591. 10.1021/acs.jmedchem.5b01153

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  ChengJ.McCorvyJ. D.GiguereP. M.ZhuH.KenakinT.RothB. L.et al (2016b). Design and Discovery of Functionally Selective Serotonin 2C (5-HT2C) Receptor Agonists. J. Med. Chem.59 (21), 9866–9880. 10.1021/acs.jmedchem.6b01194

  • CrossRef
  • Google Scholar

• 28

  ChevillardF.StotaniS.KarawajczykA.HristevaS.PardonE.SteyaertJ.et al (2019). Interrogating Dense Ligand Chemical Space with a Forward-Synthetic Library. Proc. Natl. Acad. Sci. U. S. A.116 (23), 11496–11501. 10.1073/pnas.1818718116

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  ChevillardF.KelemenÁ.BakerJ. G.AranyodiV. A.BalzerF.KolbV.et al (2021). Fragment Evolution for GPCRs: The Role of Secondary Binding Sites in Optimization. Chem. Commun.57 (81), 10516–10519. 10.1039/D1CC04636E

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  ChienE. Y.LiuW.ZhaoQ.KatritchV.HanG. W.HansonM. A.et al (2010). Structure of the Human Dopamine D3 Receptor in Complex with a D2/D3 Selective Antagonist. Science330 (6007), 1091–1095. 10.1126/science.1197410

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  ChristopherJ. A.AvesS. J.BennettK. A.DoréA. S.ErreyJ. C.JazayeriA.et al (2015). Fragment and Structure-Based Drug Discovery for a Class C GPCR: Discovery of the MGlu5 Negative Allosteric Modulator HTL14242 (3-Chloro-5-[6-(5-Fluoropyridin-2-Yl)Pyrimidin-4-Yl]Benzonitrile). J. Med. Chem.58 (16), 6653–6664. 10.1021/acs.jmedchem.5b00892

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  ChristopoulosA. (2014). Advances in G Protein-Coupled Receptor Allostery: From Function to Structure. Mol. Pharmacol.86 (5), 463–478. 10.1124/mol.114.094342

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  CongreveM.OswaldC.MarshallF. H. (2017). Applying Structure-Based Drug Design Approaches to Allosteric Modulators of GPCRs. Trends Pharmacol. Sci.38 (9), 837–847. 10.1016/j.tips.2017.05.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  CongreveM.OswaldC.MarshallF. H. (2017). Applying Structure-Based Drug Design Approaches to Allosteric Modulators of GPCRs. Trends Pharmacol. Sci.38, 837–847. 10.1016/j.tips.2017.05.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  ConnP. J.JonesC. K.LindsleyC. W. (2009). Subtype-Selective Allosteric Modulators of Muscarinic Receptors for the Treatment of CNS Disorders. Trends Pharmacol. Sci.30 (3), 148–155. 10.1016/j.tips.2008.12.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  CordeauxY.IJzermanA. P.HillS. J. (2004). Coupling of the Human A1 Adenosine Receptor to Different Heterotrimeric G Proteins: Evidence for Agonist-specific G Protein Activation. Br. J. Pharmacol.143 (6), 705–714. 10.1038/sj.bjp.0705925

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  DecherchiS.CavalliA. (2020). Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation. Chem. Rev.120 (23), 12788–12833. 10.1021/acs.chemrev.0c00534

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  DenzingerK.NguyenT. N.NoonanT.WolberG.BermudezM. (2020). Biased Ligands Differentially Shape the Conformation of the Extracellular Loop Region in 5-HT2B Receptors. Ijms21 (24), 9728. 10.3390/ijms21249728

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  Draper-JoyceC. J.MichinoM.VermaR. K.Klein HerenbrinkC.ShonbergJ.KopinathanA.et al (2018). The Structural Determinants of the Bitopic Binding Mode of a Negative Allosteric Modulator of the Dopamine D2 Receptor. Biochem. Pharmacol.148, 315–328. 10.1016/j.bcp.2018.01.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  DrorR. O.PanA. C.ArlowD. H.BorhaniD. W.MaragakisP.ShanY.et al (2011). Pathway and Mechanism of Drug Binding to G-Protein-Coupled Receptors. Proc. Natl. Acad. Sci. U. S. A.108 (32), 13118–13123. 10.1073/pnas.1104614108

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  EgloffP.HillenbrandM.KlenkC.BatyukA.HeineP.BaladaS.et al (2014). Structure of Signaling-Competent Neurotensin Receptor 1 Obtained by Directed Evolution in Escherichia Coli. Proc. Natl. Acad. Sci. U S A.111 (6), E655–E662. 10.1073/pnas.1317903111

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  EgyedA.Domány-KovácsK.KoványiB.HortiF.KurkóD.KissD. J.et al (2020). Controlling Receptor Function from the Extracellular Vestibule of G-Protein Coupled Receptors. Chem. Commun. (Camb)56 (91), 14167–14170. 10.1039/D0CC05532H

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  EgyedA.KelemenÁ. A.VassM.VisegrádyA.TheeS. A.WangZ.et al (2021). Controlling the Selectivity of Aminergic GPCR Ligands from the Extracellular Vestibule. Bioorg. Chem.111, 104832. 10.1016/j.bioorg.2021.104832

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  FelderC. C.JoyceK. E.BrileyE. M.MansouriJ.MackieK.BlondO.et al (1995). Comparison of the Pharmacology and Signal Transduction of the Human Cannabinoid CB1 and CB2 Receptors. Mol. Pharmacol.48 (3), 443–450.

  • Google Scholar
  • View reference in article

• 45

  FerruzN.DoerrS.Vanase-FrawleyM. A.ZouY.ChenX.MarrE. S.et al (2018). Dopamine D3 Receptor Antagonist Reveals a Cryptic Pocket in Aminergic GPCRs. Sci. Rep.8, 897. 10.1038/s41598-018-19345-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  FrankA.KissD. J.KeserűG. M.StarkH. (2018). Binding Kinetics of Cariprazine and Aripiprazole at the Dopamine D3 Receptor. Sci. Rep.8, 12509. 10.1038/s41598-018-30794-y

  • CrossRef
  • Google Scholar

• 47

  FreiJ. N.BroadhurstR. W.BostockM. J.SoltA.JonesA. J. Y.GabrielF.et al (2020). Conformational Plasticity of Ligand-Bound and Ternary GPCR Complexes Studied by 19F NMR of the β1-adrenergic Receptor. Nat. Commun.11 (1), 669. 10.1038/s41467-020-14526-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  FronikP.GaiserB. I.Sejer PedersenD. (2017). Bitopic Ligands and Metastable Binding Sites: Opportunities for G Protein-Coupled Receptor (GPCR) Medicinal Chemistry. J. Med. Chem.60 (10), 4126–4134. 10.1021/acs.jmedchem.6b01601

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  GaiserB. I.DanielsenM.Marcher-RørstedE.Røpke JørgensenK.WróbelT. M.FrykmanM.et al (2019). Probing the Existence of a Metastable Binding Site at the β2-Adrenergic Receptor with Homobivalent Bitopic Ligands. J. Med. Chem.62 (17), 7806–7839. 10.1021/acs.jmedchem.9b00595

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  GaoY.PetersonS.MasriB.HouglandM. T.AdhamN.GyertyánI.et al (2014). Cariprazine Exerts Antimanic Properties and Interferes with Dopamine D2 Receptor β-arrestin Interactions. Pharmacol. Res. Perspect.3 (1), e00073. 10.1002/prp2.73

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  GaoZ. G.InoueA.JacobsonK. A. (2018). On the G Protein-Coupling Selectivity of the Native A2B Adenosine Receptor. Biochem. Pharmacol.151, 201–213. 10.1016/j.bcp.2017.12.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  GhanakotaP.CarlsonH. A. (2016). Moving beyond Active-Site Detection: MixMD Applied to Allosteric Systems. J. Phys. Chem. B120 (33), 8685–8695. 10.1021/acs.jpcb.6b03515

  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  GoodfordP. J. (1985). A Computational Procedure for Determining Energetically Favorable Binding Sites on Biologically Important Macromolecules. J. Med. Chem.28 (7), 849–857. 10.1021/jm00145a002

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  GreenerJ. G.FilippisI.SternbergM. J. E. (2017). Predicting Protein Dynamics and Allostery Using Multi-Protein Atomic Distance Constraints. Structure25 (3), 546–558. 10.1016/j.str.2017.01.008

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  GrundmannM.BenderE.SchambergerJ.EitnerF. (2021). Pharmacology of Free Fatty Acid Receptors and Their Allosteric Modulators. Int. J. Mol. Sci.22 (4), 1763. 10.3390/ijms22041763

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  Guixà-GonzálezR.AlbasanzJ. L.Rodriguez-EspigaresI.PastorM.SanzF.Martí-SolanoM.et al (2017). Membrane Cholesterol Access into a G-Protein-Coupled Receptor. Nat. Commun.8 (1), 14505. 10.1038/ncomms14505

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  GuoD.PanA. C.DrorR. O.MockingT.LiuR.HeitmanL. H.et al (2016). Molecular Basis of Ligand Dissociation from the Adenosine A2A Receptor. Mol. Pharmacol.89 (5), 485–491. 10.1124/mol.115.102657

  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  HagaK.KruseA. C.AsadaH.Yurugi-KobayashiT.ShiroishiM.ZhangC.et al (2012). Structure of the Human M2 Muscarinic Acetylcholine Receptor Bound to an Antagonist. Nature482 (7386), 547–551. 10.1038/nature10753

  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  HalgrenT. A. (2009). Identifying and Characterizing Binding Sites and Assessing Druggability. J. Chem. Inf. Model.49 (2), 377–389. 10.1021/ci800324m

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  HauserA. S.AttwoodM. M.Rask-AndersenM.SchiöthH. B.GloriamD. E. (2017). Trends in GPCR Drug Discovery: New Agents, Targets and Indications. Nat. Rev. Drug Discov.16 (12), 829–842. 10.1038/nrd.2017.178

  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  HoareS. R. J.TewsonP. H.QuinnA. M.HughesT. E.BridgeL. J. (2020). Analyzing Kinetic Signaling Data for G-Protein-Coupled Receptors. Sci. Rep.10 (1), 12263. 10.1038/s41598-020-67844-3

  • CrossRef
  • Google Scholar

• 62

  HolzeJ.BermudezM.PfeilE. M.KaukM.BödefeldT.IrmenM.et al (2020). Ligand-Specific Allosteric Coupling Controls G-Protein-Coupled Receptor Signaling. ACS Pharmacol. Transl. Sci.3 (5), 859–867. 10.1021/acsptsci.0c00069

  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  HowlettA. C.ChampionT. M.WilkenG. H.MechoulamR. (1990). Stereochemical Effects of 11-OH-Delta 8-Tetrahydrocannabinol-Dimethylheptyl to Inhibit Adenylate Cyclase and Bind to the Cannabinoid Receptor. Neuropharmacology29 (2), 161–165. 10.1016/0028-3908(90)90056-w

  • CrossRef
  • Google Scholar

• 64

  HuangB.SchroederM. (2006). LIGSITEcsc: Predicting Ligand Binding Sites Using the Connolly Surface and Degree of Conservation. BMC Struct. Biol.6, 19. 10.1186/1472-6807-6-19

  • CrossRef
  • Google Scholar

• 65

  HuangM.SongK.LiuX.LuS.ShenQ.WangR.et al (2018). AlloFinder: A Strategy for Allosteric Modulator Discovery and Allosterome Analyses. Nucleic Acids Res.46, W451–W458. 10.1093/nar/gky374

  • CrossRef
  • Google Scholar

• 66

  HuangS.XuP.TanY.YouC.ZhangY.JiangY.et al (2021). Structural Basis for Recognition of Anti-migraine Drug Lasmiditan by the Serotonin Receptor 5-Ht1f-G Protein Complex. Cell Res31 (9), 1036–1038. 10.1038/s41422-021-00527-4

  • CrossRef
  • Google Scholar

• 67

  HuangW.LuS.HuangZ.LiuX.MouL.LuoY.et al (2013). Allosite: A Method for Predicting Allosteric Sites. Bioinformatics29 (18), 2357–2359. 10.1093/bioinformatics/btt399

  • CrossRef
  • Google Scholar

• 68

  HuangX. P.KarpiakJ.KroezeW. K.ZhuH.ChenX.MoyS. S.et al (2015). Allosteric Ligands for the Pharmacologically Dark Receptors GPR68 and GPR65. Nature527 (7579), 477–483. 10.1038/nature15699

  • CrossRef
  • Google Scholar

• 69

  JakubikJ.El-FakahanyE. E. (2020). Current Advances in Allosteric Modulation of Muscarinic Receptors. Biomolecules10 (2), 325. 10.3390/biom10020325

  • CrossRef
  • Google Scholar

• 70

  JakubíkJ.El-FakahanyE. E. (2021). Allosteric Modulation of GPCRs of Class A by Cholesterol. Ijms22 (4), 1953. 10.3390/ijms22041953

  • CrossRef
  • Google Scholar

• 71

  JamesT. (2018). “Cheminformatics in the Service of GPCR Drug Discovery,” in Computational Methods for GPCR Drug Discovery. Editor HeifetzA. (New York, NY: Springer), 395–411. Methods in Molecular Biology. 10.1007/978-1-4939-7465-8_20

  • CrossRef
  • Google Scholar

• 72

  JumperJ.EvansR.PritzelA.GreenT.FigurnovM.RonnebergerO.et al (2021). Highly Accurate Protein Structure Prediction with AlphaFold. Nature596 (7873), 583–589. 10.1038/s41586-021-03819-2

  • CrossRef
  • Google Scholar

• 73

  KampenS.Duy VoD.ZhangX.PanelN.YangY.JaitehM.et al (2021). Structure‐Guided Design of G‐Protein‐Coupled Receptor Polypharmacology. Angew. Chem. Int. Ed.60 (33), 18022–18030. 10.1002/anie.202101478

  • CrossRef
  • Google Scholar

• 74

  KapurA.ZhaoP.SharirH.BaiY.CaronM. G.BarakL. S.et al (2009). Atypical Responsiveness of the Orphan Receptor GPR55 to Cannabinoid Ligands. J. Biol. Chem.284 (43), 29817–29827. 10.1074/jbc.M109.050187

  • CrossRef
  • Google Scholar

• 75

  KeserűG. M.ErlansonD. A.FerenczyG. G.HannM. M.MurrayC. W.PickettS. D. (2016). Design Principles for Fragment Libraries: Maximizing the Value of Learnings from Pharma Fragment-Based Drug Discovery (FBDD) Programs for Use in Academia. J. Med. Chem.59, 8189–8206. 10.1021/acs.jmedchem.6b00197

  • CrossRef
  • Google Scholar

• 76

  KimK.CheT.PanovaO.DiBertoJ. F.LyuJ.KrummB. E.et al (2020). Structure of a Hallucinogen-Activated Gq-Coupled 5-HT2A Serotonin Receptor. Cell182 (6), 1574–e19. e19. 10.1016/j.cell.2020.08.024

  • CrossRef
  • Google Scholar

• 77

  KissB.HorváthA.NémethyZ.SchmidtE.LaszlovszkyI.BugovicsG.et al (2010). Cariprazine (RGH-188), a Dopamine D(3) Receptor-Preferring, D(3)/D(2) Dopamine Receptor Antagonist-Partial Agonist Antipsychotic Candidate: In Vitro and Neurochemical Profile. J. Pharmacol. Exp. Ther.333 (1), 328–340. 10.1124/jpet.109.160432

  • CrossRef
  • Google Scholar

• 78

  KlingR. C.BurchardtC.EinsiedelJ.HübnerH.GmeinerP. (2019). Structure-Based Exploration of an Allosteric Binding Pocket in the NTS1 Receptor Using Bitopic NT(8-13) Derivatives and Molecular Dynamics Simulations. J. Mol. Model.25 (7), 193. 10.1007/s00894-019-4064-x

  • CrossRef
  • Google Scholar

• 79

  KokhD. B.WadeR. C. (2021). G Protein-Coupled Receptor-Ligand Dissociation Rates and Mechanisms from τRAMD Simulations. J. Chem. Theor. Comput.17 (10), 6610–6623. 10.1021/acs.jctc.1c00641

  • CrossRef
  • Google Scholar

• 80

  KooistraA. J.MordalskiS.Pándy-SzekeresG.EsguerraM.MamyrbekovA.MunkC.et al (2021). GPCRdb in 2021: Integrating GPCR Sequence, Structure and Function. Nucleic Acids Res.49 (D1), D335–D343. 10.1093/nar/gkaa1080

  • CrossRef
  • Google Scholar

• 81

  KorczynskaM.ClarkM. J.ValantC.XuJ.MooE. V.AlboldS.et al (2018). Structure-Based Discovery of Selective Positive Allosteric Modulators of Antagonists for the M2 Muscarinic Acetylcholine Receptor. Proc. Natl. Acad. Sci. U S A.115 (10), E2419–E2428. 10.1073/pnas.1718037115

  • CrossRef
  • Google Scholar

• 82

  KozakovD.GroveL. E.HallD. R.BohnuudT.MottarellaS. E.LuoL.et al (2015). The FTMap Family of Web Servers for Determining and Characterizing Ligand-Binding Hot Spots of Proteins. Nat. Protoc.10 (5), 733–755. 10.1038/nprot.2015.043

  • CrossRef
  • Google Scholar

• 83

  KruseA. C.HuJ.PanA. C.ArlowD. H.RosenbaumD. M.RosemondE.et al (2012). Structure and Dynamics of the M3 Muscarinic Acetylcholine Receptor. Nature482 (7386), 552–556. 10.1038/nature10867

  • CrossRef
  • Google Scholar

• 84

  KruseA. C.RingA. M.ManglikA.HuJ.HuK.EitelK.et al (2013). Activation and Allosteric Modulation of a Muscarinic Acetylcholine Receptor. Nature504 (7478), 101–106. 10.1038/nature12735

  • CrossRef
  • Google Scholar

• 85

  KumarV.BonifaziA.EllenbergerM. P.KeckT. M.PommierE.RaisR.et al (2016). Highly Selective Dopamine D3 Receptor (D3R) Antagonists and Partial Agonists Based on Eticlopride and the D3R Crystal Structure: New Leads for Opioid Dependence Treatment. J. Med. Chem.59 (16), 7634–7650. 10.1021/acs.jmedchem.6b00860

  • CrossRef
  • Google Scholar

• 86

  Lamim RibeiroJ. M.ProvasiD.FilizolaM. (2020). A Combination of Machine Learning and Infrequent Metadynamics to Efficiently Predict Kinetic Rates, Transition States, and Molecular Determinants of Drug Dissociation from G Protein-Coupled Receptors. J. Chem. Phys.153 (12), 124105. 10.1063/5.0019100

  • CrossRef
  • Google Scholar

• 87

  LaneJ. R.SextonP. M.ChristopoulosA. (2013). Bridging the Gap: Bitopic Ligands of G-Protein-Coupled Receptors. Trends Pharmacol. Sci.34 (1), 59–66. 10.1016/j.tips.2012.10.003

  • CrossRef
  • Google Scholar

• 88

  Le GuillouxV.SchmidtkeP.TufferyP. (2009). Fpocket: An Open Source Platform for Ligand Pocket Detection. BMC Bioinformatics10, 168. 10.1186/1471-2105-10-168

  • CrossRef
  • Google Scholar

• 89

  LeeB.TaylorM.GriffinS. A.McInnisT.SumienN.MachR. H.et al (2021). Evaluation of Substituted N-Phenylpiperazine Analogs as D3 vs. D2 Dopamine Receptor Subtype Selective Ligands. Molecules26 (11), 3182. 10.3390/molecules26113182

  • CrossRef
  • Google Scholar

• 90

  LeeY.WarneT.NehméR.PandeyS.Dwivedi-AgnihotriH.ChaturvediM.et al (2020). Molecular Basis of β-arrestin Coupling to Formoterol-Bound β1-adrenoceptor. Nature583 (7818), 862–866. 10.1038/s41586-020-2419-1

  • CrossRef
  • Google Scholar

• 91

  LeoL. M.AboodM. E. (2021). CB1 Cannabinoid Receptor Signaling and Biased Signaling. Molecules26 (17), 5413. 10.3390/molecules26175413

  • CrossRef
  • Google Scholar

• 92

  LindenJ.ThaiT.FiglerH.JinX.RobevaA. S. (1999). Characterization of Human A(2B) Adenosine Receptors: Radioligand Binding, Western Blotting, and Coupling to G(q) in Human Embryonic Kidney 293 Cells and HMC-1 Mast Cells. Mol. Pharmacol.56 (4), 705–713.

  • Google Scholar

• 93

  LiuJ. J.HorstR.KatritchV.StevensR. C.WüthrichK. (2012). Biased Signaling Pathways in β2-adrenergic Receptor Characterized by 19F-NMR. Science335 (6072), 1106–1110. 10.1126/science.1215802

  • CrossRef
  • Google Scholar

• 94

  LiuW.WackerD.GatiC.HanG. W.JamesD.WangD.et al (2013). Serial Femtosecond Crystallography of G Protein-Coupled Receptors. Science342 (6165), 1521–1524. 10.1126/science.1244142

  • CrossRef
  • Google Scholar

• 95

  LuS.ZhangJ. (2017). Designed Covalent Allosteric Modulators: An Emerging Paradigm in Drug Discovery. Drug Discov. Today22 (2), 447–453. 10.1016/j.drudis.2016.11.013

  • CrossRef
  • Google Scholar

• 96

  LuY.LiuH.YangD.ZhongL.XinY.ZhaoS.et al (2021). Affinity Mass Spectrometry-Based Fragment Screening Identified a New Negative Allosteric Modulator of the Adenosine A2A Receptor Targeting the Sodium Ion Pocket. ACS Chem. Biol.16 (6), 991–1002. 10.1021/acschembio.0c00899

  • CrossRef
  • Google Scholar

• 97

  LuW.KosticM.ZhangT.CheJ.PatricelliM. P.ChouchaniE. T.et al (2021). Fragment-Based Covalent Ligand Discovery. RSC Chem. Biol.2 (2), 354–367. 10.1039/D0CB00222D

  • CrossRef
  • Google Scholar

• 98

  LudermanK. D.ConroyJ. L.FreeR. B.SouthallN.FerrerM.Sanchez-SotoM.et al (2018). Identification of Positive Allosteric Modulators of the D1 Dopamine Receptor that Act at Diverse Binding Sites. Mol. Pharmacol.94 (4), 1197–1209. 10.1124/mol.118.113175

  • CrossRef
  • Google Scholar

• 99

  LudermanK. D.JainP.Benjamin FreeR.ConroyJ. L.AubéJ.SibleyD. R.et al (2021). Development of Pyrimidone D1 Dopamine Receptor Positive Allosteric Modulators. Bioorg. Med. Chem. Lett.31, 127696. 10.1016/j.bmcl.2020.127696

  • CrossRef
  • Google Scholar

• 100

  MagarkarA.SchnappG.ApelA. K.SeeligerD.TautermannC. S. (2019). Enhancing Drug Residence Time by Shielding of Intra-protein Hydrogen Bonds: A Case Study on CCR2 Antagonists. ACS Med. Chem. Lett.10 (3), 324–328. 10.1021/acsmedchemlett.8b00590

  • CrossRef
  • Google Scholar

• 101

  MaoQ.QinW. Z.ZhangA.YeN. (2020). Recent Advances in Dopaminergic Strategies for the Treatment of Parkinson's Disease. Acta Pharmacol. Sin.41 (4), 471–482. 10.1038/s41401-020-0365-y

  • CrossRef
  • Google Scholar

• 102

  MasureelM.ZouY.PicardL. P.van der WesthuizenE.MahoneyJ. P.RodriguesJ. P. G. L. M.et al (2018). Structural Insights into Binding Specificity, Efficacy and Bias of a β2AR Partial Agonist. Nat. Chem. Biol.14 (11), 1059–1066. 10.1038/s41589-018-0145-x

  • CrossRef
  • Google Scholar

• 103

  McCorvyJ. D.ButlerK. V.KellyB.RechsteinerK.KarpiakJ.BetzR. M.et al (2018). Structure-inspired Design of β-arrestin-biased Ligands for Aminergic GPCRs. Nat. Chem. Biol.14 (2), 126–134. 10.1038/nchembio.2527

  • CrossRef
  • Google Scholar

• 104

  McNeillS. M.BaltosJ. A.WhiteP. J.MayL. T. (2021). Biased Agonism at Adenosine Receptors. Cell. Signal.82, 109954. 10.1016/j.cellsig.2021.109954

  • CrossRef
  • Google Scholar

• 105

  MelchiorreC.CassinelliA.QuagliaW. (1987). Differential Blockade of Muscarinic Receptor Subtypes by Polymethylene Tetraamines. Novel Class of Selective Antagonists of Cardiac M-2 Muscarinic Receptors. J. Med. Chem.30 (1), 201–204. 10.1021/jm00384a034

  • CrossRef
  • Google Scholar

• 106

  MiaoY.BhattaraiA.NguyenA. T. N.ChristopoulosA.MayL. T. (2018). Structural Basis for Binding of Allosteric Drug Leads in the Adenosine A1 Receptor. Sci. Rep.8 (1), 16836. 10.1038/s41598-018-35266-x

  • CrossRef
  • Google Scholar

• 107

  MiaoY.GoldfeldD. A.MooE. V.SextonP. M.ChristopoulosA.McCammonJ. A.et al (2016). Accelerated Structure-Based Design of Chemically Diverse Allosteric Modulators of a Muscarinic G Protein-Coupled Receptor. Proc. Natl. Acad. Sci. U S A.113 (38), E5675–E5684. 10.1073/pnas.1612353113

  • CrossRef
  • Google Scholar

• 108

  MichinoM.BoatengC. A.DonthamsettiP.YanoH.BakareO. M.BonifaziA.et al (2017). Toward Understanding the Structural Basis of Partial Agonism at the Dopamine D3 Receptor. J. Med. Chem.60 (2), 580–593. 10.1021/acs.jmedchem.6b01148

  • CrossRef
  • Google Scholar

• 109

  MielnikC. A.LamV. M.RossR. A. (2021). CB1 Allosteric Modulators and Their Therapeutic Potential in CNS Disorders. Prog. Neuropsychopharmacol. Biol. Psychiatry106, 110163. 10.1016/j.pnpbp.2020.110163

  • CrossRef
  • Google Scholar

• 110

  NagiriC.KobayashiK.TomitaA.KatoM.KobayashiK.YamashitaK.et al (2021). Cryo-EM Structure of the β3-adrenergic Receptor Reveals the Molecular Basis of Subtype Selectivity. Mol. Cel81 (15), 3205–e5. 10.1016/j.molcel.2021.06.024

  • CrossRef
  • Google Scholar

• 111

  NewmanA. H.BattitiF. O.BonifaziA. (20162020). 2016 Philip S. Portoghese Medicinal Chemistry Lectureship: Designing Bivalent or Bitopic Molecules for G-Protein Coupled Receptors. The Whole Is Greater Than the Sum of its Parts. J. Med. Chem.63 (5), 1779–1797. 10.1021/acs.jmedchem.9b01105

  • CrossRef
  • Google Scholar

• 112

  NewmanA. H.BeumingT.BanalaA. K.DonthamsettiP.PongettiK.LaBountyA.et al (2012). Molecular Determinants of Selectivity and Efficacy at the Dopamine D3 Receptor. J. Med. Chem.55 (15), 6689–6699. 10.1021/jm300482h

  • CrossRef
  • Google Scholar

• 113

  OrgovánZ.FerenczyG. G.KeserűG. M. (2019). Fragment-Based Approaches for Allosteric Metabotropic Glutamate Receptor (MGluR) Modulators. Ctmc19 (19), 1768–1781. 10.2174/1568026619666190808150039

  • CrossRef
  • Google Scholar

• 114

  OveringtonJ. P.Al-LazikaniB.HopkinsA. L. (2006). How Many Drug Targets Are There. Nat. Rev. Drug Discov.5 (12), 993–996. 10.1038/nrd2199

  • CrossRef
  • Google Scholar

• 115

  PedersenM. F.WróbelT. M.Märcher-RørstedE.PedersenD. S.MøllerT. C.GabrieleF.et al (2020). Biased Agonism of Clinically Approved μ-opioid Receptor Agonists and TRV130 Is Not Controlled by Binding and Signaling Kinetics. Neuropharmacology166, 107718. 10.1016/j.neuropharm.2019.107718

  • CrossRef
  • Google Scholar

• 116

  PengY.McCorvyJ. D.HarpsøeK.LansuK.YuanS.PopovP.et al (2018). 5-HT2C Receptor Structures Reveal the Structural Basis of GPCR Polypharmacology. Cell172 (4), 719–e14. 10.1016/j.cell.2018.01.001

  • CrossRef
  • Google Scholar

• 117

  PottertonA.HusseiniF. S.SoutheyM. W. Y.BodkinM. J.HeifetzA.CoveneyP. V.et al (2019). Ensemble-Based Steered Molecular Dynamics Predicts Relative Residence Time of A2A Receptor Binders. J. Chem. Theor. Comput.15 (5), 3316–3330. 10.1021/acs.jctc.8b01270

  • CrossRef
  • Google Scholar

• 118

  PottertonA.HeifetzA.Townsend-NicholsonA. (2022). Predicting Residence of GPCR Ligands with Machine Learning. Methods Mol. Biol. Clifton NJ2390, 191–205. 10.1007/978-1-0716-1787-8_8

  • CrossRef
  • Google Scholar

• 119

  ProkopS.Ábrányi-BaloghP.BartiB.VámosiM.ZöldiM.BarnaL.et al (2021). PharmacoSTORM Nanoscale Pharmacology Reveals Cariprazine Binding on Islands of Calleja Granule Cells. Nat. Commun.12 (1), 6505. 10.1038/s41467-021-26757-z

  • CrossRef
  • Google Scholar

• 120

  RafaelF.JosemaC.EnricI. C. (2020). The Kinetic Component in Drug Discovery: Using the Most Basic Pharmacological Concepts to Advance in Selecting Drugs to Combat CNS Diseases. Curr. Neuropharmacol.18 (3), 250–257. 10.2174/1570159X17666191001144309

  • CrossRef
  • Google Scholar

• 121

  RaschkaS.KaufmanB. (2020). Machine Learning and AI-Based Approaches for Bioactive Ligand Discovery and GPCR-Ligand Recognition. Methods180, 89–110. 10.1016/j.ymeth.2020.06.016

  • CrossRef
  • Google Scholar

• 122

  RiddyD. M.CookA. E.ShacklefordD. M.PierceT. L.MocaerE.Mannoury la CourC.et al (2019). Drug-Receptor Kinetics and Sigma-1 Receptor Affinity Differentiate Clinically Evaluated Histamine H3 Receptor Antagonists. Neuropharmacology144, 244–255. 10.1016/j.neuropharm.2018.10.028

  • CrossRef
  • Google Scholar

• 123

  RothB. L.ShefflerD. J.KroezeW. K. (2004). Magic Shotguns versus Magic Bullets: Selectively Non-selective Drugs for Mood Disorders and Schizophrenia. Nat. Rev. Drug Discov.3 (4), 353–359. 10.1038/nrd1346

  • CrossRef
  • Google Scholar

• 124

  SalmasoV.JacobsonK. A. (2020). In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors. Biomolecules10 (6), 812. 10.3390/biom10060812

  • CrossRef
  • Google Scholar

• 125

  SarkarP.ChattopadhyayA. (2020). Cholesterol Interaction Motifs in G Protein-Coupled Receptors: Slippery Hot Spots. Wiley Interdiscip. Rev. Syst. Biol. Med.12 (4), e1481. 10.1002/wsbm.1481

  • CrossRef
  • Google Scholar

• 126

  SchmidtkeP.Le GuillouxV.MaupetitJ.TufféryP. (2010). Fpocket: Online Tools for Protein Ensemble Pocket Detection and Tracking. Nucleic Acids Res.38, W582–W589. 10.1093/nar/gkq383

  • CrossRef
  • Google Scholar

• 127

  SchrammS.AgnettaL.BermudezM.GerweH.IrmenM.HolzeJ.et al (2019). Novel BQCA- and TBPB-Derived M1 Receptor Hybrid Ligands: Orthosteric Carbachol Differentially Regulates Partial Agonism. ChemMedChem14 (14), 1349–1358. 10.1002/cmdc.201900283

  • CrossRef
  • Google Scholar

• 128

  ShaikA. B.BoatengC. A.BattitiF. O.BonifaziA.CaoJ.ChenL.et al (2021). Structure Activity Relationships for a Series of Eticlopride-Based Dopamine D2/D3 Receptor Bitopic Ligands. J. Med. Chem.64 (20), 15313–15333. 10.1021/acs.jmedchem.1c01353

  • CrossRef
  • Google Scholar

• 129

  ShimamuraT.ShiroishiM.WeyandS.TsujimotoH.WinterG.KatritchV.et al (2011). Structure of the Human Histamine H1 Receptor Complex with Doxepin. Nature475 (7354), 65–70. 10.1038/nature10236

  • CrossRef
  • Google Scholar

• 130

  SriramK.InselP. A. (2018). G Protein-Coupled Receptors as Targets for Approved Drugs: How Many Targets and How Many Drugs. Mol. Pharmacol.93 (4), 251–258. 10.1124/mol.117.111062

  • CrossRef
  • Google Scholar

• 131

  StausD. P.HuH.RobertsonM. J.KleinhenzA. L. W.WinglerL. M.CapelW. D.et al (2020). Structure of the M2 Muscarinic Receptor-β-Arrestin Complex in a Lipid Nanodisc. Nature579 (7798), 297–302. 10.1038/s41586-020-1954-0

  • CrossRef
  • Google Scholar

• 132

  SuM.ZhuL.ZhangY.PaknejadN.DeyR.HuangJ.et al (2020). Structural Basis of the Activation of Heterotrimeric Gs-Protein by Isoproterenol-Bound β1-Adrenergic Receptor. Mol. Cel80 (1), 59–e4. 10.1016/j.molcel.2020.08.001

  • CrossRef
  • Google Scholar

• 133

  SykesD. A.StoddartL. A.KilpatrickL. E.HillS. J. (2019). Binding Kinetics of Ligands Acting at GPCRs. Mol. Cel. Endocrinol.485, 9–19. 10.1016/j.mce.2019.01.018

  • CrossRef
  • Google Scholar

• 134

  SzlenkC. T.GcJ. B.NatesanS. (2019). Does the Lipid Bilayer Orchestrate Access and Binding of Ligands to Transmembrane Orthosteric/Allosteric Sites of G Protein-Coupled Receptors. Mol. Pharmacol.96 (5), 527–541. 10.1124/mol.118.115113

  • CrossRef
  • Google Scholar

• 135

  TanL.ZhouQ.YanW.SunJ.KozikowskiA. P.ZhaoS.et al (2020). Design and Synthesis of Bitopic 2-Phenylcyclopropylmethylamine (PCPMA) Derivatives as Selective Dopamine D3 Receptor Ligands. J. Med. Chem.63 (9), 4579–4602. 10.1021/acs.jmedchem.9b01835

  • CrossRef
  • Google Scholar

• 136

  TautermannC. S.KiechleT.SeeligerD.DiehlS.WexE.BanholzerR.et al (2013). Molecular Basis for the Long Duration of Action and Kinetic Selectivity of Tiotropium for the Muscarinic M3 Receptor. J. Med. Chem.56 (21), 8746–8756. 10.1021/jm401219y

  • CrossRef
  • Google Scholar

• 137

  ThalD. M.SunB.FengD.NawaratneV.LeachK.FelderC. C.et al (2016). Crystal Structures of the M1 and M4 Muscarinic Acetylcholine Receptors. Nature531 (7594), 335–340. 10.1038/nature17188

  • CrossRef
  • Google Scholar

• 138

  ThomasT.FangY.YurievE.ChalmersD. K. (2016). Ligand Binding Pathways of Clozapine and Haloperidol in the Dopamine D2 and D3 Receptors. J. Chem. Inf. Model.56 (2), 308–321. 10.1021/acs.jcim.5b00457

  • CrossRef
  • Google Scholar

• 139

  ValantC.Robert LaneJ.SextonP. M.ChristopoulosA. (2012). The Best of Both Worlds? Bitopic Orthosteric/allosteric Ligands of G Protein-Coupled Receptors. Annu. Rev. Pharmacol. Toxicol.52 (1), 153–178. 10.1146/annurev-pharmtox-010611-134514

  • CrossRef
  • Google Scholar

• 140

  van der VeldenW. J. C.HeitmanL. H.RosenkildeM. M. (2020). Perspective: Implications of Ligand-Receptor Binding Kinetics for Therapeutic Targeting of G Protein-Coupled Receptors. ACS Pharmacol. Transl. Sci.3 (2), 179–189. 10.1021/acsptsci.0c00012

  • CrossRef
  • Google Scholar

• 141

  VassM.Agai-CsongorE.HortiF.KeserűG. M. (2014). Multiple Fragment Docking and Linking in Primary and Secondary Pockets of Dopamine Receptors. ACS Med. Chem. Lett.5 (9), 1010–1014. 10.1021/ml500201u

  • CrossRef
  • Google Scholar

• 142

  VermaR. K.AbramyanA. M.MichinoM.FreeR. B.SibleyD. R.JavitchJ. A.et al (2018). The E2.65A Mutation Disrupts Dynamic Binding Poses of SB269652 at the Dopamine D2 and D3 Receptors. PLOS Comput. Biol.14 (1), e1005948. 10.1371/journal.pcbi.1005948

  • CrossRef
  • Google Scholar

• 143

  VuckovicZ.GentryP. R.BerizziA. E.HirataK.VargheseS.ThompsonG.et al (2019). Crystal Structure of the M5 Muscarinic Acetylcholine Receptor. Proc. Natl. Acad. Sci. U. S. A.116 (51), 26001–26007. 10.1073/pnas.1914446116

  • CrossRef
  • Google Scholar

• 144

  WackerD.WangC.KatritchV.HanG. W.HuangX. P.VardyE.et al (2013). Structural Features for Functional Selectivity at Serotonin Receptors. Science340 (6132), 615–619. 10.1126/science.1232808

  • CrossRef
  • Google Scholar

• 145

  WackerD.WangS.McCorvyJ. D.BetzR. M.VenkatakrishnanA. J.LevitA.et al (2017). Crystal Structure of an LSD-Bound Human Serotonin Receptor. Cell168 (3), 377–e12. 10.1016/j.cell.2016.12.033

  • CrossRef
  • Google Scholar

• 146

  WakefieldA. E.MasonJ. S.VajdaS.KeserűG. M. (2019). Analysis of Tractable Allosteric Sites in G Protein-Coupled Receptors. Sci. Rep.9 (1), 6180. 10.1038/s41598-019-42618-8

  • CrossRef
  • Google Scholar

• 147

  WangC.JiangY.MaJ.WuH.WackerD.KatritchV.et al (2013). Structural Basis for Molecular Recognition at Serotonin Receptors. Science340 (6132), 610–614. 10.1126/science.1232807

  • CrossRef
  • Google Scholar

• 148

  WangH.ReineckeB. A.ZhangY. (2020). Computational Insights into the Molecular Mechanisms of Differentiated Allosteric Modulation at the Mu Opioid Receptor by Structurally Similar Bitopic Modulators. J. Comput. Aided Mol. Des.34 (8), 879–895. 10.1007/s10822-020-00309-x

  • CrossRef
  • Google Scholar

• 149

  WangY.YuZ.XiaoW.LuS.ZhangJ. (2021). Allosteric Binding Sites at the Receptor-Lipid Bilayer Interface: Novel Targets for GPCR Drug Discovery. Drug Discov. Today26 (3), 690–703. 10.1016/j.drudis.2020.12.001

  • CrossRef
  • Google Scholar

• 150

  WhiteJ. F.NoinajN.ShibataY.LoveJ.KlossB.XuF.et al (2012). Structure of the Agonist-Bound Neurotensin Receptor. Nature490 (7421), 508–513. 10.1038/nature11558

  • CrossRef
  • Google Scholar

• 151

  WoldE. A.ChenJ.CunninghamK. A.ZhouJ. (2019). Allosteric Modulation of Class A GPCRs: Targets, Agents, and Emerging Concepts. J. Med. Chem.62 (1), 88–127. 10.1021/acs.jmedchem.8b00875

  • CrossRef
  • Google Scholar

• 152

  WuY.ZengL.ZhaoS. (2021). Ligands of Adrenergic Receptors: A Structural Point of View. Biomolecules11 (7), 936. 10.3390/biom11070936

  • CrossRef
  • Google Scholar

• 153

  XiaR.WangN.XuZ.LuY.SongJ.ZhangA.et al (2021). Cryo-EM Structure of the Human Histamine H1 Receptor/Gq Complex. Nat. Commun.12 (1), 2086. 10.1038/s41467-021-22427-2

  • CrossRef
  • Google Scholar

• 154

  XiaoP.YanW.GouL.ZhongY. N.KongL.WuC.et al (2021). Ligand Recognition and Allosteric Regulation of DRD1-Gs Signaling Complexes. Cell184 (4), 943–e18. 10.1016/j.cell.2021.01.028

  • CrossRef
  • Google Scholar

• 155

  XuP.HuangS.ZhangH.MaoC.ZhouX. E.ChengX.et al (2021a). Structural Insights into the Lipid and Ligand Regulation of Serotonin Receptors. Nature592 (7854), 469–473. 10.1038/s41586-021-03376-8

  • CrossRef
  • Google Scholar

• 156

  XuP.HuangS.MaoC.KrummB. E.ZhouX. E.TanY.et al (2021b). Structures of the Human Dopamine D3 Receptor-Gi Complexes. Mol. Cel81 (6), 1147–e4. 10.1016/j.molcel.2021.01.003

  • CrossRef
  • Google Scholar

• 157

  XuX.KaindlJ.ClarkM. J.HübnerH.HirataK.SunaharaR. K.et al (2021). Binding Pathway Determines Norepinephrine Selectivity for the Human β1AR over β2AR. Cel Res31 (5), 569–579. 10.1038/s41422-020-00424-2

  • CrossRef
  • Google Scholar

• 158

  YanW.FanL.YuJ.LiuR.WangH.TanL.et al (2021). 2-Phenylcyclopropylmethylamine Derivatives as Dopamine D2 Receptor Partial Agonists: Design, Synthesis, and Biological Evaluation. J. Med. Chem.64 (23), 17239–17258. 10.1021/acs.jmedchem.1c01327

  • CrossRef
  • Google Scholar

• 159

  YangD.ZhouQ.LabroskaV.QinS.DarbalaeiS.WuY.et al (2021). G Protein-Coupled Receptors: Structure- and Function-Based Drug Discovery. Sig Transduct Target. Ther.6 (1), 1–27. 10.1038/s41392-020-00435-w

  • CrossRef
  • Google Scholar

• 160

  YangF.LingS.ZhouY.ZhangY.LvP.LiuS.et al (2021). Different Conformational Responses of the β2-adrenergic Receptor-Gs Complex upon Binding of the Partial Agonist Salbutamol or the Full Agonist Isoprenaline. Natl. Sci. Rev.8 (9), nwaa284. 10.1093/nsr/nwaa284

  • CrossRef
  • Google Scholar

• 161

  YinJ.ChenK. M.ClarkM. J.HijaziM.KumariP.BaiX. C.et al (2020). Structure of a D2 Dopamine Receptor-G-Protein Complex in a Lipid Membrane. Nature584 (7819), 125–129. 10.1038/s41586-020-2379-5

  • CrossRef
  • Google Scholar

• 162

  YuanD.LiuZ.KaindlJ.MaedaS.ZhaoJ.SunX.et al (2020). Activation of the α2B Adrenoceptor by the Sedative Sympatholytic Dexmedetomidine. Nat. Chem. Biol.16 (5), 507–512. 10.1038/s41589-020-0492-2

  • CrossRef
  • Google Scholar

• 163

  ZhangG.ChengJ.McCorvyJ. D.LorelloP. J.CaldaroneB. J.RothB. L.et al (2017). Discovery of N-Substituted (2-Phenylcyclopropyl)Methylamines as Functionally Selective Serotonin 2C Receptor Agonists for Potential Use as Antipsychotic Medications. J. Med. Chem.60 (14), 6273–6288. 10.1021/acs.jmedchem.7b00584

  • CrossRef
  • Google Scholar

• 164

  ZhangY.YangF.LingS.LvP.ZhouY.FangW.et al (2020). Single-particle Cryo-EM Structural Studies of the β2AR-Gs Complex Bound with a Full Agonist Formoterol. Cell Discov6 (1), 45–5. 10.1038/s41421-020-0176-9

  • CrossRef
  • Google Scholar

• 165

  ZhuangY.XuP.MaoC.WangL.KrummB.ZhouX. E.et al (2021a). Structural Insights into the Human D1 and D2 Dopamine Receptor Signaling Complexes. Cell184 (4), 931–e18. e18. 10.1016/j.cell.2021.01.027

  • CrossRef
  • Google Scholar

• 166

  ZhuangY.KrummB.ZhangH.ZhouX. E.WangY.HuangX.-P.et al (2021b). Mechanism of Dopamine Binding and Allosteric Modulation of the Human D1 Dopamine Receptor. Cel Res31 (5), 593–596. 10.1038/s41422-021-00482-0

  • CrossRef
  • Google Scholar