Rotamer Toggle Research Articles

A Novel "Activation Switch" Motif Common to All Aminergic Receptors.

Aminergic receptors are G protein-coupled receptors (GPCRs) that transduce signals from small endogenous biogenic amines to regulate intracellular signaling pathways. Agonist binding in the ligand binding pocket on the extracellular side opens and prepares a cavity on the intracellular face of the receptors to interact with and activate G proteins and β-arrestins. Here, by reviewing and analyzing all available aminergic receptor structures, we seek to identify activation-related conformational changes that are independent of the specific scaffold of the bound agonist, which we define as "activation conformational changes" (ACCs). While some common intracellular ACCs have been well-documented, identifying common extracellular ACCs, including those in the ligand binding pocket, is complicated by local adjustments to different ligand scaffolds. Our analysis shows no common ACCs at the extracellular ends of the transmembrane helices. Furthermore, the restricted access to the ligand binding pocket identified previously in some receptors is not universal. Notably, the Trp6.48 toggle switch and the Pro5.50-Ile3.40-Phe6.44 (PIF) motif at the bottom of the ligand binding pocket have previously been proposed to mediate the conformational consequences of ligand binding to the intracellular side of the receptors. Our analysis shows that common ACCs in the ligand binding pocket are associated with the PIF motif and nearby residues, including Trp6.48, but fails to support a shared rotamer toggle associated with activation. However, we identify two common rearrangements between the extracellular and middle subsegments, and propose a novel "activation switch" motif common to all aminergic receptors. This motif includes the middle subsegments of transmembrane helices 3, 5, and 6 and integrates both the PIF motif and Trp6.48.

Different conformational responses of the β2-adrenergic receptor-Gs complex upon binding of the partial agonist salbutamol or the full agonist isoprenaline.

G protein-coupled receptors (GPCRs) are responsible for most cytoplasmic signaling in response to extracellular ligands with different efficacy profiles. Various spectroscopic techniques have identified that agonists exhibiting varying efficacies can selectively stabilize a specific conformation of the receptor. However, the structural basis for activation of the GPCR-G protein complex by ligands with different efficacies is incompletely understood. To better understand the structural basis underlying the mechanisms by which ligands with varying efficacies differentially regulate the conformations of receptors and G proteins, we determined the structures of β2AR-Gαs[Formula: see text]γ bound with partial agonist salbutamol or bound with full agonist isoprenaline using single-particle cryo-electron microscopy at resolutions of 3.26Å and 3.80Å, respectively. Structural comparisons between the β2AR-Gs-salbutamol and β2AR-Gs-isoprenaline complexes demonstrated that the decreased binding affinity and efficacy of salbutamol compared with those of isoprenaline might be attributed to weakened hydrogen bonding interactions, attenuated hydrophobic interactions in the orthosteric binding pocket and different conformational changes in the rotamer toggle switch in TM6. Moreover, the observed stronger interactions between the intracellular loop 2 or 3 (ICL2 or ICL3) of β2AR and Gαs with binding of salbutamol versus isoprenaline might decrease phosphorylation in the salbutamol-activated β2AR-Gs complex. From the observed structural differences between these complexes of β2AR, a mechanism of β2AR activation by partial and full agonists is proposed to provide structural insights into β2AR desensitization.

Biased Agonist TRV027 Determinants in AT1R by Molecular Dynamics Simulations.

Functional selectivity is a phenomenon observed in G protein-coupled receptors in which intermediate active-state conformations are stabilized by mutations or ligand binding, resulting in different sets of signaling pathways. Peptides capable of selectively activating β-arrestin, known as biased agonists, have already been characterized in vivo and could correspond to a new therapeutic approach for treatment of cardiovascular diseases. Despite the potential of biased agonism, the mechanism involved in selective signaling remains unclear. In this work, molecular dynamics simulations were employed to compare the conformational profile of the angiotensin II type 1 receptor (AT1R) crystal bound to angiotensin II, bound to the biased ligand TRV027, and in the apo form. Our results show that both ligands induce changes near the NPxxY motif in transmembrane domain 7 that are related to receptor activation. However, the biased ligand does not cause the rotamer toggle alternative positioning and displays an exclusive hydrogen-bonding pattern. Our work sheds light on the biased agonism mechanism and will help in the future design of novel biased agonists for AT1R.

