Identification of residues required for ligand binding to the beta-adrenergic receptor

Molecular determinants of ligand efficacy and potency in GPCR signaling

Heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) bind to extracellular ligands and drugs and modulate intracellular responses through conformational changes. Despite their importance as drug targets, the molecular origins of pharmacological properties such as efficacy (maximum signaling response) and potency (the ligand concentration at half-maximal response) remain poorly understood for any ligand-receptor-signaling system. We used the prototypical adrenaline-β2 adrenergic receptor-G protein system to reveal how specific receptor residues decode and translate the information encoded in a ligand to mediate a signaling response. We present a data science framework to integrate pharmacological and structural data to uncover structural changes and allosteric networks relevant for ligand pharmacology. These methods can be tailored to study any ligand-receptor-signaling system, and the principles open possibilities for designing orthosteric and allosteric compounds with defined signaling properties.

Predicting the duration of action of β2-adrenergic receptor agonists: Ligand and structure-based approaches

Agonists of the β2 adrenergic receptor (ADRB2) are an important class of medications used for the treatment of respiratory diseases. They can be classified as short acting (SABA) or long acting (LABA), with each class playing a different role in patient management. In this work we explored both ligand-based and structure-based high-throughput approaches to classify β2-agonists based on their duration of action. A completely in-silico prediction pipeline using an AlphaFold generated structure was used for structure-based modelling. Our analysis identified the ligands' 3D structure and lipophilicity as the most relevant features for the prediction of the duration of action. Interaction-based methods were also able to select ligands with the desired duration of action, incorporating the bias directly in the structure-based drug discovery pipeline without the need for further processing.

Structural basis of agonist specificity of α1A-adrenergic receptor

α1-adrenergic receptors (α1-ARs) play critical roles in the cardiovascular and nervous systems where they regulate blood pressure, cognition, and metabolism. However, the lack of specific agonists for all α1 subtypes has limited our understanding of the physiological roles of different α1-AR subtypes, and led to the stagnancy in agonist-based drug development for these receptors. Here we report cryo-EM structures of α1A-AR in complex with heterotrimeric G-proteins and either the endogenous common agonist epinephrine or the α1A-AR-specific synthetic agonist A61603. These structures provide molecular insights into the mechanisms underlying the discrimination between α1A-AR and α1B-AR by A61603. Guided by the structures and corresponding molecular dynamics simulations, we engineer α1A-AR mutants that are not responsive to A61603, and α1B-AR mutants that can be potently activated by A61603. Together, these findings advance our understanding of the agonist specificity for α1-ARs at the molecular level, opening the possibility of rational design of subtype-specific agonists.

Naturally occurring mutations in G protein-coupled receptors associated with obesity and type 2 diabetes mellitus

G protein-coupled receptors (GPCRs) are the largest family of membrane receptors involved in the regulation of almost all known physiological processes. Dysfunctions of GPCR-mediated signaling have been shown to cause various diseases. The prevalence of obesity and type 2 diabetes mellitus (T2DM), two strongly associated disorders, is increasing worldwide, with tremendous economical and health burden. New safer and more efficacious drugs are required for successful weight reduction and T2DM treatment. Multiple GPCRs are involved in the regulation of energy and glucose homeostasis. Mutations in these GPCRs contribute to the development and progression of obesity and T2DM. Therefore, these receptors can be therapeutic targets for obesity and T2DM. Indeed some of these receptors, such as melanocortin-4 receptor and glucagon-like peptide 1 receptor, have provided important new drugs for treating obesity and T2DM. This review will focus on the naturally occurring mutations of several GPCRs associated with obesity and T2DM, especially incorporating recent large genomic data and insights from structure-function studies, providing leads for future investigations.

Distinctive communication networks in inactive states of β2 -adrenergic receptor: Mutual information and entropy transfer analysis

Mutual information and entropy transfer analysis employed on two inactive states of human beta-2 adrenergic receptor (β₂ -AR) unraveled distinct communication pathways. Previously, a so-called "highly" inactive state of the receptor was observed during 1.5 microsecond long molecular dynamics simulation where the largest intracellular loop (ICL3) was swiftly packed onto the G-protein binding cavity, becoming entirely inaccessible. Mutual information quantifying the degree of correspondence between backbone-C_α fluctuations was mostly shared between intra- and extra-cellular loop regions in the original inactive state, but shifted to entirely different regions in this latest inactive state. Interestingly, the largest amount of mutual information was always shared among the mobile regions. Irrespective of the conformational state, polar residues always contributed more to mutual information than hydrophobic residues, and also the number of polar-polar residue pairs shared the highest degree of mutual information compared to those incorporating hydrophobic residues. Entropy transfer, quantifying the correspondence between backbone-C_α fluctuations at different timesteps, revealed a distinctive pathway directed from the extracellular site toward intracellular portions in this recently exposed inactive state for which the direction of information flow was the reverse of that observed in the original inactive state where the mobile ICL3 and its intracellular surroundings drove the future fluctuations of extracellular regions.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

Exploring GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity. It is essential for an efficient structure-based drug design (SBDD) process. FMO enables ab initio approaches to be applied to systems that conventional quantum-mechanical (QM) methods would find challenging. The key advantage of the Fragment Molecular Orbital Method (FMO) is that it can reveal atomistic details about the individual contributions and chemical nature of each residue and water molecule toward ligand binding which would otherwise be difficult to detect without using QM methods. In this chapter, we demonstrate the typical use of FMO to analyze 19 crystal structures of β1 and β2 adrenergic receptors with their corresponding agonists and antagonists.

Astrocytic β2-adrenergic receptors mediate hippocampal long-term memory consolidation

Emotionally relevant experiences form strong and long-lasting memories by critically engaging the stress hormone/neurotransmitter noradrenaline, which mediates and modulates the consolidation of these memories. Noradrenaline acts through adrenergic receptors (ARs), of which β2-adrenergic receptors (βARs) are of particular importance. The differential anatomical and cellular distribution of βAR subtypes in the brain suggests that they play distinct roles in memory processing, although much about their specific contributions and mechanisms of action remains to be understood. Here we show that astrocytic rather than neuronal β2ARs in the hippocampus play a key role in the consolidation of a fear-based contextual memory. These hippocampal β2ARs, but not β1ARs, are coupled to the training-dependent release of lactate from astrocytes, which is necessary for long-term memory formation and for underlying molecular changes. This key metabolic role of astrocytic β2ARs may represent a novel target mechanism for stress-related psychopathologies and neurodegeneration.

Buried ionizable networks are an ancient hallmark of G protein-coupled receptor activation

