Conopeptide ρ-TIA Defines a New Allosteric Site on the Extracellular Surface of the α1B-Adrenoceptor*♦

Background: Mechanistic insight into allosteric modulation of GPCRs can facilitate antagonist design. Results: Extracellular surface residues (ECS) of the α1B-adrenoceptor at the base of extracellular loop 3 interact with the allosteric antagonist TIA. Conclusion: The identified ECS pharmacophore provides the first structural constraints for allosteric antagonist design at α1-adrenoceptors. Significance: Binding to the ECS of a GPCR can allosterically inhibit agonist signaling. The G protein-coupled receptor (GPCR) superfamily is an important drug target that includes over 1000 membrane receptors that functionally couple extracellular stimuli to intracellular effectors. Despite the potential of extracellular surface (ECS) residues in GPCRs to interact with subtype-specific allosteric modulators, few ECS pharmacophores for class A receptors have been identified. Using the turkey β1-adrenergic receptor crystal structure, we modeled the α1B-adrenoceptor (α1B-AR) to help identify the allosteric site for ρ-conopeptide TIA, an inverse agonist at this receptor. Combining mutational radioligand binding and inositol 1-phosphate signaling studies, together with molecular docking simulations using a refined NMR structure of ρ-TIA, we identified 14 residues on the ECS of the α1B-AR that influenced ρ-TIA binding. Double mutant cycle analysis and docking confirmed that ρ-TIA binding was dominated by a salt bridge and cation-π between Arg-4-ρ-TIA and Asp-327 and Phe-330, respectively, and a T-stacking-π interaction between Trp-3-ρ-TIA and Phe-330. Water-bridging hydrogen bonds between Asn-2-ρ-TIA and Val-197, Trp-3-ρ-TIA and Ser-318, and the positively charged N terminus and Glu-186, were also identified. These interactions reveal that peptide binding to the ECS on transmembrane helix 6 (TMH6) and TMH7 at the base of extracellular loop 3 (ECL3) is sufficient to allosterically inhibit agonist signaling at a GPCR. The ligand-accessible ECS residues identified provide the first view of an allosteric inhibitor pharmacophore for α1-adrenoceptors and mechanistic insight and a new set of structural constraints for the design of allosteric antagonists at related GPCRs.