An insight into antagonist binding and induced conformational dynamics of class B GPCR corticotropin-releasing factor receptor 1.

The human corticotropin-releasing factor receptor type 1 (CRF1R) is a class B G-protein-coupled receptor (GPCR), which mediates the response to stress and has been considered as a drug target for depression and anxiety. Based on the CRF1R-antagonist crystal structure, we study the binding mechanism of two distinct antagonists, CP-376395 and MTIP, and the dynamics behaviors of CRF1R induced by an antagonist binding. Key residues interacting with both antagonists and residues specifically binding to one of them are identified. Both antagonists interact with Asn283, Phe203, Met206, Leu280, Tyr316, Leu323, Leu287, Phe284, Val279, Leu319, Phe207, Gly210 and Phe362. CP-376395 specifically binds to Glu209 and Phe160, while MTIP specifically binds to Leu320, Leu213, Ile290, Phe162 and Val313. The total binding free energy of MTIP is lower than that of CP-376395; this is consistent with the experimental observation that MTIP shows higher binding affinity than CP-376395. The conformational dynamic behaviors of antagonist bound holo-CRF1R were found to be different from those of apo-CRF1R in three aspects: (i) the "ionic lock" between side chains of Arg151 in TM2 and Glu209 in TM3 was broken in apo-CRF1R, but was formed in holo-CRF1Rs; (ii) Phe203 in TM3 and Tyr327 in TM6 were in close proximity to each other in apo-CRF1R, while they were far apart resulting from the shift of TM6 in holo-CRF1Rs; and (iii) the "rotamer toggle switch", Tyr327/Leu323/Phe284, adopted different rotameric conformations in apo-CRF1R and holo-CRF1Rs. We hope that our results could be helpful in further development of the drug design of CRF1R.

W246(6.48) opens a gate for a continuous intrinsic water pathway during activation of the adenosine A2A receptor.

The question how G-protein-coupled receptors transduce an extracellular signal by a sequence of transmembrane conformational transitions into an intracellular response remains to be solved at molecular detail. Herein, we use molecular dynamics simulations to reveal distinct conformational transitions of the adenosine A2A receptor, and we found that the conserved W246(6.48) residue in transmembrane helix TM6 performs a key rotamer toggle switch. Agonist binding induces the sidechain of W246(6.48) to fluctuate between two distinct conformations enabling the diffusion of water molecules from the bulk into the center of the receptor. After passing the W246(6.48) gate, the internal water molecules induce another conserved residue, Y288(7.53), to switch to a distinct rotamer conformation establishing a continuous transmembrane water pathway. Further, structural changes of TM6 and TM7 induce local structural changes of the adjacent lipid bilayer.

W2466.48 Opens a Gate for a Continuous Intrinsic Water Pathway during Activation of the Adenosine A2A Receptor

AbstractThe question how G‐protein‐coupled receptors transduce an extracellular signal by a sequence of transmembrane conformational transitions into an intracellular response remains to be solved at molecular detail. Herein, we use molecular dynamics simulations to reveal distinct conformational transitions of the adenosine A2A receptor, and we found that the conserved W2466.48 residue in transmembrane helix TM6 performs a key rotamer toggle switch. Agonist binding induces the sidechain of W2466.48 to fluctuate between two distinct conformations enabling the diffusion of water molecules from the bulk into the center of the receptor. After passing the W2466.48 gate, the internal water molecules induce another conserved residue, Y2887.53, to switch to a distinct rotamer conformation establishing a continuous transmembrane water pathway. Further, structural changes of TM6 and TM7 induce local structural changes of the adjacent lipid bilayer.

Long time molecular dynamic simulation on the agonist binding and activation of the β2-adrenergic receptor

β2-Adrenergic receptor (β2AR) plays a critical role in mediating the effects of catecholamine hormones. Due to the flexibility of the structure of its active state, study of agonist–β2AR is usually performed by molecular dynamic (MD) simulation. In this study we show the representative characteristics of agonist binding and activation on β2AR by MD simulation. The binding site and the conformational changes in the specific regions of β2AR are reasonable which confirmed the conclusion that agonist–β2AR reached its active-like state. We have studied the disruption of non-covalent intramolecular interactions, including the conserved DRY motif, the rotamer toggle switch and the ionic lock, the cytoplasmic ends of transmembranes 5 and 6, and some water-mediated hydrogen bond network regions. We conclude that agonist induced β2AR to its active conformation, or at least the active-like conformation.