Seven-transmembrane receptors (7TMRs) have evolved in prokaryotes and eukaryotes over hundreds of millions of years. Comparative structural analysis suggests that these receptors may share a remote evolutionary origin, despite their lack of sequence similarity. Here we used structure-based computations to compare 221 7TMRs from all domains of life. Unexpectedly, we discovered that these receptors contain spatially conserved networks of buried ionizable groups. In microbial 7TMRs these networks are used to pump ions across the cell membrane in response to light. In animal 7TMRs, which include light- and ligand-activated G protein-coupled receptors (GPCRs), homologous networks were found to be characteristic of activated receptor conformations. These networks are likely relevant to receptor function because they connect the ligand-binding pocket of the receptor to the nucleotide-binding pocket of the G protein. We propose that agonist and G protein binding facilitate the formation of these electrostatic networks and promote important structural rearrangements such as the displacement of transmembrane helix-6. We anticipate that robust classification of activated GPCR structures will aid the identification of ligands that target activated GPCR structural states.

A novel role of microglial NADPH oxidase in mediating extra-synaptic function of norepinephrine in regulating brain immune homeostasis

Although the peripheral anti-inflammatory effect of norepinephrine (NE) is well documented, the mechanism by which this neurotransmitter functions as an anti-inflammatory/neuroprotective agent in the central nervous system (CNS) is unclear. This article aimed to determine the anti-inflammatory/neuroprotective effects and underlying mechanisms of NE in inflammation-based dopaminergic neurotoxicity models. In mice, NE-depleting toxin N-(2-chloroethyl)-N-ethyl-2-bromobenzylamine (DSP-4) was injected at 6 months of lipopolysaccharide (LPS)-induced neuroinflammation. It was found that NE depletion enhanced LPS-induced dopaminergic neuron loss in the substantia nigra. This piece of in vivo data prompted us to conduct a series of studies in an effort to elucidate the mechanism as to how NE affects dopamine neuron survival by using primary midbrain neuron/glia cultures. Results showed that submicromolar concentrations of NE dose-dependently protected dopaminergic neurons from LPS-induced neurotoxicity by inhibiting microglia activation and subsequent release of pro-inflammatory factors. However, NE-elicited neuroprotection was not totally abolished in cultures from β2-adrenergic receptor (β2-AR)-deficient mice, suggesting that novel pathways other than β2-AR are involved. To this end, It was found that submicromolar NE dose-dependently inhibited NADPH oxidase (NOX2)-generated superoxide, which contributes to the anti-inflammatory and neuroprotective effects of NE. This novel mechanism was indeed adrenergic receptors independent since both (+) and (-) optic isomers of NE displayed the same potency. We further demonstrated that NE inhibited LPS-induced NOX2 activation by blocking the translocation of its cytosolic subunit to plasma membranes. In summary, we revealed a potential physiological role of NE in maintaining brain immune homeostasis and protecting neurons via a novel mechanism.

Mapping functional group free energy patterns at protein occluded sites: nuclear receptors and G-protein coupled receptors

Occluded ligand-binding pockets (LBP) such as those found in nuclear receptors (NR) and G-protein coupled receptors (GPCR) represent a significant opportunity and challenge for computer-aided drug design. To determine free energies maps of functional groups of these LBPs, a Grand-Canonical Monte Carlo/Molecular Dynamics (GCMC/MD) strategy is combined with the Site Identification by Ligand Competitive Saturation (SILCS) methodology. SILCS-GCMC/MD is shown to map functional group affinity patterns that recapitulate locations of functional groups across diverse classes of ligands in the LBPs of the androgen (AR) and peroxisome proliferator-activated-γ (PPARγ) NRs and the metabotropic glutamate (mGluR) and β₂-adreneric (β₂AR) GPCRs. Inclusion of protein flexibility identifies regions of the binding pockets not accessible in crystal conformations and allows for better quantitative estimates of relative ligand binding affinities in all the proteins tested. Differences in functional group requirements of the active and inactive states of the β₂AR LBP were used in virtual screening to identify high efficacy agonists targeting β₂AR in Airway Smooth Muscle (ASM) cells. Seven of the 15 selected ligands were found to effect ASM relaxation representing a 46% hit rate. Hence, the method will be of use for the rational design of ligands in the context of chemical biology and the development of therapeutic agents.

Human α1-adrenoceptor subtype selectivity of substituted homobivalent 4-aminoquinolines

A series of ring-substituted ethyl- and heptyl-linked 4-aminoquinoline dimers were synthesized and evaluated for their affinities at the 3 human α(1)-adrenoceptor (α(1)-AR) subtypes and the human serotonin 5-HT(1A)-receptor (5-HT(1A)-R). We find that the structure-specificity profiles are different for the two series at the α(1)-AR subtypes, which suggests that homobivalent 4-aminoquinolines can be developed with α(1)-AR subtype selectivity. The 8-methyl (8-Me) ethyl-linked analogue has the highest affinity for the α(1A)-AR, 7 nM, and the greatest capacity for discriminating between α(1A)-AR and α(1B)-AR (6-fold), α(1D)-AR (68-fold), and the 5-HT(1A)-R (168-fold). α(1B)-AR selectivity was observed with the 6-methyl (6-Me) derivative of the ethyl- and heptyl-linked 4-aminoquinoline dimers and the 7-methoxy (7-OMe) derivative of the heptyl-linked analogue. These substitutions result in 4- to 80-fold selectivity for α(1B)-AR over α(1A)-AR, α(1D)-AR, and 5-HT(1A)-R. In contrast, 4-aminoquinoline dimers with selectivity for α(1D)-AR are more elusive, since none studied to date has greater affinity for the α(1D)-AR over the other two α(1)-ARs. The selectivity of the 8-Me ethyl-linked 4-aminoquinoline dimer for the α(1A)-AR, and 6-Me ethyl-linked, and the 6-Me and 7-OMe heptyl-linked 4-aminoquinoline dimers for the α(1B)-AR, makes them promising leads for drug development of α(1A)-AR or α(1B)-AR subtype selective ligands with reduced 5-HT(1A)-R affinity.

Salmeterol's extreme β2 selectivity is due to residues in both extracellular loops and transmembrane domains