The G protein-coupled receptor (GPCR) superfamily is an important drug target that includes over 1000 membrane receptors that functionally couple extracellular stimuli to intracellular effectors. Despite the potential of extracellular surface (ECS) residues in GPCRs to interact with subtype-specific allosteric modulators, few ECS pharmacophores for class A receptors have been identified. Using the turkey ␤ 1 -adrenergic receptor crystal structure, we modeled the ␣ 1B -adrenoceptor (␣ 1B -AR) to help identify the allosteric site for -conopeptide TIA, an inverse agonist at this receptor. Combining mutational radioligand binding and inositol 1-phosphate signaling studies, together with molecular docking simulations using a refined NMR structure of -TIA, we identified 14 residues on the ECS of the ␣ 1B -AR that influenced -TIA binding. Double mutant cycle analysis and docking confirmed that -TIA binding was dominated by a salt bridge and cationbetween Arg-4--TIA and Asp-327 and Phe-330, respectively, and a T-stacking-interaction between Trp-3--TIA and Phe-330. Water-bridging hydrogen bonds between Asn-2--TIA and Val-197, Trp-3--TIA and Ser-318, and the positively charged N terminus and Glu-186, were also identified. These interactions reveal that peptide binding to the ECS on transmembrane helix 6 (TMH6) and TMH7 at the base of extracellular loop 3 (ECL3) is sufficient to allosterically inhibit agonist signaling at a GPCR. The ligand-accessible ECS residues identified provide the first view of an allosteric inhibitor pharmacophore for ␣ 1 -adrenoceptors and mechanistic insight and a new set of structural constraints for the design of allosteric antagonists at related GPCRs.
G protein-coupled receptors (GPCRs) 2 play a major role in physiological responses to hormones and neurotransmitters, in addition to being responsible for vision, olfaction, and taste (1). The initial GPCR x-ray structure of bovine rhodopsin determined by Palczewski et al. (2) provided the first useful template for building homology models of mammalian GPCRs. However, the low level of sequence homology with other GPCR families and the uncertainty in aligning the variable loop regions have limited structure-based drug design using the bovine rhodopsin derived models (3). Recent x-ray structures of bovine opsin (4,5), human ␤ 2 -adrenergic (6), turkey ␤ 1 -adrenergic (7), human A 2A -adenosine (8), human dopamine D3 (9), human muscarinic M2 (10), and rat M3 (11) receptors have helped overcome many of these limitations associated with modeling class A GPCRs (3).
As the major mediators of smooth muscle contraction, the ␣ 1 -adrenoceptors are important targets for the treatment of hypertension and benign prostatic hyperplasia. ␣ 1 -Adrenoceptors are also present in neurons in the brain and spinal cord, as well as in the liver (12). Three ␣ 1 -adrenoceptor (␣ 1A -, ␣ 1B -, and ␣ 1D -AR) subtypes have been characterized both pharmacologically and at the gene level. In human tissue, the ␣ 1A -AR is the most abundant subtype in liver, heart, and cerebral cortex, with lower expression levels in bladder; the ␣ 1B -AR is also expressed in liver, heart, and cerebral cortex but not in bladder, although the ␣ 1D -AR is the main subtype expressed in bladder with lower expression levels in other tissues (12,13). However, their relative roles are poorly understood given the lack of selective ␣ 1 -AR antagonists that can distinguish between the different ␣ 1 -subtypes. This lack of discrimination is not surprising because current small molecule inhibitors act competitively at the relatively conserved norepinephrine (NE) binding pocket and fail to take advantage of the potential offered by more struc-* This work was supported by National Health and Medical Research Council turally divergent allosteric sites. Recent NMR studies of the ␤ 2 -adrenoceptor have revealed that agonists and antagonists stabilize distinct GPCR conformations, demonstrating conformational coupling between the extracellular surface (ECS) and the orthosteric binding site (14). However, GPCR structure-function studies have focused on relatively conserved 7TM residues, and the role of residues that form the ECS remains poorly defined, despite their potential to allosterically modulate GPCR signaling.
-TIA, a conopeptide from the piscivorous Conus tulipa, acts as an allosteric inhibitor of the ␣ 1 -AR (15,16), with structureactivity studies identifying Arg-4 and Trp-3 as the key contributors to -TIA affinity for the hamster ␣ 1B -AR (17). Given that -TIA is a small and highly positively charged peptide that is unlikely to cross cell membranes, we systematically mutated residues associated with the ECS of the ␣ 1B -AR. Using a refined NMR structure of -TIA and a turkey ␤ 1 -AR-derived homology model of the ␣ 1B -AR to guide docking simulations and mutagenesis studies, we established the specific molecular interactions between -TIA on the ECS of the ␣ 1B -AR that were confirmed by double mutant cycle experiments. These findings, which reveal that peptide binding to ECS residues on TMH6 and TMH7 at the base of ECL3 can inhibit ␣ 1 -AR signaling, may prove useful in guiding the design of new classes of allosteric inhibitors.
Transient Expression of ␣ 1B -ARs and Membrane Preparation-COS-1 cells (ATCC) were cultured in Dulbecco's modified Eagle's medium (DMEM) containing glutamine and 5% fetal bovine serum. Cells were transiently transfected with purified plasmid DNA encoding WT or mutant ␣ 1B -ARs, using Lipofectamine 2000 (Invitrogen) (25 g of DNA/145 cm 2 ) or FuGENE HD (Roche Applied Science) (36 g of DNA/145 cm 2 ), following the manufacturers' protocol. Cell membranes were prepared 48 h post-transfection as follows. Cells were harvested and homogenized using a Polytron homogenizer (Brinkmann Instruments) in HEM buffer (20 mM HEPES, 1.5 mM EGTA, 12.5 mM MgCl 2 , pH 7.5) containing complete protease inhibitor (Roche Diagnostics). The homogenate was centrifuged at 2000 rpm for 10 min, and the resulting supernatant was centrifuged at 14,000 rpm for 30 min. The pellet was resuspended in HEM buffer containing 10% v/v glycerol and stored at Ϫ80°C prior to use. Protein concentration was determined using the BCA protein assay kit (Pierce) following the manufacturer's protocol.
Radioligand Binding Assays-The affinity of -TIA and analogs at the ␣ 1B -AR mutants were determined using the radiolabeled ␣ 1 -AR antagonist [ 3 H]prazosin (0.5 nM) or 125 I-(Ϯ)␤-(iodo-4-hydroxyphenyl)ethylaminomethyl-tetralone ( 125 I-HEAT; 70 pM). Reactions containing radioligand, membranes from ␣ 1B -ARtransfected COS-1 cells (5 g of protein), and increasing concentrations of -TIA or analogs (10 pM-10 M) in HEM buffer were established in clear round bottom 96-well plates. Saturation binding experiments were also performed to determine the K d value for prazosin at each of the mutants. Nonspecific binding was determined in the presence of 10 M phentolamine. Each assay was performed in triplicate in a total reaction volume of 150 l. After incubation for 60 min at room temperature, the membranes were harvested onto Whatman GF/B filter mats (PerkinElmer Life Sciences) pretreated with 0.6% polyethyleneimine using a TomTec harvester. Beta Plate scintillant (PerkinElmer Life Sciences) was then applied, and the filter-bound radioactivity detected using a Wallac MicroBeta (PerkinElmer Life Sciences).
IP 1 HTRF Assay-COS-1 cells (ATCC) were cultured and transiently transfected with WT or mutant ␣ 1B -AR DNA following the manufacturer's protocol (Lipofectamine 2000, Invitrogen or FuGENE HD, and Roche Applied Science). Assays measuring IP 1 accumulation were performed 48 h post-transfection following the manufacturer's instructions (IPone HTRF assay kit, Cisbio International). In brief, increasing concentrations of NE, prazosin, or -TIA were added to 30,000 transfected cells in stimulation buffer in a white 384-well plate (Optiplate, PerkinElmer Life Sciences). The plates were incubated for 1 h at 37°C with 5% CO 2 . Cells were then lysed by the addition of HTRF reagents, the europium cryptate-labeled anti-IP 1 antibody, and the d 2 -labeled IP 1 analog and diluted in lysis buffer (IPone HTRF assay kit, Cisbio International), followed by incubation for 1 h at room temperature. The emission signals were measured at 590 and 665 nm after excitation at 340 nm using the Envision multilabel plate reader (PerkinElmer Life Sciences), and the signal was expressed as the HTRF ratio: F ϭ ((fluorescence 665 nm /fluorescence 590 nm ) ϫ 10 4 ). The WT ␣ 1B -AR was routinely run each day to control for any batch-tobatch variation in response. In addition, antagonist-induced inhibition of the basal IP 1 mediated by the constitutive active mutant E289K was also determined (three separate experiments performed in replicates of six).
Synthetic Strategy--TIA and its analogs are peptide amides and were synthesized on a Protein Technologies Symphony automated peptide synthesizer using Rink amide resin (0.1 mmol). Assembly of the peptides was performed using 2-(1Hbenzotriazol-1-yl)-1,1,3,3-tetramethyluronium hexafluorophosphate/DIEA in situ activation protocols (19) to couple the Fmoc-protected amino acid to the resin (5 eq excess, coupling time 5 min). Fmoc deprotection was performed with 30% piperidine/DMF for 1 min followed by a 2-min repeat. Washes were performed 10 times after each coupling as well as after each de-protection step. Acetylation of -TIA was performed by using a 10-fold excess of acetic acid anhydride/DIEA in DMF. The success of the acetylation was monitored by the quantitative ninhydrin test (20). After chain assembly and final Fmoc deprotection, the peptide resins were washed with methanol and dichloromethane and dried in a stream of N 2 . The cleavage of peptide resin was performed at room temperature (RT) in TFA/ H 2 O/triisopropyl silane/ethane dithiol (87.5:5:5:2.5) for 3 h. Cold diethyl ether (30 ml) was then added to the filtered cleavage mixture, and the peptide was precipitated. The precipitate was collected by centrifugation and subsequently washed with further cold diethyl ether to remove scavengers. The final product was dissolved in 50% aqueous acetonitrile and lyophilized to yield a white solid product. The crude reduced peptide was examined by reversed phase-HPLC for purity and the correct molecular weight confirmed by electrospray mass spectrometry (ES-MS).
Disulfide Bond Formation of Peptides 1-5-Pure reduced peptides (1 mg/ml) were oxidized by stirring at RT in 30% isopropyl alcohol, 0.1 M NH 4 HCO 3 at pH 7.0 for 16 h. The solutions were subsequently diluted to an isopropyl alcohol concentration Ͻ5%, prior to reversed phase-HPLC purification. The obtained single main oxidized product with the cysteine connectivity (1-3, 2-4) (17) was purified to Ͼ95% purity and lyophilized.
HPLC Analysis and Purification-Analytical HPLC runs were performed using a Shimadzu HPLC system LC10A with a dual wavelength UV detector set at 214 and 254 nm. A reversed phase C-18 column (Zorbax 300-SB C-18; 4.6 ϫ 50 mm) with a flow rate of 2 ml/min was used. Gradient elution was performed with the following buffer systems: A, 0.05% TFA in water, and B, 0.043% TFA in 90% acetonitrile in water, from 0% B to 80% B in 20 min. The crude peptides and oxidized peptides were purified by semi-preparative HPLC on a Shimadzu HPLC system LC8A associated with a reversed phase C-18 column (Vydac C18, 25 ϫ 10 mm) at a flow rate of 5 ml/min with a 1%/min gradient of 5% B to 50% B. The purity of the final product was evaluated by analytical HPLC (Zorbax 300SB C-18: 4.6 ϫ 100 mm) at a flow rate of 1 ml/min and a 1.67%/min gradient of B (5-45%). Purities of synthesized peptides were all Ͼ95%.
Peptide Concentration-Peptide concentrations used in in vitro screening were calculated based on peak size detected at 214 nm by HPLC. Peak size was calibrated using a peptide standard (in this case 1) with known peptide content established by amino acid analysis. Molecular extinction coefficients are calculated for the standard and the peptide of interest applying increments established by Buck et al. (21). Using Lambert Beer Law, the peptide concentration was calculated based on absorptions of standard and sample using calculated extinction coefficients. Quantitation was confirmed using averaged IC 50 values for -TIA displacement of [ 3 H]prazosin determined for each batch of -TIA.
ES-MS-Electrospray mass spectra were collected in line during analytical HPLC runs on an Applied Biosystems API-150 spectrometer operating in the positive ion mode with a declustering potential of 20, a focusing potential of 220, and Turbo spray of 350°. Masses between 300 and 2200 atomic mass units were detected (step 0.2 atomic mass units, dwell 0.3 ms).
NMR Analysis-Samples of -TIA prepared for NMR analysis included 1 mg/ml peptide in 95% H 2 O, 5% D 2 O at pH 4.5. All data, including two-dimensional TOCSY, NOESY, and DQF-COSY, were acquired on a Bruker Avance 600 MHz spectrometer at 298 K with 4k data points in the direct dimension and 512 increments in the indirect dimension (22). The data were processed using Topspin 2.1 with the free induction decay multiplied by a 90°phase-shifted sine square window function, and zero filled to 1k data points in the indirect dimension prior to Fourier transform. The data were analyzed using the program CARA and assigned using standard two-dimensional homonuclear sequential assignment methods. Inter-proton distance restraints were derived from cross-peak intensities in a NOESY spectrum recorded with a mixing time of 250 ms. Backbone dihedral angles were derived from 3 J H␣HN coupling constants, with residues Ͻ5 Hz (Trp-3, Arg-4, Cys-5, Ala-10, Arg-12, Arg-13, Lys-16, Lys-17, and Cys-19) being restrained to Ϫ60 Ϯ 30°a nd residues Ͼ8 Hz (Asn-2, Leu-7, Ile-8, Cys-11, Asn-14, and His-15) being restrained to Ϫ120 Ϯ 30°. Side chain 1 angles and stereo specific assignments of H␤ protons were derived from a combination of 3 J H␣H␤ coupling constants and H␣-H␤ and HN-H␤ NOESY peak intensities. Residues found to be in a t2g3 conformation (Arg-4, Cys-5, Cys-6, Leu-7, His-15, and Cys-19) were restrained to Ϫ60 Ϯ 30°, and residues in a g2t3 conformation (Cys-11 and Arg-12) were restrained to Ϫ180 Ϯ 30°. No residues could be confirmed to be in a g2g3 conformation. The amide protons of Cys-5, Cys-6, Leu-7, and Ile-8, Cys-11, and Cys-19 have previously been reported to be slow exchanging. During preliminary structure calculations, backbone carbonyl groups acting as hydrogen bond acceptors were unambiguously identified for all these amides, and thus restraints for these hydrogen bonds could be introduced in the calculations (Cys-5 3 Asn-2, Cys-6 3 Trp-3, Leu-7 3 Arg-4, Ile-8 3 Cys-5, Cys-11 3 Ile-8, and Cys-19 3 His-15). Structures were calculated using restrained simulated annealing within the program CNS using the topallhdg5.3.pro/ parallhdg5.3.pro force field. Initial structures were calculated using torsion angle dynamics with a high temperature phase comprising 4000 steps of 0.015 ps at 50,000 K and a cooling phase during which the temperature was lowered from 50,000 to 0 K through 4000 steps of 0.015 ps, followed by 5000 steps of Powell energy minimization. In a second step, these structures were refined through Cartesian dynamics in a water shell. This refinement involves the following: (i) a heating phase where the temperature was increased to 500 K in steps of 100 K, each comprising 50 steps of 0.005-ps dynamics; (ii) a high temperature phase comprising 2500 steps of 0.005-ps dynamics at 500 K; (iii) a cooling phase where the temperature was lowered in steps of 100 K, each comprising 2500 steps of 0.005-ps dynamics, and finally (iv) an energy minimization phase comprising 3000 steps of Powell minimization. In the final round, 50 structures were calculated, and the 20 structures with lowest overall energy were chosen to represent the solution structure of -TIA. These structures were analyzed and visualized using the programs MolProbity (23) and MOLMOL (24). The coordi-nates of the structure and chemical shifts have been submitted to the PDB and BMRB databases and given the ID numbers 2LR9 and RCSB102733, respectively.
Modeling and Docking Simulations-A molecular homology model of the inactive ␣ 1B -AR was built using the crystal structure of turkey ␤ 1 -adrenoceptor (PDB code 2VT4 (7)) as the template in the program Modeler 9 Version 2 (25). This template was chosen because of its high degree of sequence identity to ␣ 1B -AR as determined by BLAST (26) and PHYRE (27). All sequence alignments were generated using Clust-alW (28) and further corrected by hand based on the secondary structure predicted by PHYRE (27) to maximize the  4DAJ). The three-dimensional structural alignments were performed using program TOPOFIT (55) with a joint distance cutoff of 3 Å between the C␣ atoms of the residues in the ␣ 1B -AR model and each crystal structure. The alignment represents the degree of conservation of residues spatially. The poorly aligned regions indicate residues are not topologically close, although the sequence may be conserved. accuracy of sequence alignment and protein sequence identity prior to modeling (Fig. 1). The structural models were validated using the on-line server Verify3D (29) and the Ramachandran plot available from ProFunc (30) database.
Blind (global) docking was applied initially to search for possible binding sites and used to assess the quality of our homology model by comparing the literature with our experimental data. In blind docking simulation, all side chains on -TIA (PDB code 2LR9), prazosin (DrugBank), and the entire ECS on ␣ 1B -AR were set as active residues using Haddock (31). The passive residues, which are solvent-exposed and surround the active residues, were defined by Haddock. No inter-or intramolecular restraints were defined. Semi-flexible and fully flexible regions were automatically detected by Haddock, and sampling parameters were set as default. The blind docking solutions for -TIA and prazosin were examined manually in PyMOL to define the predicted ligandbinding sites. The predicted binding sites for -TIA and prazosin from blind docking were in agreement with our experimental data, suggesting that the ␣ 1B -AR homology model was predictive.
Flexible solvated docking simulations were performed using Haddock (31), with experimentally determined restraints used to refine the binding mode. Specifically, Trp-3 and Arg-4 on -TIA (17) and Ser-318, Asp-327, and Phe-330 on ␣ 1B -AR were set as active residues, and the neighboring solvent-exposed residues were set as passive restraints and semi-and fully flexible fragments were automatically detected by Haddock. Sampling and cluster parameters were default values, except that 2000 structures were calculated for rigid body docking, with the best 400 solutions then subjected to semi-flexible simulated annealing and explicit solvent refinement. The Haddock score, minimal energy, existing experimental data, and the multiplicity of the same docking orientation (r.m.s.d. Ͻ7.5 Å) were used as selection criteria for choosing possible binding modes. The final 200 best solutions generated by Haddock were confirmed by visual inspection using PyMOL. Comparing the three approached solutions derived from Haddock provided the best fit for the available -TIA-␣ 1B -AR structureactivity data. Slight outward conformational movements at TMH6 and TMH7 and ECL3 (1.5-2.2 Å) were observed in the model following -TIA docking, although the structure of -TIA was little altered upon docking (Fig. 2). Additional ␣ 1B -AR mutants (D327A/F330A and S318A/F330A) were modeled following residue replacement and docked against -TIA and its W3A and R4A analogs, as described above, to provide additional support for the predicted intermolecular interactions.
Statistics and Data Analysis-Data are presented as means Ϯ S.E. of results obtained from 2 to 52 separate experiments, as indicated. Sigmoidal curves for the calculation of the half-maximal excitatory concentration (EC 50 ) and half-maximal inhibitory concentration (IC 50 ) values were fitted to individual data points by nonlinear regression using the software package Prism (GraphPad Software). For multiple comparisons, one-way analysis of variance was used with post hoc t tests performed by Dunnett's method. Values of p Ͻ 0.05 were considered significant.