GPCR A2AAR Agonist Binding and Induced Conformation Changes of Functional Switches

Agonist binding of A2A adenosine receptor (A2AAR) shows protective effects against inflammatory and immune. Efforts are exerted in understanding the general mechanism and developing A2AAR selectively binding agonists. Using molecular dynamics (MD) simulations, we have studied the interactions between A2AAR and its agonist (adenosine), and analyzed the induced dynamic behaviors of the receptor. Key residues interacting with adenosine are identified: A632.61, I662.64, V843.32, L853.33, T883.36, F1685.29, M1775.38, L2496.51, H2506.52, and N2536.55 interacting with adenosine with affinities larger than 0.5 kcal/mol. Moreover, no interaction between adenosine and L1675.28 is observed, which supports our previous findings that L1675.28 is an antagonist specific binding reside. The dynamic behaviors of agonist bound A2AAR are found to be different from apo-A2AAR in three typical functional switches: (i) tight “ionic lock” forms in adenosine-A2AAR, but it is in equilibrium between formation and breakage in apo-A2AAR; (ii) the “rotamer toggle switch”, T883.36/F2426.44/W2466.48, adopted different rotameric conformations in adenosine-A2AAR and apo-A2AAR; (iii) adenosine-A2AAR has a flexible intracellular loop 2 (IC2) and α-helical IC3, while apo-A2AAR preferred α-helical IC2 and flexible IC3. Our results indicate that agonist binding induced different conformational rearrangements of these characteristic functional switches in adenosine-A2AAR and apo-A2AAR.

Lipid Receptor S1P1 Activation Scheme Concluded from Microsecond All-Atom Molecular Dynamics Simulations

Sphingosine 1-phosphate (S1P) is a lysophospholipid mediator which activates G protein–coupled sphingosine 1-phosphate receptors and thus evokes a variety of cell and tissue responses including lymphocyte trafficking, endothelial development, integrity, and maturation. We performed five all-atom 700 ns molecular dynamics simulations of the sphingosine 1-phosphate receptor 1 (S1P1) based on recently released crystal structure of that receptor with an antagonist. We found that the initial movements of amino acid residues occurred in the area of highly conserved W2696.48 in TM6 which is close to the ligand binding location. Those residues located in the central part of the receptor and adjacent to kinks of TM helices comprise of a transmission switch. Side chains movements of those residues were coupled to the movements of water molecules inside the receptor which helped in the gradual opening of intracellular part of the receptor. The most stable parts of the protein were helices TM1 and TM2, while the largest movement was observed for TM7, possibly due to the short intracellular part starting with a helix kink at P7.50, which might be the first helix to move at the intracellular side. We show for the first time the detailed view of the concerted action of the transmission switch and Trp (W6.48) rotamer toggle switch leading to redirection of water molecules flow in the central part of the receptor. That event is a prerequisite for subsequent changes in intracellular part of the receptor involving water influx and opening of the receptor structure.

Antagonist binding and induced conformational dynamics of GPCR A2A adenosine receptor

The A2A adenosine receptor (A2AAR) is a unique G-protein coupled receptor (GPCR), because besides agonist, its antagonist could also lead to therapeutic relevance. Based on A2AAR-antagonist crystal structure, we have studied the binding mechanism of two distinct antagonists, ZM241385 and KW6002, and dynamic behaviors of A2AAR induced by antagonist binding. Key residues interacting with both antagonists and residues specifically binding to one of them are identified. ZM241385 specifically bound to S67(2.65), M177(5.38), and N253(6.55), while KW6002 binds to F62(2.60), A81(3.29), and H264(7.29). Moreover, interactions with L167(5.28) are found for both antagonists, which were not reported in agonist binding. The dynamic behaviors of antagonist bound holo-A2AARs were found to be different from the apo-A2AAR in three typical functional switches, (i) the "ionic lock" was in equilibrium between formation and breakage in apo-A2AAR, but stayed broken in holo-A2AARs; (ii) the "rotamer toggle switch," T88(3.36)/F242(6.44)/W246(6.48), adopted different rotameric conformations in apo-A2AAR and holo-A2AARs; (iii) apo-A2AAR preferred α-helical intracellular loop (IC)2 and flexible IC3, while holo-A2AARs had a flexible IC2 and α-helical IC3. Our results indicated that antagonist binding induced different conformational rearrangements of these characteristic functional switches in apo-A2AAR and holo-A2AARs.