Salmeterol is a long-acting β2-agonist, widely used as an inhaled treatment of asthma and chronic obstructive pulmonary disease. It has very high β2-affinity (log KD -8.95) and is very selective for the β2-adrenoceptor (1000-fold selectivity over the β1-adrenoceptor). This study used a mutagenesis approach to determine the exact amino acids in the human β2-adrenoceptor responsible for this very high selectivity. Wild-type β2- and β1-adrenoceptors, chimeric β2/β1-adrenoceptors, and receptors with single-point mutations were transfected into Chinese hamster ovary-K1 cells, and affinity and function were studied using [(3)H]CGP 12177 [(-)-4-(3-tert-butylamino-2-hydroxypropoxy)-benzimidazol-2-one] whole-cell binding and [(3)H]cAMP accumulation. Extracellular loop 3 (and specifically amino acid K305) had the largest single effect by reducing salmeterol's affinity for the β2-adrenoceptor by 31-fold. H296 in transmembrane 6 also had a major effect (18-fold reduction in salmeterol affinity). Combining these, in the double mutant β2-H296K-K305D, reduced salmeterol's affinity by 275-fold, to within 4-fold of that of the β1-adrenoceptor, without affecting the affinity or selectivity of other β2-agonists (salbutamol, formoterol, fenoterol, clenbuterol, or adrenaline). Another important amino acid was Y308 in transmembrane 7, although this also affected the affinity and selectivity of other agonists. F194 in extracellular loop 2 and R304 in extracellular loop 3 also had minor effects. None of these mutations (including the double mutant β2-H296K-K305D) affected the efficacy or duration of action of salmeterol. This suggests that the high affinity and selectivity of salmeterol are due to specific amino acids within the receptor itself, but that the duration of action is at least in part due to other factors, for example lipophilicity.

Exploring the determinants of trace amine-associated receptor 1's functional selectivity for the stereoisomers of amphetamine and methamphetamine

Amphetamines are widely abused drugs that interfere with dopamine transport and storage. Recently, however, another mechanism of action was identified: stereoselective activation of the GαS protein-coupled trace amine-associated receptor 1 (TAAR1). To identify structural determinants of this stereoselectivity, we functionally evaluated six mutant receptors in vitro and then used homology modeling and dynamic simulation to predict drug affinities. Converting Asp102 to Ala rendered mouse and rat TAAR1 (mTAAR1 and rTAAR1, respectively) insensitive to β-phenylethylamine, amphetamine (AMPH), and methamphetamine (METH). Mutating Met268 in rTAAR1 to Thr shifted the concentration-response profiles for AMPH and METH isomers rightward an order of magnitude, whereas replacing Thr268 with Met in mTAAR1 resulted in profiles leftward shifted 10-30-fold. Replacing Asn287 with Tyr in rTAAR1 produced a mouselike receptor, while the reciprocal mTAAR1 mutant was rTAAR1-like. These results confirm TAAR1 is an AMPH/METH receptor in vitro and establish residues 102 (3.32) and 268 (6.55) as major contributors to AMPH/METH binding with residue 287 (7.39) determining species stereoselectivity.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Orthosteric binding of ρ-Da1a, a natural peptide of snake venom interacting selectively with the α1A-adrenoceptor

ρ-Da1a is a three-finger fold toxin from green mamba venom that is highly selective for the α_1A-adrenoceptor. This toxin has atypical pharmacological properties, including incomplete inhibition of ³H-prazosin or ¹²⁵I-HEAT binding and insurmountable antagonist action. We aimed to clarify its mode of action at the α_1A-adrenoceptor. The affinity (pKi 9.26) and selectivity of ρ-Da1a for the α_1A-adrenoceptor were confirmed by comparing binding to human adrenoceptors expressed in eukaryotic cells. Equilibrium and kinetic binding experiments were used to demonstrate that ρ-Da1a, prazosin and HEAT compete at the α_1A-adrenoceptor. ρ-Da1a did not affect the dissociation kinetics of ³H-prazosin or ¹²⁵I-HEAT, and the IC₅₀ of ρ-Da1a, determined by competition experiments, increased linearly with the concentration of radioligands used, while the residual binding by ρ-Da1a remained stable. The effect of ρ-Da1a on agonist-stimulated Ca²⁺ release was insurmountable in the presence of phenethylamine- or imidazoline-type agonists. Ten mutations in the orthosteric binding pocket of the α_1A-adrenoceptor were evaluated for alterations in ρ-Da1a affinity. The D106^3.32A and the S188^5.42A/S192^5.46A receptor mutations reduced toxin affinity moderately (6 and 7.6 times, respectively), while the F86^2.64A, F288^6.51A and F312^7.39A mutations diminished it dramatically by 18- to 93-fold. In addition, residue F86^2.64 was identified as a key interaction point for ¹²⁵I-HEAT, as the variant F86^2.64A induced a 23-fold reduction in HEAT affinity. Unlike the M1 muscarinic acetylcholine receptor toxin MT7, ρ-Da1a interacts with the human α_1A-adrenoceptor orthosteric pocket and shares receptor interaction points with antagonist (F86^2.64, F288^6.51 and F312^7.39) and agonist (F288^6.51 and F312^7.39) ligands. Its selectivity for the α_1A-adrenoceptor may result, at least partly, from its interaction with the residue F86^2.64, which appears to be important also for HEAT binding.

The N-terminal region of the dopamine D2 receptor, a rhodopsin-like GPCR, regulates correct integration into the plasma membrane and endocytic routes

BACKGROUND AND PURPOSE

Functional roles of the N-terminal region of rhodopsin-like GPCR family remain unclear. Using dopamine D(2) and D(3) receptors as a model system, we probed the roles of the N-terminal region in the signalling, intracellular trafficking of receptor proteins, and explored the critical factors that determine the functionality of the N-terminal region.

EXPERIMENTAL APPROACH

The N-terminal region of the D(2) receptor was gradually shortened or switched with that of the D(3) receptor or a non-specific sequence (FLAG), or potential N-terminal glycosylation sites were mutated. Effects of these manipulations on surface expression, internalization, post-endocytic behaviours and signalling were determined.

KEY RESULTS

Shortening the N-terminal region of the D(2) receptor enhanced receptor internalization and impaired surface expression and signalling; ligand binding, desensitization and down-regulation were not affected but their association with a particular microdomain, caveolae, was disrupted. Replacement of critical residues within the N-terminal region with the FLAG epitope failed to restore surface expression but partially restored the altered internalization and signalling. When the N-terminal regions were switched between D(2) and D(3) receptors, cell surface expression pattern of each receptor was switched. Mutations of potential N-terminal glycosylation sites inhibited surface expression but enhanced internalization of D(2) receptors.

CONCLUSIONS AND IMPLICATIONS

Shortening of N-terminus or mutation of glycosylation sites located within the N-terminus enhanced receptor internalization but impaired the surface expression of D(2) receptors. The N-terminal region of the D(2) receptor, in a sequence-specific manner, controls the receptor's conformation and integration into the plasma membrane, which determine its subcellular localization, intracellular trafficking and signalling properties.

Reversed binding of a small molecule ligand in homologous chemokine receptors - differential role of extracellular loop 2

BACKGROUND AND PURPOSE