Characterization of -TIA and Prazosin Affinity at ECS
Mutants of the ␣ 1B -AR-Given -TIA is a relatively hydrophilic noncompetitive inhibitor, we hypothesized that this peptide does not partition into or across the cell membrane and likely binds to ECS residues above and separate from the orthosteric site that is deeper in the transmembrane helical bundle of ␣ 1 -ARs. To determine the contribution of individual ECS residues to -TIA binding, we constructed 53 alanine-substituted point mutations of the ␣ 1B -AR across ECL1-3 and the extracellular portions of associated transmembrane helices ( Fig. 3 and Table 1). Wild type and mutant ␣ 1B -ARs were expressed in COS-1 cells and first tested for their ability to bind the inverse agonist prazosin. Saturation binding experiments generated prazosin K d values for all mutants except E186A and W307A (Fig. 3A). The E186A mutant showed readily detectable 125 I-HEAT binding allowing any potency shift for -TIA to be determined, whereas the W307A mutant failed to bind detectable FIGURE 2. Structural comparison of ␣ 1B -AR molecular models and template turkey ␤ 1 -AR (PDB code 2VT4). A, superimposition of ␣ 1B -AR model (white) to the template turkey ␤ 1 -AR crystal structure (blue) (r.m.s.d. ϭ 0.334 Å), showing that ECL2 is the only region that significantly differs from the template. B, comparison of -TIA-bound (yellow) and prazosin-bound (magenta) conformations to the initial ␣ 1B -AR model (white), showing that the orthosteric ligand binding does not perturb the overall structure, which is in agreement with the literature. Allosteric peptide binding introduces slightly conformational movements (Ͻ1Å between C␣ atoms) in TMH6, TMH7, ECL2, and ECL3 extracellularly and in TMH5 intracellularly. C, superimposition of the docked -TIA (orange) to its NMR structures (green) shows that the conformation of -TIA was little altered during docking. levels of either 125 I-HEAT or [ 3 H]prazosin. Saturation binding experiments indicated that the G109A and G114A mutants increased the affinity of prazosin 9-and 8-fold, respectively, compared with that of WT (K d ϭ 0.76 Ϯ 0.05 nM; n ϭ 29), whereas the remaining mutations either did not significantly affect or slightly decreased (4-, 5-, and 7-fold changes for mutants G183A, F310A, and F311A, respectively) the K d value for prazosin (Fig. 3A). In contrast, -TIA affinity was significantly affected by 14 receptor mutants (Fig. 3B). Ten mutants V107A (4-fold), G109A (7-fold), Y110A (9-fold), G114A (5-fold), C118A (60-fold), C195A (70-fold), F310A (6-fold), F311A (15-fold), K324A (6-fold), and K331A (6-fold) enhanced affinity, whereas four mutants G183A (4-fold), S318A (3-fold), D327A (4-fold), and F330A (10-fold) reduced affinity compared with the WT receptor (IC 50 ϭ 27.8 Ϯ 1.1 nM; n ϭ 52) (Fig. 3B).
Inverse Agonist Activity of -TIA-To further characterize the pharmacological properties of -TIA, we determined whether it acted as an inverse agonist or a neutral antagonist using the constitutively active E289K mutant of the ␣ 1B -AR (32). Using the IP 1 assay, we confirmed the E289K mutant was constitutively active (220% over basal) and inhibited both by prazosin (EC 50 ϭ 11.6 Ϯ 7.5 nM; n ϭ 3), a known inverse agonist (33), and -TIA (EC 50 ϭ 19.0 Ϯ 5.8 nM; n ϭ 3). In contrast, basal signaling was unaffected by prazosin or -TIA in WT ␣ 1B -AR.
Molecular Model of the ␣ 1B -Adrenoceptor-To better understand the mechanism of allosteric inhibition by -TIA, we constructed a homology model for the ␣ 1B -adrenoceptor using the crystal structure of the turkey ␤ 1 -adrenoceptor (7) as a template using the alignment shown in Fig. 1. An overlay of the ␣ 1B -AR model with turkey ␤ 1 -AR is shown in Fig. 2  were identified within ECL2, which is likely further stabilized by the conserved disulfide bond between Cys-195 and Cys-118, anchoring ECL2 to TMH3 (Fig. 6A). Remarkably, the alanine substitutions C118A or C195A improved -TIA affinity 60-and 70-fold, respectively, but produced a dramatic 286-and 1645fold reduction in the potency of NE to signal to IP 1 , while having no significant effect on prazosin binding (Figs. 3 and 4). Refinement of the Solution NMR Structure of -TIA-Although an initial solution NMR structure of -TIA has been reported (15), detailed analysis revealed inconsistencies between this structure and the NMR data that might compromise docking studies. Thus, we further refined the NMR structure of -TIA using an explicit solvent force field and additional restraints, including 1 dihedral angles. New NMR data were recorded at 600 MHz, and complete resonance assignments were achieved using standard homonuclear sequential assignment methods. In all cases, the resonance assignments were consistent with those previously reported (15). Structural restraints for -TIA included inter-proton distances, backbone and side chain dihedral angles, and hydrogen bonds, which were based on NOESY cross-peak intensities, coupling constants derived from one-dimensional 1 H and DQF-COSY spectra, and amide exchange data. These restraints were used as input for the generation of a refined -TIA NMR structure through sim-   ulated annealing and energy minimization in a water shell. The derived structural family shown in Fig. 5 (see Table 2 for structural statistics) is of high quality in terms of resolution of both backbone and side chains, and it is in excellent agreement with both experimental data and ideal covalent geometry, as evident from minimal violations and overall MolProbity scores of Ͼ90. The overall fold, which is dominated by two helical regions, is consistent with the previously reported structure (15), but some significant differences are seen in the N-terminal region critical for high affinity. In particular, in the new structure the Trp-3 and Arg-4 side chains stack against each other, which is consistent with the observed upfield shifts of the Arg-4 H␤ and H␥ resonances. Finally, the structure shows an improved correlation between the hydrogen bonding pattern, which includes a well defined 3 10 helical region between residues 3 and 7, and the amide exchange data. Based on these differences, we have used the new structure of -TIA (PDB code 2LR9) to establish its molecular complementarity with surface residues on the ␣ 1B -AR. Blind Docking of Prazosin-To validate our model, we performed a blind docking simulation of prazosin binding to our ␣ 1B -AR model. Prazosin was predicted to form extensive interactions with Trp-307, Phe-310, and Phe-311 in TMH6 and a strong hydrophobic interaction between the methyl tail with Phe-334 in TMH7 of the ␣ 1B -AR (Fig. 6, A  and B). The final docking pose chosen was consistent with these observations, the F310A and F311A mutants increased prazosin K d 6-and 7-fold (4.4 Ϯ 0.31, n ϭ 3, and 5.4 Ϯ 0.92 nM, n ϭ 3, respectively) compared with WT (0.76 Ϯ 0.05 nM, n ϭ 29) (Fig. 3), although the F334A mutant also slightly increased prazosin K d (2.1 Ϯ 0.22, n ϭ 3), although this change did not reach significance (Fig. 3A). Unfortunately, the W307A caused a dramatic drop in receptor expression levels (ϳ1% of WT) that precluded determination of prazosin affinity.
Experimentally Guided Docking of -TIA-Interactions identified from our mutational studies were used to guide the docking of -TIA with the program Haddock. This docking revealed that the binding site for -TIA was located on top of ECL3 and TMH7 (Fig. 7, A and B). The final docking pose for -TIA binding to the ␣ 1B -AR reveals it binds high on the ␣ 1B -AR, with the N terminus facing down toward ECL3 of the receptor (Fig. 7A).
-TIA binding is predicted to be dominated by a salt bridge between Arg-4--TIA and Asp-327, a cation-interaction between Arg-4--TIA and Phe-330, and a predicted T-stacking-interaction between Trp-3--TIA and Phe-330, consistent with our mutational and previous SAR studies (17). Three hydrogen bonds bridged by water molecules were also observed between N-terminal positive charge of -TIA and Glu-186, Asn-2--TIA and the backbone Val-197, and Trp-3--TIA and Ser-318. In this conformation, Asn-2--TIA forms a hydrogen bond with the backbone of Val-197, which could not be validated via mutational approaches. Surprisingly, ECL2 contributed little to -TIA binding, despite its high electronegativity that was expected to complement the positive charges on -TIA. Docking simulation with or without solvation yielded similar solutions (r.m.s.d. Ͻ 3 Å), except that the N-terminal positive charge of -TIA formed a direct ionic interaction with Glu-186 under nonsolvated conditions, whereas a weaker hydrogen bond through water to Glu-186 was predicted under solvated conditions. To better understand the role of water in allosteric binding, the N-terminal acetylated -TIA analog and the E186A mutant were constructed and tested. As Glu-199 is in close proximity to Glu-186 and may compensate a lost interaction in the E186A mutant, the E199A and E186A/E199A mutants were also constructed and tested. A significant 7-fold enhancement in affinity was obtained for acetylated -TIA at the ␣ 1B -AR (IC 50 ϭ 3.6 Ϯ 0.49 nM; n ϭ 3), whereas a significant 5-fold decrease in -TIA affinity was obtained at the E186A/ E199A mutant (IC 50 ϭ 139 Ϯ 22 nM; n ϭ 4). In contrast, acetylated -TIA at the E186A/E199A mutant had intermediate potency (IC 50 ϭ 9.4 Ϯ 2.2 nM; n ϭ 3), and -TIA affinity at mutants E186A (IC 50 ϭ 20.6 Ϯ 2.1 nM; n ϭ 4) and E199A (IC 50 ϭ 27.5 Ϯ 1.9 nM; n ϭ 3) was unchanged from WT receptor (IC 50 ϭ 27.8 Ϯ 1.1 nM; n ϭ 52). Together, these results indicate only a modest interaction between the N-terminal positive charge and Glu-186, supporting the solvated docking model that we used to generate the final docking pose for -TIA.