Ligand-Perturbed Allosteric Communication within the Human A2A Adenosine Receptor

The human A2A adenosine receptor (A2AR) is a member of Class A G protein-coupled receptors (GPCRs). Activation of A2A adenosine receptor leads to the stimulation of adenylyl cyclase (AC), and thereby increases the formation of cyclic adenosine monophosphate (cAMP). It has been recognized that A2A adenosine receptor is a potential therapeutic target for the treatment of Parkinson's disease, Huntington's disease, schizophrenia, and other neurodegenerative diseases. In 2008, the first human A2A adenosine receptor was crystalized in complex with its high affinity selective antagonist ZM241385. In 2011, six more A2A adenosine receptor structures were published. Three of them co-crystallized with its agonists (UK432097, adenosine, and NECA) and three of the others co-crystallized with its antagonists (ZM241385, xanthine amine congener, and caffeine). In this study, the changes of allosteric communication inside the A2A adenosine receptor induced by these agonists and antagonists are examined by molecular dynamics simulations. We found that agonists are able to switch the torsion angle of W2466.48, and the correlated motions within the A2A adenosine receptor are also perturbed by this rotamer toggle switch. In addition to conventional dynamical cross-correlation matrix analysis, we also investigated the time-lagged correlated motions to see the persistence time of autocorrelation and cross-correlation, and to examine the causality relationships between different dynamics states of the receptor.1. Huang, N.-K.; Lin, J.-H. ….Chern Y. A new drug design targeting the adenosinergic system for Huntington's disease. PLoS. ONE 2011, 6, e20934.2. Chen, J.-B.;. ….Lin, J.-H., Fang J.-M. Design and synthesis of novel dual-action compounds targeting adenosine A2A receptor for neuroprotection. ChemMedChem 2011, 6, 1390-1400.

G-protein-coupled receptor dynamics: dimerization and activation models compared with experiment

Our previously derived models of the active state of the β2-adrenergic receptor are compared with recently published X-ray crystallographic structures of activated GPCRs (G-protein-coupled receptors). These molecular dynamics-based models using experimental data derived from biophysical experiments on activation were used to restrain the receptor to an active state that gave high enrichment for agonists in virtual screening. The β2-adrenergic receptor active model and X-ray structures are in good agreement over both the transmembrane region and the orthosteric binding site, although in some regions the active model is more similar to the active rhodopsin X-ray structures. The general features of the microswitches were well reproduced, but with minor differences, partly because of the unexpected X-ray results for the rotamer toggle switch. In addition, most of the interacting residues between the receptor and the G-protein were identified. This analysis of the modelling has also given important additional insight into GPCR dimerization: re-analysis of results on photoaffinity analogues of rhodopsin provided additional evidence that TM4 (transmembrane helix 4) resides at the dimer interface and that ligands such as bivalent ligands may pass between the mobile helices. A comparison, and discussion, is also carried out between the use of implicit and explicit solvent for active-state modelling.

The role of helices 5 and 6 on the human β1-adrenoceptor activation mechanism

The human β1-adrenoceptor (β1AR) is a G-protein-coupled receptor (GPCR) involved in sympathetic system regulation through agonist-induced activation. The conserved CWXP-motif in helix 6 (rotamer toggle switch) is one of the most important activation switches in Class A GPCRs. In order to investigate how the agonist binding disturbs this switch, we carried out molecular dynamics simulations of a hβ1AR model in the apo and R-noradrenaline-bound forms. The results show that the agonist binding changes the β1-angle distribution of Cys336, Trp337 and Phe341 residues and increases the helix 6 bending. Overall, we provide a functional hβ1AR model, showing how the rotamer toggle switch mechanism works at atomic level.

Action of Molecular Switches in GPCRs - Theoretical and Experimental Studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called “molecular switches” buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homo- and heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Structural insights into human GPCR protein OA1: a computational perspective