The majority of small molecule compounds targeting chemokine receptors share a similar pharmacophore with a centrally located aliphatic positive charge and flanking aromatic moieties. Here we describe a novel piperidine-based compound with structural similarity to previously described CCR8-specific agonists, but containing a unique phenyl-tetrazol moiety which, in addition to activity at CCR8 was also active at CCR1.

EXPERIMENTAL APPROACH

Single point mutations were introduced in CCR1 and CCR8, and their effect on small molecule ligand-induced receptor activation was examined through inositol trisphosphate (IP(3) ) accumulation. The molecular interaction profile of the agonist was verified by molecular modeling.

KEY RESULTS

The chemokine receptor conserved glutamic acid in TM-VII served as a common anchor for the positively charged amine in the piperidine ring. However, whereas the phenyl-tetrazol group interacted with TyrIV:24 (Tyr(172) ) and TyrIII:09 (Tyr(114) ) in the major binding pocket (delimited by TM-III to VII) of CCR8, it also interacted with TrpII:20 (Trp(90) ) and LysII:24 (Lys(94) ) in the minor counterpart (delimited TM-I to III, plus TM-VII) in CCR1. A straightening of TM-II by Ala-substitution of ProII:18 confirmed its unique role in CCR1. The extracellular loop 2 (ECL-2) contributed directly to the small molecule binding site in CCR1, whereas it contributed to efficacy, but not potency in CCR8.

CONCLUSION AND IMPLICATIONS

Despite high ligand potency and efficacy and receptor similarity, this dual-active and bitopic compound binds oppositely in CCR1 and CCR8 with different roles of ECL-2, thereby expanding and diversifying the influence of extracellular receptor regions in drug action.

A prospective cross-screening study on G-protein-coupled receptors: lessons learned in virtual compound library design

We present the systematic prospective evaluation of a protein-based and a ligand-based virtual screening platform against a set of three G-protein-coupled receptors (GPCRs): the β-2 adrenoreceptor (ADRB2), the adenosine A(2A) receptor (AA2AR), and the sphingosine 1-phosphate receptor (S1PR1). Novel bioactive compounds were identified using a consensus scoring procedure combining ligand-based (frequent substructure ranking) and structure-based (Snooker) tools, and all 900 selected compounds were screened against all three receptors. A striking number of ligands showed affinity/activity for GPCRs other than the intended target, which could be partly attributed to the fuzziness and overlap of protein-based pharmacophore models. Surprisingly, the phosphodiesterase 5 (PDE5) inhibitor sildenafil was found to possess submicromolar affinity for AA2AR. Overall, this is one of the first published prospective chemogenomics studies that demonstrate the identification of novel cross-pharmacology between unrelated protein targets. The lessons learned from this study can be used to guide future virtual ligand design efforts.

Involvement of the first transmembrane segment of human α(2) -adrenoceptors in the subtype-selective binding of chlorpromazine, spiperone and spiroxatrine

BACKGROUND AND PURPOSE

Some large antagonist ligands (ARC239, chlorpromazine, prazosin, spiperone, spiroxatrine) bind to the human α(2A) -adrenoceptor with 10- to 100-fold lower affinity than to the α(2B)- and α(2C)-adrenoceptor subtypes. Previous mutagenesis studies have not explained this subtype selectivity.

EXPERIMENTAL APPROACH

The possible involvement of the extracellular amino terminus and transmembrane domain 1 (TM1) in subtype selectivity was elucidated with eight chimaeric receptors: six where TM1 and the N-terminus were exchanged between the α(2)-adrenoceptor subtypes and two where only TM1 was exchanged. Receptors were expressed in CHO cells and tested for ligand binding with nine chemically diverse antagonist ligands. For purposes of interpretation, molecular models of the three human α(2)-adrenoceptors were constructed based on the β(2)-adrenoceptor crystal structure.

KEY RESULTS

The affinities of three antagonists (spiperone, spiroxatrine and chlorpromazine) were significantly improved by TM1 substitutions of the α(2A)-adrenoceptor, but reciprocal effects were not seen for chimaeric receptors based on α(2B)- and α(2C)-adrenoceptors. Molecular docking of these ligands suggested that binding occurs in the orthosteric ligand binding pocket.

CONCLUSIONS AND IMPLICATIONS

TM1 is involved in determining the low affinity of some antagonist ligands at the human α(2A)-adrenoceptor. The exact mechanism is not known, but the position of TM1 at a large distance from the binding pocket indicates that TM1 does not participate in specific side-chain interactions with amino acids within the binding pocket of the receptor or with ligands bound therein. Instead, molecular models suggest that TM1 has indirect conformational effects related to the charge distribution or overall shape of the binding pocket.

Ligand modulation of the Epstein-Barr virus-induced seven-transmembrane receptor EBI2: identification of a potent and efficacious inverse agonist

The Epstein-Barr virus-induced receptor 2 (EBI2) is a constitutively active seven-transmembrane receptor, which was recently shown to orchestrate the positioning of B cells in the follicle. To date, no ligands, endogenously or synthetic, have been identified that modulate EBI2 activity. Here we describe an inverse agonist, GSK682753A, which selectively inhibited the constitutive activity of EBI2 with high potency and efficacy. In cAMP-response element-binding protein-based reporter and guanosine 5'-3-O-(thio)triphosphate (GTPγS) binding assays, the potency of this compound was 2.6-53.6 nm, and its inhibitory efficacy was 75%. In addition, we show that EBI2 constitutively activated extracellular signal-regulated kinase (ERK) in a pertussis toxin-insensitive manner. Intriguingly, GSK682753A inhibited ERK phosphorylation, GTPγS binding, and cAMP-response element-binding protein activation with similar potency. Overexpression of EBI2 profoundly potentiated antibody-stimulated ex vivo proliferation of murine B cells compared with WT cells, whereas this was equivalently reduced for EBI2-deficient B cells. Inhibition of EBI2 constitutive activity suppressed the proliferation in all cases. Importantly, the suppression was of much higher potency (32-fold) in WT or EBI2-overexpressing B cells compared with EBI2-deficient counterparts. Finally, we screened GSK682753A against an EBI2 mutant library to determine putative molecular binding determinants in EBI2. We identified Phe(111) at position III:08/3.32 as being crucial for GSK682753A inverse agonism because Ala substitution resulted in a >500-fold decrease in IC(50). In conclusion, we present the first ligand targeting EBI2. In turn, this molecule provides a useful tool for further characterization of EBI2 as well as serving as a potent lead compound.

Analysis of multiple compound-protein interactions reveals novel bioactive molecules

The authors use machine learning of compound-protein interactions to explore drug polypharmacology and to efficiently identify bioactive ligands, including novel scaffold-hopping compounds for two pharmaceutically important protein families: G-protein coupled receptors and protein kinases.