DISCUSSION
Allosteric ligands have the potential to be more selective modulators of class A GPCRs due to the structural divergence typically observed beyond the structurally and functionally conserved transmembrane region (3,34,35). In this study, we performed mutational and docking studies to determine the sites on the ␣ 1B -AR, where the orthosteric antagonist prazosin and the allosteric antagonist -TIA bind. Our studies establish that -TIA binds to the ECS at the base of ECL3 to inhibit GPCR signaling at a site that is distinct from the prazosin-binding site. Despite the differences in their binding sites, both -TIA and prazosin acted as inverse agonists at the ␣ 1B -AR.
Docking studies indicated that prazosin binds to the ␣ 1B -AR through extensive hydrophobic interactions with Trp-307, Phe-310, and Phe-311 in TMH6, and Phe-334 in TMH7. With the exception of Phe-334, which produced only a modest reduction in prazosin affinity, direct involvement of these resi-dues in prazosin binding was supported by mutational studies. The significant reduction in the affinity of prazosin for the F310A-␣ 1B -AR and F311A-␣ 1B -AR mutants was observed previously in the corresponding alanine mutants of closely related subtypes (36 -38). The W307A mutant, which forms part of the conserved CWXP motif and is proposed to function as an activation switch (39), reduced NE potency to signal to IP 1 but almost abolished specific prazosin binding, indicating this residue also contributes to prazosin binding. Because the hydrophobic residues identified as interacting with prazosin are highly conserved across all class A GPCRs, it appears that prazosin selectivity for ␣ 1 -ARs is determined by either the specific shape of this binding pocket or by adjacent residues in ECL2 and ECL3 that define the shape of the entrance to the orthosteric binding pocket and are highly divergent across class A GPCRs (3).
Based on its size, polarity, and noncompetitive pharmacology, -TIA was predicted to bind to the ECS at a site that was FIGURE 8. Docking of -TIA analogs and -TIA to the ␣ 1B -AR and mutants. The 10 lowest energy docking complexes for each simulation are shown. All panels are prepared in the same orientation for ease of comparison, with water molecules not shown for clarity. A, docking of W3A--TIA to ␣ 1B -AR showed ϳ25°rotation along the y axis (left and right) and/or z axis (middle), compared with that of -TIA in ␣ 1B -AR. An additional salt bridge is observed between the positively charged N termini of W3A--TIA and Glu-186 of ␣ 1B -AR. B, docking complexes of -TIA to ␣ 1B -AR-S318A/F330A showed a similar orientation to that of W3A--TIA to ␣ 1B -AR. A potential salt bridge is predicted between the positively charged N termini of W3A--TIA and Glu-186 or Glu-199 of ␣ 1B -AR. C, docking complexes of W3A--TIA and ␣ 1B -AR-S318A/F330A showed a different binding mode (ϳ180°rotation along the y axis) to that of -TIA in ␣ 1B -AR, with Arg-4 of -TIA-W3A forming a salt bridge with either Glu-199 or Glu-194, instead of Asp-327. D, R4A--TIA to ␣ 1B -AR showed rotation angles of ϳ20°and ϳ50°along the y and/or z axis, respectively, with the N terminus forming a salt bridge with either Glu-186 or Glu-194 on ␣ 1B -AR.
topologically distinct from the prazosin-binding site. Blind docking of -TIA to ␣ 1B -AR indicated a major docking cluster dominated by interactions with Asp-327. A refined docking simulation incorporating the structural restraints obtained from mutational studies and the SAR of -TIA (17) using Haddock revealed that -TIA formed a salt bridge with Asp-327 in TMH7, T-stacking-and cation-interactions with Phe-330 in TMH7, and a hydrogen bond with Ser-318 in TMH6. These interactions reveal that the allosteric binding site for -TIA is located on the ECS of the ␣ 1 -AR at the base of ECL3. These interactions were supported by mutational studies that confirmed -TIA and prazosin-binding sites were separate and distinct. Phe-330 is conserved across the ␣ 1 -ARs, whereas position 327 is either a negatively charged Asp in ␣ 1B -AR or Glu in ␣ 1Aand ␣ 1D -AR, explaining the shallow selectivity profile for -TIA at ␣ 1 -ARs (17). Interestingly, Phe-330 on the ␣ 1 -AR is at an equivalent position to the Trp in muscarinic receptors that contribute to allosteric small molecule inhibitor (40 -42) and three-fingered snake toxin MT7 (41) binding (Fig. 9), suggesting that this position may play a central role in the allosteric inhibition of GPCRs.
The identification of the intramolecular cationinteraction between Trp-3 and Arg-4 within -TIA, and the combination of intermolecular interactions between Trp-3 and Arg-4 on -TIA with Phe-330 and/or Asp-327 on ␣ 1B -AR complicate the interpretation of mutational data investigating the influence of these residues on -TIA affinity. For the double mutant S318A/ F330A, the interactions between ␣ 1B -AR and Trp-3--TIA were removed, resulting in a significant 65-fold drop in affinity. A comparable reduction in affinity (35-fold) was also observed for W3A--TIA at the WT ␣ 1B -AR, as observed previously (17). However, interactions between S318A/F330A and W3A--TIA decreased nearly 10-fold more than expected (590-fold) for a strictly complementary interaction. The adjacent Arg-4 in -TIA was also predicted to interact strongly with the ␣ 1B -AR through Asp-327 and Phe-330. Here, the single mutants D327A and F330A decreased -TIA affinity only 4-and 10-fold, respectively, although the double mutant D327A/F330A decreased the affinity ϳ280-fold, comparable with the ϳ440-fold loss in affinity for R4A-TIA at the WT ␣ 1B -AR. Again, R4A--TIA binding to the double mutant D327A/F330A caused a larger than expected drop in affinity (ϳ1845-fold). To address these discrepancies, docking studies of these mutant cycle experiments were performed that revealed ϳ20°shifts in the binding orientation of -TIA (Fig. 8). These larger than expected drops in -TIA affinity do not appear to arise from broader structural effects, because the double mutants E186A/E199A, S318A/ F330A, and D327A/F330A that affected -TIA affinity had no effect on the affinity of prazosin, with only the E186A/E199A mutant decreasing NE potency. Thus, we conclude that compensatory shifts in binding for the single mutant interactions resulted in an underestimate of the strength of individual binding interactions.
N-terminal truncated -TIA significantly (ϳ5-fold) reduced -TIA affinity for the ␣ 1B -AR (17), suggesting that the N-terminal positive charge might also contribute to -TIA binding. However, our solvated docking studies identified only a hydrogen bond bridged by water between Glu-186 and the N-terminal positive charge. This relatively weak hydrogen-bonding interaction was supported by studies on acetylated -TIA, which maintained WT or better affinity at E186A-␣ 1B -AR, the adjacent E199A-␣ 1B -AR mutant, and E186A/E199A-␣ 1B -AR double mutants.
In addition to identifying three closely positioned residues that contribute directly to -TIA binding, 11 more widely distributed residues were found that either enhanced the affinity of -TIA (V107A, G109A, Y110A, G114A, C118A, C195A, F310A, F311A, K324A, and K331A) or reduced -TIA affinity (Gly-183). These residues either had no effect, reduced (Gly-183, F310A, and F311A), or enhanced (Gly-109 and Gly-114) prazosin affinity and were predicted to play a structural role as they were topologically distant from the predicted allosteric binding site. To further characterize the effect of these mutations on the ␣ 1B -AR conformation, we determined their effect on NE signaling. V107A, Gly-109, K324A, and K331A had no effect, whereas Y110A, Gly-114, C118A, Gly-183, C195A, F310A, and F311A reduced NE potency to signal to IP 1 . Phe-310 has previously been reported to bind NE, and the F311A mutation has been suggested to cause global changes in receptor conformation, including folding and expression (36), consistent with the reduced expression observed for the Phe-311 mutant and reduced NE potency for the Phe-310 mutant. From our model, Lys-331 in TMH7 is predicted to form a salt bridge with Glu-106 in TMH2 and thus contribute to the ␣ 1B -AR structure by stabilizing TMH2 and TMH7. It appears that disrupting this salt bridge in ␣ 1B -AR favors conformations with enhanced -TIA affinity. However, the involvement of Lys-324 (ECL3) in -TIA binding remains unclear, and the increased -TIA affinity and decreased NE potency at the Y110A (ECL1) mutant cannot be easily explained from our model. Two glycine residues in ECL1 (Gly-109 and Gly-114) increased both -TIA and prazosin affinity when mutated to alanine. Interestingly, the G114A mutant reduced NE potency, whereas the G109A mutant had no effect. In contrast, the glycine mutant G183A in TMH4 showed decreased affinity for -TIA, prazosin, and NE, despite its distant position from the binding cavity. Because these five mutations are at positions distant from the binding pocket for -TIA and orthosteric ligands, it is likely that mutating these positions influences the structure of the ␣ 1B -AR to affect ligand binding at both the allosteric and orthosteric binding sites.
To better understand the role of extracellular surface residues in ligand binding and signaling, we compared the secondary structure of our model to the GPCR crystal structures determined (27). Although ECL2 varies widely across class A GPCR families, it invariably presents as a highly ordered hairpin-like conformation stabilized by a network of intramolecular interactions (4 -11, 43, 44) despite the lacking of an extensive well defined secondary structure. Consistent with this view obtained directly from x-ray crystallography, our homology model also showed a highly ordered hairpin-like conformation stabilized by a network of intramolecular interactions. The stability and integrity of ECL2 have been shown to affect ligand binding and receptor activation in many GPCRs by altering the molecular shape and size of the entrance of ligand-binding site (14,(45)(46)(47)(48)(49)(50)(51)(52)(53)(54). Interestingly, disrupting the disulfide bond between Cys-118 and Cys-195 significantly enhanced -TIA affinity ϳ60-fold, had no effect on prazosin affinity, and reduced NE potency ϳ300-fold. Surprisingly, [ 3 H]prazosin binding revealed that the C118A mutant had only ϳ1% of the expression level of the C195A mutant, indicating that these positions are not equivalent. A decrease in agonist potency of up to 4400-fold has been observed previously in other class A GPCRs when this disulfide bond was disrupted by mutations or exposure to reducing agents (45,(47)(48)(49)(50)(51)(52)(53)(54).
In conclusion, this study provides the first analysis of an allosteric binding site on ␣ 1 -ARs. Through extensive pharmacological characterization and computational simulations, we determined that -TIA binds to an allosteric modulatory site on TMH6 and TMH7 at the base of ECL3. -TIA binding was dominated by a salt bridge and cation-interaction between Arg-4--TIA and Asp-327 and Phe-330 in TMH7, respectively, and a hydrogen bond and T-stacking-interaction between Trp-3--TIA and Phe-330 in TMH7. The ECS pharmacophore identified in this study provides the first set of structural constraints for the design of novel selective allosteric antagonists acting at the ␣ 1 -ARs. A number of ECS mutations also affected NE potency, confirming the important role played by the ECS in regulating agonist signaling. The D(E)XXF diad in the ␣ 1B -AR varies across related GPCRs (Fig. 9), being mostly charged at the first position (DEKRH) and typically an aromatic residue at the second position (FWY), revealing opportunities to exploit this allosteric site for the development of new classes of subtype-selective allosteric inhibitors of GPCRs.