Human ocular albinism type 1 protein (OA1)-a member of the G-protein coupled receptor (GPCR) superfamily-is an integral membrane glycoprotein expressed exclusively by intracellular organelles known as melanocytes, and is responsible for the proper biogenesis of melanosomes. Mutations in the Oa1 gene are responsible for the disease ocular albinism. Despite its clinical importance, there is a lack of in-depth understanding of its structure and mechanism of activation due to the absence of a crystal structure. In the present study, homology modeling was applied to predicting OA1 structure following thorough sequence analysis and secondary structure predictions. The predicted model had the signature residues and motifs expected of GPCRs, and was used for carrying out molecular docking studies with an endogenous ligand, L-DOPA and an antagonist, dopamine; the results agreed quite well with the available experimental data. Finally, three sets of explicit molecular dynamics simulations were carried out in lipid bilayer, the results of which not only confirmed the stability of the predicted model, but also helped witness some differences in structural features such as rotamer toggle switch, helical tilts and hydrogen bonding pattern that helped distinguish between the agonist- and antagonist-bound receptor forms. In place of the typical "D/ERY"-motif-mediated "ionic lock", a hydrogen bond mediated by the "DAY" motif was observed that could be used to distinguish the agonist and antagonist bound forms of OA1. In the absence of a crystal structure, this study helped to shed some light on the structural features of OA1, and its behavior in the presence of an agonist and an antagonist, which might be helpful in the future drug discovery process for ocular albinism.

Modeling of ligand binding to G protein coupled receptors: cannabinoid CB1, CB2 and adrenergic β2AR

Cannabinoid and adrenergic receptors belong to the class A (similar to rhodopsin) G protein coupled receptors. Docking of agonists and antagonists to CB(1) and CB(2) cannabinoid receptors revealed the importance of a centrally located rotamer toggle switch and its possible participation in the mechanism of agonist/antagonist recognition. The switch is composed of two residues, F3.36 and W6.48, located on opposite transmembrane helices TM3 and TM6 in the central part of the membranous domain of cannabinoid receptors. The CB(1) and CB(2) receptor models were constructed based on the adenosine A(2A) receptor template. The two best scored conformations of each receptor were used for the docking procedure. In all poses (ligand-receptor conformations) characterized by the lowest ligand-receptor intermolecular energy and free energy of binding the ligand type matched the state of the rotamer toggle switch: antagonists maintained an inactive state of the switch, whereas agonists changed it. In case of agonists of β(2)AR, the (R,R) and (S,S) stereoisomers of fenoterol, the molecular dynamics simulations provided evidence of different binding modes while preserving the same average position of ligands in the binding site. The (S,S) isomer was much more labile in the binding site and only one stable hydrogen bond was created. Such dynamical binding modes may also be valid for ligands of cannabinoid receptors because of the hydrophobic nature of their ligand-receptor interactions. However, only very long molecular dynamics simulations could verify the validity of such binding modes and how they affect the process of activation.

Induced Effects of Sodium Ions on Dopaminergic G-Protein Coupled Receptors

G-protein coupled receptors, the largest family of proteins in the human genome, are involved in many complex signal transduction pathways, typically activated by orthosteric ligand binding and subject to allosteric modulation. Dopaminergic receptors, belonging to the class A family of G-protein coupled receptors, are known to be modulated by sodium ions from an allosteric binding site, although the details of sodium effects on the receptor have not yet been described. In an effort to understand these effects, we performed microsecond scale all-atom molecular dynamics simulations on the dopaminergic D2 receptor, finding that sodium ions enter the receptor from the extracellular side and bind at a deep allosteric site (Asp2.50). Remarkably, the presence of a sodium ion at this allosteric site induces a conformational change of the rotamer toggle switch Trp6.48 which locks in a conformation identical to the one found in the partially inactive state of the crystallized human β2 adrenergic receptor. This study provides detailed quantitative information about binding of sodium ions in the D2 receptor and reports a possibly important sodium-induced conformational change for modulation of D2 receptor function.

Comparison of the pharmacological properties of human and rat histamine H 3-receptors