We have demonstrated that machine learning of multiple compound–protein interactions is useful for efficient ligand screening and for assessing drug polypharmacology.

This approach successfully identified novel scaffold-hopping compounds for two pharmaceutically important protein families: G-protein-coupled receptors and protein kinases.

These bioactive compounds were not detected by existing computational ligand-screening methods in comparative studies.

The results of this study indicate that data derived from chemical genomics can be highly useful for exploring chemical space, and this systems biology perspective could accelerate drug discovery processes.

The discovery of novel bioactive molecules advances our systems-level understanding of biological processes and is crucial for innovation in drug development. Perturbations of biological systems by chemical probes provide broader applications not only for analysis of complex systems but also for intentional manipulations of these systems. Nevertheless, the lack of well-characterized chemical modulators has limited their use. Recently, chemical genomics has emerged as a promising area of research applicable to the exploration of novel bioactive molecules, and researchers are currently striving toward the identification of all possible ligands for all target protein families (Wang et al, 2009). Chemical genomics studies have shown that patterns of compound–protein interactions (CPIs) are too diverse to be understood as simple one-to-one events. There is an urgent need to develop appropriate data mining methods for characterizing and visualizing the full complexity of interactions between chemical space and biological systems. However, no existing screening approach has so far succeeded in identifying novel bioactive compounds using multiple interactions among compounds and target proteins.

High-throughput screening (HTS) and computational screening have greatly aided in the identification of early lead compounds for drug discovery. However, the large number of assays required for HTS to identify drugs that target multiple proteins render this process very costly and time-consuming. Therefore, interest in using in silico strategies for screening has increased. The most common computational approaches, ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS; Oprea and Matter, 2004; Muegge and Oloff, 2006; McInnes, 2007; Figure 1A), have been used for practical drug development. LBVS aims to identify molecules that are very similar to known active molecules and generally has difficulty identifying compounds with novel structural scaffolds that differ from reference molecules. The other popular strategy, SBVS, is constrained by the number of three-dimensional crystallographic structures available. To circumvent these limitations, we have shown that a new computational screening strategy, chemical genomics-based virtual screening (CGBVS), has the potential to identify novel, scaffold-hopping compounds and assess their polypharmacology by using a machine-learning method to recognize conserved molecular patterns in comprehensive CPI data sets.

The CGBVS strategy used in this study was made up of five steps: CPI data collection, descriptor calculation, representation of interaction vectors, predictive model construction using training data sets, and predictions from test data (Figure 1A). Importantly, step 1, the construction of a data set of chemical structures and protein sequences for known CPIs, did not require the three-dimensional protein structures needed for SBVS. In step 2, compound structures and protein sequences were converted into numerical descriptors. These descriptors were used to construct chemical or biological spaces in which decreasing distance between vectors corresponded to increasing similarity of compound structures or protein sequences. In step 3, we represented multiple CPI patterns by concatenating these chemical and protein descriptors. Using these interaction vectors, we could quantify the similarity of molecular interactions for compound–protein pairs, despite the fact that the ligand and protein similarity maps differed substantially. In step 4, concatenated vectors for CPI pairs (positive samples) and non-interacting pairs (negative samples) were input into an established machine-learning method. In the final step, the classifier constructed using training sets was applied to test data.

To evaluate the predictive value of CGBVS, we first compared its performance with that of LBVS by fivefold cross-validation. CGBVS performed with considerably higher accuracy (91.9%) than did LBVS (84.4%; Figure 1B). We next compared CGBVS and SBVS in a retrospective virtual screening based on the human β2-adrenergic receptor (ADRB2). Figure 1C shows that CGBVS provided higher hit rates than did SBVS. These results suggest that CGBVS is more successful than conventional approaches for prediction of CPIs.

We then evaluated the ability of the CGBVS method to predict the polypharmacology of ADRB2 by attempting to identify novel ADRB2 ligands from a group of G-protein-coupled receptor (GPCR) ligands. We ranked the prediction scores for the interactions of 826 reported GPCR ligands with ADRB2 and then analyzed the 50 highest-ranked compounds in greater detail. Of 21 commercially available compounds, 11 showed ADRB2-binding activity and were not previously reported to be ADRB2 ligands. These compounds included ligands not only for aminergic receptors but also for neuropeptide Y-type 1 receptors (NPY1R), which have low protein homology to ADRB2. Most ligands we identified were not detected by LBVS and SBVS, which suggests that only CGBVS could identify this unexpected cross-reaction for a ligand developed as a target to a peptidergic receptor.

The true value of CGBVS in drug discovery must be tested by assessing whether this method can identify scaffold-hopping lead compounds from a set of compounds that is structurally more diverse. To assess this ability, we analyzed 11 500 commercially available compounds to predict compounds likely to bind to two GPCRs and two protein kinases. Functional assays revealed that nine ADRB2 ligands, three NPY1R ligands, five epidermal growth factor receptor (EGFR) inhibitors, and two cyclin-dependent kinase 2 (CDK2) inhibitors were concentrated in the top-ranked compounds (hit rate=30, 15, 25, and 10%, respectively). We also evaluated the extent of scaffold hopping achieved in the identification of these novel ligands. One ADRB2 ligand, two NPY1R ligands, and one CDK2 inhibitor exhibited scaffold hopping (Figure 4), indicating that CGBVS can use this characteristic to rationally predict novel lead compounds, a crucial and very difficult step in drug discovery. This feature of CGBVS is critically different from existing predictive methods, such as LBVS, which depend on similarities between test and reference ligands, and focus on a single protein or highly homologous proteins. In particular, CGBVS is useful for targets with undefined ligands because this method can use CPIs with target proteins that exhibit lower levels of homology.

In summary, we have demonstrated that data mining of multiple CPIs is of great practical value for exploration of chemical space. As a predictive model, CGBVS could provide an important step in the discovery of such multi-target drugs by identifying the group of proteins targeted by a particular ligand, leading to innovation in pharmaceutical research.

The discovery of novel bioactive molecules advances our systems-level understanding of biological processes and is crucial for innovation in drug development. For this purpose, the emerging field of chemical genomics is currently focused on accumulating large assay data sets describing compound–protein interactions (CPIs). Although new target proteins for known drugs have recently been identified through mining of CPI databases, using these resources to identify novel ligands remains unexplored. Herein, we demonstrate that machine learning of multiple CPIs can not only assess drug polypharmacology but can also efficiently identify novel bioactive scaffold-hopping compounds. Through a machine-learning technique that uses multiple CPIs, we have successfully identified novel lead compounds for two pharmaceutically important protein families, G-protein-coupled receptors and protein kinases. These novel compounds were not identified by existing computational ligand-screening methods in comparative studies. The results of this study indicate that data derived from chemical genomics can be highly useful for exploring chemical space, and this systems biology perspective could accelerate drug discovery processes.