Ligand pharmacology of histamine H 3-receptors is species-dependent. In previous studies, two amino acids in transmembrane domain 3 (TM III) were shown to play a significant role. In this study, we characterized human and rat histamine H 3-receptors (hH 3R and rH 3R, respectively), co-expressed with mammalian G proteins in Sf9 insect cell membranes. We compared a series of imidazole-containing H 3R ligands in radioligand binding and steady-state GTPase assays. H 3Rs similarly coupled to Gα i/o-proteins. Affinities and potencies of the agonists histamine, N α-methylhistamine and R-(α)-methylhistamine were in the same range. Imetit was only a partial agonist. The pharmacology of imetit and proxifan was similar at both species. However, impentamine was more potent and efficacious at rH 3R. The inverse agonists ciproxifan and thioperamide showed higher potency but lower efficacy at rH 3R. Clobenpropit was not species-selective. Strikingly, imoproxifan was almost full agonist at hH 3R, but an inverse agonist at rH 3R. Imoproxifan was docked into the binding pocket of inactive and active hH 3R- and rH 3R-models and molecular dynamic simulations were performed. Imoproxifan bound to hH 3R and rH 3R in E-configuration, which represents the trans-isomer of the oxime-moiety as determined in crystallization studies, and stabilized active hH 3R-, but inactive rH 3R-conformations. Large differences in electrostatic surfaces between TM III and TM V cause differential orientation of the oxime-moiety of imoproxifan, which then differently interacts with the rotamer toggle switch Trp 6.48 in TM VI. Collectively, the substantial species differences at H 3Rs are explained at a molecular level by the use of novel H 3R active-state models.

In silico analysis of the histaprodifen induced activation pathway of the guinea-pig histamine H1-receptor

The binding of (partial) agonists in the binding pocket of biogenic amine receptors induces a conformational change from the inactive to the active state of the receptors. There is only little knowledge about the binding pathways of ligands into binding pocket on molecular level. So far, it was not possible with molecular dynamic simulations to observe the ligand binding and receptor activation. Furthermore, there is nearly nothing known, in which state of ligand binding, the receptor gets activated. The aim of this study was to get more detailed insight into the process of ligand binding and receptor activation. With the recently developed LigPath algorithm, we scanned the potential energy surface of the binding process of dimeric histaprodifen, a partial agonist at the histamine H(1)-receptor, into the guinea pig histamine H(1)-receptor, taking also into account the receptor activation. The calculations exhibited large conformational changes of Trp(6.48) and Phe(6.55) during ligand binding and receptor activation. Additionally, conformational changes were also observed for Phe(6.52), Tyr(6.51) and Phe(6.44). Conformational changes of Trp(6.48) and Phe(6.52) are discussed in literature as rotamer toggle switch in context with receptor activation. Additionally, the calculations indicate that the binding of dimeric histaprodifen, accompanied by receptor activation is energetically preferred. In general, this study gives new, theoretical insights onto ligand binding and receptor activation on molecular level.

Structural Requirements for Inverse Agonism and Neutral Antagonism of Indole-, Benzimidazole-, and Thienopyrrole-Derived Histamine H4 Receptor Ligands

The human histamine H(4) receptor (hH(4)R), coexpressed with Galpha(i2) and Gbeta(1)gamma(2) in Sf9 insect cells, is highly constitutively active, and thioperamide [THIO; N-cyclohexyl-4-(imidazol-4-yl)-1-piperidinecarbothioamide] is one of the most efficacious hH(4)R inverse agonists. High constitutive hH(4)R activity may have pathophysiological implications in which case inverse agonists may behave differently than neutral antagonists. To learn more about the structural requirements for hH(4)R inverse agonism, we investigated 25 compounds (indole, benzimidazole, and thienopyrrole derivatives) structurally related to the standard antagonist JNJ-7777120 [1-[(5-chloro-1H-indol-2-yl)carbonyl]-4-methyl-piperazine]. We characterized the compounds in radioligand binding assays by using [(3)H]histamine ([(3)H]HA) and in steady-state GTPase assays in the presence (antagonist mode) and absence (inverse agonist mode) of the agonist HA, yielding the following results: 1) Twenty-two compounds were inverse agonists (efficacy: 15-62% of the THIO effect), and only three compounds (12%) showed neutral antagonism. Thus, inverse agonism is far more common than neutral antagonism. 2) The inverse agonistic efficacy of the R5-monosubstituted indole-derived compounds increased with the volume of R5. R5 may interact with Trp(6.48) of the rotamer toggle switch and stabilize the inactive receptor conformation. 3) A subset of compounds showed large differences between the K(i) value from [(3)H]HA competition binding and the EC(50) value from steady-state GTPase assays, whereas the K(b) values were closer to the K(i) values. Thus, the two-state model should be extended to a model comprising a constitutively active hH(4)R state, which can be discriminated by inverse agonists from a structurally distinct HA-stabilized active state.