Remote control of neuronal signaling

A significant challenge for neuroscientists is to determine how both electrical and chemical signals affect the activity of cells and circuits and how the nervous system subsequently translates that activity into behavior. Remote, bidirectional manipulation of those signals with high spatiotemporal precision is an ideal approach to addressing that challenge. Neuroscientists have recently developed a diverse set of tools that permit such experimental manipulation with varying degrees of spatial, temporal, and directional control. These tools use light, peptides, and small molecules to primarily activate ion channels and G protein-coupled receptors (GPCRs) that in turn activate or inhibit neuronal firing. By monitoring the electrophysiological, biochemical, and behavioral effects of such activation/inhibition, researchers can better understand the links between brain activity and behavior. Here, we review the tools that are available for this type of experimentation. We describe the development of the tools and highlight exciting in vivo data. We focus primarily on designer GPCRs (receptors activated solely by synthetic ligands, designer receptors exclusively activated by designer drugs) and microbial opsins (e.g., channelrhodopsin-2, halorhodopsin, Volvox carteri channelrhodopsin) but also describe other novel techniques that use orthogonal receptors, caged ligands, allosteric modulators, and other approaches. These tools differ in the direction of their effect (activation/inhibition, hyperpolarization/depolarization), their onset and offset kinetics (milliseconds/minutes/hours), the degree of spatial resolution they afford, and their invasiveness. Although none of these tools is perfect, each has advantages and disadvantages, which we describe, and they are all still works in progress. We conclude with suggestions for improving upon the existing tools.

Structural insights into adrenergic receptor function and pharmacology

It has been over 50years since Sir James Black developed the first beta adrenergic receptor (βAR) blocker to treat heart disease. At that time, the concept of cell surface receptors was relatively new and not widely accepted, and most of the tools currently used to characterize plasma membrane receptors had not been developed. There has been remarkable progress in receptor biology since then, including the development of radioligand binding assays, the biochemical characterization of receptors as discrete membrane proteins, and the cloning of the first G-protein-coupled receptors (GPCRs), which led to the identification of other members of the large family of GPCRs. More recently, progress in GPCR structural biology has led to insights into the three-dimensional structures of βARs in both active and inactive states. Despite all of this progress, the process of developing a drug for a particular GPCR target has become more complex, time-consuming and expensive.

Internalization dissociates β2-adrenergic receptors

G protein-coupled receptors (GPCRs) self-associate as dimers or higher-order oligomers in living cells. The stability of associated GPCRs has not been extensively studied, but it is generally thought that these receptors move between the plasma membrane and intracellular compartments as intact dimers or oligomers. Here we show that β(2)-adrenergic receptors (β(2)ARs) that self-associate at the plasma membrane can dissociate during agonist-induced internalization. We use bioluminescence-resonance energy transfer (BRET) to monitor movement of β(2)ARs between subcellular compartments. BRET between β(2)ARs and plasma membrane markers decreases in response to agonist activation, while at the same time BRET between β(2)ARs and endosome markers increases. Energy transfer between β(2)ARs is decreased in a similar manner if either the donor- or acceptor-labeled receptor is mutated to impair agonist binding and internalization. These changes take place over the course of 30 minutes, persist after agonist is removed, and are sensitive to several inhibitors of arrestin- and clathrin-mediated endocytosis. The magnitude of the decrease in BRET between donor- and acceptor-labeled β(2)ARs suggests that at least half of the receptors that contribute to the BRET signal are physically segregated by internalization. These results are consistent with the possibility that β(2)ARs associate transiently with each other in the plasma membrane, or that β(2)AR dimers or oligomers are actively disrupted during internalization.

Structure of a nanobody-stabilized active state of the β(2) adrenoceptor

G protein coupled receptors (GPCRs) exhibit a spectrum of functional behaviors in response to natural and synthetic ligands. Recent crystal structures provide insights into inactive states of several GPCRs. Efforts to obtain an agonist-bound active-state GPCR structure have proven difficult due to the inherent instability of this state in the absence of a G protein. We generated a camelid antibody fragment (nanobody) to the human β₂ adrenergic receptor (β₂AR) that exhibits G protein-like behavior, and obtained an agonist-bound, active-state crystal structure of the receptor-nanobody complex. Comparison with the inactive β₂AR structure reveals subtle changes in the binding pocket; however, these small changes are associated with an 11Å outward movement of the cytoplasmic end of transmembrane segment 6, and rearrangements of transmembrane segments 5 and 7 that are remarkably similar to those observed in opsin, an active form of rhodopsin. This structure provides insights into the process of agonist binding and activation.

Hydrophobicity profiles in G protein-coupled receptor transmembrane helical domains

The lack of a crystallographically derived structure for all but three G (TP [guanosine triphosphate]-binding) protein-coupled receptor (GPCRs) proteins necessitates the use of computationally derived methods to determine their structures. Computational methodologies allow a mechanistic glimpse into GPCR-ligand interactions at a molecular level to better understand the initial steps leading to a protein's biologic functions, ie, protecting the ligands that activate, deactivate, or inhibit the protein, stabilizing protein structure in the membrane's lipid bilayer, and ensuring that the hydrophilic environment within the GPCR-binding pocket is maintained. Described here is a formalism that quantifies the amphiphilic nature of a helix, by determining the effective hydrophobicity (or hydrophilicity) at specific positions around it. This formalism will enable computational protein modelers to position helices so that the functional aspects of GPCRs are adequately represented in the model. Hydro-Eff®, an online tool, allows users to calculate effective helical hydrophobicities.

Ultrastructural characterization of noradrenergic axons and Beta-adrenergic receptors in the lateral nucleus of the amygdala

Norepinephrine (NE) is thought to play a key role in fear and anxiety, but its role in amygdala-dependent Pavlovian fear conditioning, a major model for understanding the neural basis of fear, is poorly understood. The lateral nucleus of the amygdala (LA) is a critical brain region for fear learning and regulating the effects of stress on memory. To understand better the cellular mechanisms of NE and its adrenergic receptors in the LA, we used antibodies directed against dopamine beta-hydroxylase (DβH), the synthetic enzyme for NE, or against two different isoforms of the beta-adrenergic receptors (βARs), one that predominately recognizes neurons (βAR 248) and the other astrocytes (βAR 404), to characterize the microenvironments of DβH and βAR. By electron microscopy, most DβH terminals did not make synapses, but when they did, they formed both asymmetric and symmetric synapses. By light microscopy, βARs were present in both neurons and astrocytes. Confocal microscopy revealed that both excitatory and inhibitory neurons express βAR248. By electron microscopy, βAR 248 was present in neuronal cell bodies, dendritic shafts and spines, and some axon terminals and astrocytes. When in dendrites and spines, βAR 248 was frequently concentrated along plasma membranes and at post-synaptic densities of asymmetric (excitatory) synapses. βAR 404 was expressed predominately in astrocytic cell bodies and processes. These astrocytic processes were frequently interposed between unlabeled terminals or ensheathed asymmetric synapses. Our findings provide a morphological basis for understanding ways in which NE may modulate transmission by acting via synaptic or non-synaptic mechanisms in the LA.

The importance of valine 114 in ligand binding in beta(2)-adrenergic receptor

G-protein coupled receptors (GPCRs) are transmembrane signaling molecules, with a majority of them performing important physiological roles. beta(2)-Adrenergic receptor (beta(2)-AR) is a well-studied GPCRs that mediates natural responses to the hormones adrenaline and noradrenaline. Analysis of the ligand-binding region of beta(2)-AR using the recently solved high-resolution crystal structures revealed a number of highly conserved amino acids that might be involved in ligand binding. However, detailed structure-function studies on some of these residues have not been performed, and their role in ligand binding remains to be elucidated. In this study, we have investigated the structural and functional role of a highly conserved residue valine 114, in hamster beta(2)-AR by site-directed mutagenesis. We replaced V114 in hamster beta(2)-AR with a number of amino acid residues carrying different functional groups. In addition to the complementary substitutions V114I and V114L, the V114C and V114E mutants also showed significant ligand binding and agonist dependent G-protein activation. However, the V114G, V114T, V114S, and V114W mutants failed to bind ligand in a specific manner. Molecular modeling studies were conducted to interpret these results in structural terms. We propose that the replacement of V114 influences not only the interaction of the ethanolamine side-chains but also the aryl-ring of the ligands tested. Results from this study show that the size and orientation of the hydrophobic residue at position V114 in beta(2)-AR affect binding of both agonists and antagonists, but it does not influence the receptor expression or folding.

Amino acid residues critical for endoplasmic reticulum export and trafficking of platelet-activating factor receptor

Several residues are conserved in the transmembrane domains (TMs) of G-protein coupled receptors. Here we demonstrate that a conserved proline, Pro(247), in TM6 of platelet-activating factor receptor (PAFR) is required for endoplasmic reticulum (ER) export and trafficking after agonist-induced internalization. Alanine-substituted mutants of the conserved residues of PAFRs, including P247A, were retained in the ER. Because a PAFR antagonist, Y-24180, acted as a pharmacological chaperone to rescue ER retention, this retention is due to misfolding of PAFR. Methylcarbamyl (mc)-PAF, a PAFR agonist, did not increase the cell surface expression of P247A, even though another ER-retained mutant, D63A, was effectively trafficked. Signaling and accumulation of the receptors in the early endosomes were observed in the mc-PAF-treated P247A-expressing cells, suggesting that P247A was trafficked to the cell surface by mc-PAF, and thereafter disappeared from the surface due to aberrant trafficking, e.g. enhanced internalization, deficiency in recycling, and/or accelerated degradation. The aberrant trafficking was confirmed with a sortase-A-mediated method for labeling cell surface proteins. These results demonstrate that the conserved proline in TM6 is crucial for intracellular trafficking of PAFR.

Mutation of important amino acid residue of Asp104Lys in human beta(1)-adrenergic receptor triggers functional and constitutive inactivation

Based on our previous molecular modeling and radioligand binding study, we have demonstrated that aspartic acid of 104 in transmembrane helix (TMH) II of beta(1)-adrenergic receptor (beta(1)-AR) is important for functional characteristics of these receptors. We have also showed that mutation of negatively charged aspartic acid to neutral charged alanine exhibited constitutive activity of beta(1)-AR. However, the mutation of negatively charged aspartic acid to positively charged lysine is still remained to be examined, which is very important to know for fully understanding the characteristics of beta(1)-AR. At the present study, we mutated aspartic acid to lysine (Asp104Lys) residue in human beta(1)-AR. This resultant mutant (Asp104Lys) markedly reduced the binding affinity of isoproterenol and (-)-epinephrine. On the other hand, antagonist binding with this mutant was similar to the wild type receptor. Isoproterenol at its saturation concentrations produced lower amount of intracellular cyclic adenosine-3',5' cyclic monophosphate (cAMP) in HEK-293 cells expressing Asp104Lys mutant receptor as compared to cells expressing wild type receptor. Moreover, cAMP accumulation of Asp104Lys mutant was unchanged in the presence or absence of isoproterenol. Therefore, it has been demonstrated that Asp104Lys mutation in the human beta(1)-AR differentially affects the binding of antagonist and exhibits a functional uncoupling of G-protein-coupled receptors. Thus, we may suggest that mutation of negatively charged aspartic acid to positively charged lysine as well as neutral charged alanine may help to understand the mechanism of the activation or inactivation of beta(1)-AR by its conformational changes and this finding would be helpful for clarifying the functional responses mediated by beta(1)-AR.

Structural basis for ligand binding and specificity in adrenergic receptors: implications for GPCR-targeted drug discovery

Crystal structures of engineered human beta 2-adrenergic receptors (ARs) in complex with an inverse agonist ligand, carazolol, provide three-dimensional snapshots of the disposition of seven transmembrane helices and the ligand-binding site of an important G protein-coupled receptor (GPCR). As expected, beta 2-AR shares substantial structural similarities with rhodopsin, the dim-light photoreceptor of the rod cell. However, although carazolol and the 11- cis-retinylidene moiety of rhodopsin are situated in the same general binding pocket, the second extracellular (E2) loop structures are quite distinct. E2 in rhodopsin shows beta-sheet structure and forms part of the chromophore-binding site. In the beta 2-AR, E2 is alpha-helical and seems to be distinct from the receptor's active site, allowing a potential entry pathway for diffusible ligands. The structures, together with extensive structure-activity relationship (SAR) data from earlier studies, provide insight about possible structural determinants of ligand specificity and how the binding of agonist ligands might alter receptor conformation. We review key features of the new beta 2-AR structures in the context of recent complementary work on the conformational dynamics of GPCRs. We also report 600 ns molecular dynamics simulations that quantified beta 2-AR receptor mobility in a membrane bilayer environment and show how the binding of an agonist ligand, adrenaline (epinephrine), causes conformational changes to the ligand-binding pocket and neighboring helices.

Selective structure-based virtual screening for full and partial agonists of the beta2 adrenergic receptor

The recently solved high-resolution X-ray structure of the beta2 adrenergic receptor has been challenged for its ability to discriminate inverse agonists/antagonists from partial/full agonists. Whereas the X-ray structure of the ground state receptor was unsuitable to distinguish true ligands with different functional effects, modifying this structure to reflect early conformational events in receptor activation led to a receptor model able to selectively retrieve full and partial agonists by structure-based virtual screening. The use of a topological scoring function based on molecular interaction fingerprints was shown to be mandatory to properly rank docking poses and achieve acceptable enrichments for partial and full agonists only.

Agonist binding, agonist affinity and agonist efficacy at G protein-coupled receptors

Measurements of affinity and efficacy are fundamental for work on agonists both in drug discovery and in basic studies on receptors. In this review I wish to consider methods for measuring affinity and efficacy at G protein coupled receptors (GPCRs). Agonist affinity may be estimated in terms of the dissociation constant for agonist binding to a receptor using ligand binding or functional assays. It has, however, been suggested that measurements of affinity are always contaminated by efficacy so that it is impossible to separate the two parameters. Here I show that for many GPCRs, if receptor/G protein coupling is suppressed, experimental measurements of agonist affinity using ligand binding (K(obs)) provide quite accurate measures of the agonist microscopic dissociation constant (KA). Also in pharmacological functional studies, good estimates of agonist dissociation constants are possible. Efficacy can be quantitated in several ways based on functional data (maximal effect of the agonist (E(max)), ratio of agonist dissociation constant to concentration of agonist giving half maximal effect in functional assay (K(obs)/EC50), a combined parameter E(max)K(obs)/EC50). Here I show that E(max)K(obs)/EC50 provides the best assessment of efficacy for a range of agonists across the full range of efficacy for full to partial agonists. Considerable evidence now suggests that ligand efficacy may be dependent on the pathway used to assess it. The efficacy of a ligand may, therefore, be multidimensional. It is still, however, necessary to have accurate measures of efficacy in different pathways.

Asp125 and Thr130 in transmembrane domain 3 are major sites of alpha1b-adrenergic receptor antagonist binding

Site-directed mutagenesis was used to investigate the molecular interactions involved in prazosin binding to the human alpha(1b)-adrenergic receptor (alpha(1b)-AR) receptor. Based on molecular modeling studies, Thr130 and Asp125 in transmembrane region III of the alpha(1b)-AR receptor were found to interact with prazosin. Thr130 and Asp125 were mutated to alanine (Ala) and expressed in HEK293 cells. The radioligand [(3)H]prazosin did not show any binding to Asp125Ala mutant of alpha(1b)-AR. Therefore, it was not possible to find any prazosin affinity to the mutant using the radioligand [(3)H]prazosin. The mutation also abolished phenylephrine-stimulated inositol phosphate (IP) formation of [(3)H]myo-inositol. On the other hand, the Thr130Ala mutant showed reduced binding affinity for [(3)H]prazosin (dissociation constant, K(d) 674.27 pM versus 90.27 pM for the wild-type receptor) and had reduced affinity for both tamsulosin and prazosin (11-fold and 9-fold, respectively). However, the Thr130Ala mutant receptor retained the ability to stimulate the formation of [(3)H]myo-inositol. The results provide direct evidence that Asp125 and Thr130 are responsible for the interactions between alpha(1b)-AR receptor and radioligand [(3)H]prazosin as well as tamsulosin.

GPCR engineering yields high-resolution structural insights into beta2-adrenergic receptor function

The beta2-adrenergic receptor (beta2AR) is a well-studied prototype for heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) that respond to diffusible hormones and neurotransmitters. To overcome the structural flexibility of the beta2AR and to facilitate its crystallization, we engineered a beta2AR fusion protein in which T4 lysozyme (T4L) replaces most of the third intracellular loop of the GPCR ("beta2AR-T4L") and showed that this protein retains near-native pharmacologic properties. Analysis of adrenergic receptor ligand-binding mutants within the context of the reported high-resolution structure of beta2AR-T4L provides insights into inverse-agonist binding and the structural changes required to accommodate catecholamine agonists. Amino acids known to regulate receptor function are linked through packing interactions and a network of hydrogen bonds, suggesting a conformational pathway from the ligand-binding pocket to regions that interact with G proteins.

Two model system of the alpha1A-adrenoceptor docked with selected ligands

In this study, we have developed a two model system to mimic the active and inactive states of a G-protein coupled receptor specifically the alpha1A adrenergic receptor. We have docked two agonists, epinephrine (phenylamine type) and oxymetazoline (imidazoline type), as well as two antagonists, prazosin and 5-methylurapidil, into two alpha1A receptor models, active and inactive. The best docking complexes for both agonists had hydrophilic interactions with D106, while neither antagonist did. Prazosin and oxymetazoline had hydrophobic interactions with F308 and F312. We predict from our study that the active state is stabilized by the interaction of F193 with I114, L197, V278, F281, and V282. The active state is further stabilized by the interaction of F312 with L75, V79, and L80. We also predict that the inactive state of the receptor is stabilized by the interaction of F312 with W102, F288, and M292.

The second extracellular loop of alpha2A-adrenoceptors contributes to the binding of yohimbine analogues

BACKGROUND AND PURPOSE

Rodent alpha(2A)-adrenoceptors bind the classical alpha(2)-antagonists yohimbine and rauwolscine with lower affinity than the human alpha(2A)-adrenoceptor. A serine-cysteine difference in the fifth transmembrane helix (TM; position 5.43) partially explains this, but all determinants of the interspecies binding selectivity are not known. Molecular models of alpha(2A)-adrenoceptors suggest that the second extracellular loop (XL2) folds above the binding cavity and may participate in antagonist binding.

EXPERIMENTAL APPROACH

Amino acids facing the binding cavity were identified using molecular models: side chains of residues 5.43 in TM5 and xl2.49 and xl2.51 in XL2 differ between the mouse and human receptors. Reciprocal mutations were made in mouse and human alpha(2A)-adrenoceptors at positions 5.43, xl2.49 and xl2.51, and tested with a set of thirteen chemically diverse ligands in competition binding assays.

KEY RESULTS

Reciprocal effects on the binding of yohimbine and rauwolscine in human and mouse alpha(2A)-adrenoceptors were observed for mutations at 5.43, xl2.49 and xl2.51. The binding profile of RS-79948-197 was reversed only by the XL2 substitutions.

CONCLUSIONS AND IMPLICATIONS

Positions 5.43, xl2.49 and xl2.51 are major determinants of the species preference for yohimbine and rauwolscine of the human versus mouse alpha(2A)-adrenoceptors. Residues at positions xl2.49 and xl2.51 determine the binding preference of RS-79948-197 for the human alpha(2A)-adrenoceptor. Thus, XL2 is involved in determining the species preferences of alpha(2A)-adrenoceptors of human and mouse for some antagonists.