Abstract
Binding of the vasodilator peptides adrenomedullin (AM) and calcitonin gene–related peptide (CGRP) to the class B G protein–coupled receptor calcitonin receptor–like receptor (CLR) is modulated by receptor activity–modifying proteins (RAMPs). RAMP1 favors CGRP, whereas RAMP2 and RAMP3 favor AM. Crystal structures of peptide-bound RAMP1/2-CLR extracellular domain (ECD) heterodimers suggested RAMPs alter ligand preference through direct peptide contacts and allosteric modulation of CLR. Here, we probed this dual mechanism through rational structure-guided design of AM and CGRP antagonist variants. Variants were characterized for binding to purified RAMP1/2-CLR ECD and for antagonism of the full-length CGRP (RAMP1:CLR), AM1 (RAMP2:CLR), and AM2 (RAMP3:CLR) receptors. Short nanomolar affinity AM(37–52) and CGRP(27–37) variants were obtained through substitutions including AM S45W/Q50W and CGRP K35W/A36S designed to stabilize their β-turn. K46L and Y52F substitutions designed to exploit RAMP allosteric effects and direct peptide contacts, respectively, yielded AM variants with selectivity for the CGRP receptor over the AM1 receptor. AM(37–52) S45W/K46L/Q50W/Y52F exhibited nanomolar potency at the CGRP receptor and micromolar potency at AM1. A 2.8-Å resolution crystal structure of this variant bound to the RAMP1-CLR ECD confirmed that it bound as designed. CGRP(27–37) N31D/S34P/K35W/A36S exhibited potency and selectivity comparable to the traditional antagonist CGRP(8–37). Giving this variant the ability to contact RAMP2 through the F37Y substitution increased affinity for AM1, but it still preferred the CGRP receptor. These potent peptide antagonists with altered selectivity inform the development of AM/CGRP-based pharmacological tools and support the hypothesis that RAMPs alter CLR ligand selectivity through allosteric effects and direct peptide contacts.
Introduction
Adrenomedullin (AM) and calcitonin gene–related peptide (CGRP) are related vasodilator peptides that have overlapping and distinct functions in human physiology and pathophysiology (Russell et al., 2014; Kato and Kitamura, 2015; Klein and Caron, 2015). AM signaling is critical during development, as mice lacking AM or its receptor components were embryonic lethal due to cardiovascular and lymphatic defects (Caron and Smithies, 2001; Shindo et al., 2001; Dackor et al., 2006; Fritz-Six et al., 2008). AM is produced by endothelial and vascular smooth muscle cells. Circulating AM levels are elevated in various diseases, including myocardial infarction and heart failure, where it has cardioprotective effects. Additional AM functions include lymphatic maintenance and anti-inflammatory activity. CGRP is a neuropeptide released from perivascular sensory nerves where it acts as a microvascular vasodilator (Russell et al., 2014). CGRP has cardioprotective effects and is involved in migraine headache pathogenesis through its unique role in the trigeminovascular system (Russell et al., 2014; Russo, 2015; Edvinsson, 2017).
There is considerable interest in developing therapeutics targeting AM and CGRP signaling. Agonism of AM signaling holds promise for cardiovascular disorders, lymphedema, and inflammatory bowel disease, and antagonism of AM signaling may be of value for cancer (Kato and Kitamura, 2015; Klein and Caron, 2015). Agonism of CGRP signaling shows promise for cardiovascular disease (Aubdool et al., 2017), and antagonism of CGRP signaling by small-molecule antagonists of the CGRP receptor or monoclonal antibodies against the CGRP peptide or its receptor showed significant benefit in clinical trials for migraine (Edvinsson, 2015; Tso and Goadsby, 2017). Nonetheless, our understanding of AM and CGRP signaling is still limited because AM and CGRP share three highly similar heterodimeric receptors. These are composed of a common class B G protein–coupled receptor subunit, the calcitonin receptor–like receptor (CLR), and one of three variable receptor activity–modifying protein (RAMP1–3) subunits that alter CLR ligand selectivity (McLatchie et al., 1998; Hay and Pioszak, 2016). RAMP1 confers a preference for CGRP, and its complex with CLR is designated the CGRP receptor. RAMP2 and RAMP3 favor AM, and their complexes with CLR are designated the AM1 and AM2 receptors, respectively. Understanding the mechanistic basis for RAMP modulation of CLR pharmacology is crucial to harness the potential of this system for novel therapeutics.
CLR has an N-terminal extracellular domain (ECD) and a 7-transmembrane helical bundle domain. RAMPs have an N-terminal ECD and a single-transmembrane helix. AM and CGRP bind the heterodimeric receptor complexes through a “two-domain” model (Barwell et al., 2012). The C-terminal half of the peptide binds the ECD complex, and the N-terminal half binds and activates the CLR 7-transmembrane domain. N-terminally truncated peptides thus act as competitive antagonists (Chiba et al., 1989). CGRP(8–37) and AM(22–52) antagonists have traditionally been used to discriminate the CGRP and AM1/2 receptors, but these are only moderately selective. Antagonists that discriminate the AM1 and AM2 receptors are lacking. Consequently, it is not always clear which receptor mediates a given peptide effect.
The peptide selectivity profiles of the CGRP and AM1 receptors were recapitulated by purified soluble ECD heterodimers engineered as tethered RAMP ECD-CLR ECD fusion constructs (Moad and Pioszak, 2013). Crystal structures of tethered RAMP1-CLR ECD and RAMP2-CLR ECD with bound antagonist fragments of the high-affinity CGRP N31D/S34P/K35F variant (Rist et al., 1998) and AM, respectively, revealed that the peptides occupy a shared binding site on CLR and make limited contact with the RAMPs (Booe et al., 2015). The RAMPs augment a pocket occupied by the peptide C-terminal residue such that RAMP1 W84 contacts CGRP F37 and RAMP2 E101 hydrogen bonds with AM Y52. The two structures also revealed subtle conformational differences in CLR, which suggested a possible allosteric component to RAMP function.
We hypothesized that the mechanism of RAMP modulation of CLR ligand selectivity could be explored with AM and CGRP variants bearing amino acid substitutions designed to probe the putative RAMP allosteric effects and direct peptide contacts. Guided by the structures (Booe et al., 2015), we designed short high-affinity antagonist variants and then used these as scaffolds in which to add substitutions probing the RAMP mechanism. We characterized the binding of the variants to tethered RAMP-CLR ECD proteins and their antagonism at intact receptors, and used X-ray crystallography to confirm that an AM variant designed to have altered preference for the CGRP receptor bound as designed. Our results inform efforts to develop novel AM/CGRP-based pharmacological tools and therapeutics and support a dual mechanism for RAMP modulation of CLR ligand selectivity comprising both allosteric and direct peptide contact components.
Materials and Methods
Reagents.
General laboratory chemicals were from Sigma-Aldrich (St. Louis, MO). Dulbecco’s modified Eagle’s medium (DMEM) 4.5 g/l glucose, l-glutamine was from Lonza (Basel, Switzerland). Fatty-acid-free bovine serum albumin (BSA) was from PAA Laboratories. Fetal bovine serum was from Life Technologies (Carlsbad, CA). Sodium malonate, pH 6.0, was from Hampton Research (Aliso Viejo, CA). COS-7 cells were from American Type Culture Collection (Manassas, VA) (CRL-1651).
Plasmids.
All plasmids used the human RAMP and CLR sequences. The pETDuet-1 plasmids for bacterial expression of maltose-binding protein (MBP)-RAMP1.24-111-(GSA)3-CLR.29-144-H6 and MBP-RAMP1.24-111-(GS)5-CLR.29-144-H6 were previously described (Moad and Pioszak, 2013; Booe et al., 2015). The pETDuet-1–based plasmid encoding MBP-RAMP2.55-140[L106R]-(GS)5-CLR.29-144-H6 was constructed using Gibson Assembly (New England Biolabs, Ipswich, MA) by generating overlapping fragments to replace the (GS)2-(RS)-(GS)2 linker with (GS)5 in the previously reported plasmid MBP-RAMP2.55-140-(GS)2-(RS)-(GS)2-CLR.29-144-H6 (Moad and Pioszak, 2013). The L106R mutation in RAMP2 was introduced by site-directed mutagenesis using the QuikChange II kit (Agilent, Santa Clara, CA) as previously described (Booe et al., 2015). Primer sequences to generate Gibson Assembly fragments and for site-directed mutagenesis are available from the authors upon request. All plasmid sequences were verified by automated DNA sequencing by the University of Oklahoma Health Sciences Center Laboratory for Molecular Biology and Cytometry Research core facility. For transient transfections in COS-7 cells, expression plasmids containing untagged RAMP1-, RAMP2-, RAMP3-, or CLR-encoding sequences in the pcDNA3.1(+) vector were purchased from the Missouri University of Science and Technology cDNA resource center (Rolla, MO), and plasmid DNA was isolated using a Midi kit (Macherey-Nagel) according to the manufacturer’s instructions.
Purified Proteins.
The MBP-RAMP1.24-111-(GS)5-CLR.29-144-H6, MBP-RAMP1.24-111-(GSA)3-CLR.29-144-H6, and MBP-RAMP2.55-140[L106R]-(GS)5-CLR.29-144-H6 fusion proteins were expressed in Escherichia coli and purified as previously described (Moad and Pioszak, 2013; Booe et al., 2015). The numbers after RAMP1/2 and CLR indicate the amino acid residues included in the constructs. The RAMP2 L106R substitution prevents dimerization of the tethered construct; this mutation had no effect on the function of the intact AM1 receptor complex in cells (Booe et al., 2015). Purified proteins were dialyzed to storage buffer containing 50 mM Tris-HCl (pH 7.5), 50% (v/v) glycerol, and 150 mM NaCl and stored at −80°C. The control MBP-thrombin cut site-RSPO4-H6 fusion protein was previously described (Warner et al., 2015). Protein concentrations were determined by Bradford assay with a BSA standard curve, and the concentrations are stated in terms of the fusion protein monomer.
Synthetic Peptides.
High-performance liquid chromatography–purified custom synthetic peptides, including human AM(13–52) and all truncated CGRP and AM variants used in this study, were purchased from RS Synthesis (Louisville, KY). Human αCGRP(1–37) was purchased from Bachem (Bubendorf, Switzerland). Lyophilized peptides were reconstituted in sterile ultrapure water, and aliquots were stored at −80°C. Peptide concentrations were determined by UV absorbance at 280 nm with molar absorptivity calculated based on tyrosine, tryptophan, and cystine residues. Fluorescein isothiocyanate (FITC)-aminohexanoic acid (Ahx)-AM(37–52) S45W/Q50W concentration was determined by visible absorbance at 495 nm using the molar absorptivity of the FITC molecule (ε = 72,000 M−1 cm−1 at pH 8.0). The concentration of human αCGRP(1–37) was determined using the reported peptide content from amino acid analysis by Bachem. For peptides that lacked UV-absorbing residues, we assumed 80% peptide content for the lyophilized powder. The biotinylated peptides used in the AlphaLISA assays were previously described (Moad and Pioszak, 2013). All peptides were amidated (−CONH2) at their C terminus. The sequences of all peptides used in this study are provided in Supplemental Table S1.
AlphaLISA Luminescent Proximity Competition Peptide-Binding Assay.
This assay was performed as previously described (Lee et al., 2016), except that a 10-μg/ml final bead concentration was used and 0.3% (v/v) Triton X-100 was included in the reaction buffer for the MBP-RAMP2 ECD[L106R]-(GS)5-CLR ECD fusion protein. In brief, a master mix containing biotin-AM(19–52) C21A or biotin-αCGRP(8–37) and purified MBP-RAMP1 or −2 ECD-(GS)5-CLR ECD and streptavidin-coated donor and anti-MBP antibody–coated acceptor beads (PerkinElmer, Waltham, MA) was mixed with unlabeled competitor peptides at the indicated concentrations and incubated for 5 hours in the dark before reading luminescence in a Polarstar Omega plate reader equipped with AlphaLISA filters (BMG Labtech, Ortenberg, Germany). The competition binding curves were fit by nonlinear regression to the “log (inhibitor) versus response (three parameters)” equation in GraphPad Prism 5.0f (GraphPad Software, La Jolla, CA) to determine log IC50 values. No parameters were shared among data sets and no weighting or constraints were used.
Fluorescence Polarization/Anisotropy Peptide-Binding Assay.
This assay was performed as previously described (Lee et al., 2017), except using FITC-Ahx-AM(37–52) S45W/Q50W as the probe, and Tween-20 was omitted from the reaction buffer. Saturation binding used 7 nM FITC-Ahx-AM(37–52) S45W/Q50W and the indicated concentrations of MBP-RAMP1/2 ECD-(GS)5-CLR ECD fusion proteins with a 1-hour incubation at room temperature. Anisotropy values were corrected for the slight fluorescence intensity enhancement observed upon receptor binding (Supplemental Fig. S1E) as previously described (Lee et al., 2017). Competition binding assays used 7 nM FITC-peptide and 15 nM MBP-RAMP1 ECD-(GS)5-CLR ECD or 110 nM MBP-RAMP2 ECD-(GS)5-CLR ECD with the indicated concentrations of competitor peptides and a 2-hour incubation at room temperature. Anisotropy values were not corrected for the slight decrease in total fluorescence intensity observed upon probe displacement because determination of the competitor KI value is insensitive to minor intensity changes (Roehrl et al., 2004). Equilibrium dissociation constants were determined for the labeled probe from the saturation experiments (log KD) and for the unlabeled peptides from competition experiments (log KI) by nonlinear regression fitting of the binding curves to the exact analytical equations of Roehrl et al. (2004) using user-defined equations in GraphPad Prism version 5.0f as previously described (Lee et al., 2017). Roehrl et al. (2004) eqs. 6 and 39 were used for saturation binding, and eqs. 17 and 39 were used for competition binding. These equations are formulated in terms of the total probe, receptor, and competitor concentrations. Where applicable, the anisotropy values for the free and bound probe were constrained to be the same for multiple curves within the same experiment. In limited cases where a single competitor was assessed in an experiment and the full curve could not be obtained, the anisotropy value for the free probe (bottom of the curve) was constrained to be >0. No weighting was used. Fluorescence anisotropy detection used a PolarStar Omega plate reader equipped with fluorescence polarization (FP) optics and FITC filters (BMG Labtech).
cAMP Signaling Assay.
COS-7 cells (American Type Culture Collection) were cultured in DMEM with 10% fetal bovine serum and incubated at 37°C in a humidified 5% CO2 incubator. At 90%–100% confluency, cells were trypsinized, counted, and seeded in 96-well plates (Corning, Corning, NY) at 20,000 cells/well. Twenty-four hours after seeding, a plasmid encoding human RAMP1, RAMP2, or RAMP3 was cotransfected with a plasmid encoding CLR at a 1:1 ratio (125 ng/well RAMP1/2/3 and 125 ng/well CLR) using polyethyleneimine as previously described (Bailey and Hay, 2006). Forty-eight hours after transfection, the culture medium was aspirated, and cells were incubated with cAMP assay medium [serum-free DMEM supplemented with 1 mM IBMX and 0.1% (w/v) fatty-acid-free BSA] for 30 minutes at 37°C in the presence or absence of the indicated concentration of antagonist peptide. At the 30-minute time point, the cells were stimulated with the indicated concentrations of AM(13–52) or αCGRP(1–37) agonist in cAMP assay buffer in the absence or presence of antagonist at the same concentration used for the preincubation and incubated a further 30 minutes at 37°C. The medium was then aspirated, and the cells were lysed with ice-cold 6% perchloric acid. The plate was shaken on an orbital shaker for 15 minutes at room temperature, and the lysates were neutralized by addition of potassium bicarbonate and Na-HEPES, pH 7.4. cAMP in the neutralized lysates was quantitated using a LANCE cAMP kit (PerkinElmer) according to the manufacturer’s instructions with a 4-hour assay incubation time. LANCE TR-FRET detection used a PolarStar Omega plate reader equipped with an advanced TRF optic head and LANCE filters (BMG Labtech). The 665/620-nm emission ratio was used in the data analysis, and cAMP concentration was determined by interpolation from a cAMP standard curve. Concentration-response curves were fitted by nonlinear regression in GraphPad Prism version 5.0f to the “Gaddum/Schild EC50 shift” equation with the Schild and Hill slopes constrained to 1 to determine antagonist potency (apparent pKB). No weighting was used. Antagonist concentrations were chosen to give a significant shift of the concentration-response curve to enable accurate determination of the apparent pKB values.
Crystallization, X-Ray Diffraction Data Collection, Structure Solution, and Refinement.
Purified MBP-RAMP1.24-111-(GSA)3-CLR.29-144-H6 fusion protein was mixed with synthetic AM(37–52) S45W/K46L/Q50W/Y52F at a 1:1.2 protein:peptide molar ratio and concentrated to 30 mg/ml as previously described (Booe et al., 2015). Crystals were grown by the hanging drop vapor diffusion method in 15% (w/v) polyethylene glycol 3350, 0.1 M sodium malonate (pH 6.0), 50 mM potassium sodium tartrate, and 1% (w/v) cadaverine. The best crystals were obtained by microseeding using a Seed Bead kit (Hampton Research). Due to their fragile nature, the crystals were cross-linked before cryoprotection by incubating the drop for 15 minutes over 25% glutaraldehyde. The cross-linked crystals were directly transferred into a cryoprotection solution comprising crystallization solution +8% (v/v) 2-methyl-2,4-pentanediol and flash frozen in liquid nitrogen. Diffraction data were collected remotely at LS-CAT beamline 21-ID-G (λ = 0.9786 Å) at the Advanced Photon Source (Argonne, IL). A diffraction data set from a single crystal was processed using HKL2000 version 712 (HKL Research, Charlottesville, VA) (Otwinowski and Minor, 1997). The structure was solved by molecular replacement (MR) with Phaser version 2.6.0 (McCoy et al., 2007) using MBP with maltose removed (Protein Data Bank: 3C4M) and the RAMP1 ECD-CLR ECD heterodimer with MBP and peptide removed (Protein Data Bank: 4RWG) as search inputs. The MR solution contained four molecules of MBP and four molecules of RAMP1 ECD-CLR ECD heterodimer in the asymmetric unit. The MR solution was rigid body refined with Phenix.refine version 1.10.1-2155 (Adams et al., 2010) by treating MBP, RAMP1 ECD, and CLR ECD as rigid bodies. The model was completed by iterative rounds of manual model building using COOT (Emsley et al., 2010) and NCS and TLS restrained refinement using Phenix version 1.10.1-2155 (Adams et al., 2010). Automatic NCS restraints, stereochemistry weight optimization, and B-factor weight optimization were used.
Structural Analysis and Modeling of Peptide Variants.
Structural analysis was performed in PyMOL (Schrödinger). Structural superimpositions were performed using the align command based on Cα atom positions of CLR ECD. For modeling of peptide variants, in silico mutagenesis was performed with the mutagenesis wizard.
Statistical Analysis.
The binding and signaling assays were performed at least three independent times (on different days) with duplicate samples. Means from the independent replicates are reported as ±S.E.M. Statistical comparison of pKI values from the FP assays and apparent pKB values from the signaling assays were performed in GraphPad Prism version 5.0f. Unpaired two-tailed t test was used to compare the pKI values of an individual peptide variant at the RAMP1-CLR ECD and RAMP2-CLR ECD complexes. One-way analysis of variance followed by Tukey’s post-hoc test was used to compare the apparent pKB values for an individual peptide variant at the three intact receptor complexes. Similarly, analysis of variance with Tukey’s post-hoc test was used to compare the pKI or apparent pKB values of all AM or CGRP peptide variants at an individual receptor complex (e.g., RAMP1-CLR ECD or RAMP1:CLR). Significance was determined as P < 0.05, P < 0.01, or P < 0.001. All statistical comparisons are summarized in Supplemental Tables S2 and S3 and selected comparisons are shown in Tables 1 and 2 and in scatter plot format in Supplemental Figs. S4 and S7.
pKI values for AM and CGRP variants determined by FP peptide-binding assay
Apparent pKB values for AM and CGRP antagonist variants determined by cell-based signaling assay in COS-7 cells
Results
Rational Design of Short High-Affinity AM and CGRP Antagonist Variants.
The two-domain agonist-binding mechanism enables N-terminal peptide truncation to be used to create competitive antagonists such as the traditional AM(22–52) and CGRP(8–37) antagonists that have nanomolar affinity for their receptors (Fig. 1, A and G). Unfortunately, further truncation to the ECD-binding regions, AM(37–52) and CGRP(27–37), weakens affinity into the micromolar range (Moad and Pioszak, 2013). We first set out to develop high-affinity AM(37–52) and CGRP(27–37) variants because we reasoned that these would serve as ideal scaffolds in which to add additional amino acid substitutions designed to probe RAMP allosteric effects and direct peptide contacts.
Enhancing AM affinity through rational structure-based design. (A) “Two-domain” peptide-binding model for RAMP:CLR complexes. (B) Crystal structures of AM-bound MBP-RAMP2-CLR ECD (Protein Data Bank 4RWF) and CGRPmut-bound MBP-RAMP1-CLR ECD (PDB 4RWG) in cartoon representation. MBP is omitted for clarity. Dotted lines represent disordered linkers connecting the ECDs. (C) Detailed view centered on the pocket occupied by the C-terminal AM Y52 and CGRPmut F37. (D) Model of AM S45W/Q50W bound to RAMP2-CLR ECD. Two conformations are modeled for Q50W. (E and F) Competition AlphaLISA peptide-binding assays with purified MBP-RAMP1-(GS)5-CLR ECD (E) or MBP-RAMP2[L160R]-(GS)5-CLR ECD fusion proteins (F) and indicated competitor peptides. Each panel is representative of two to three independent experiments performed in duplicate. Individual data points for each technical replicate are shown. (G) Amino acid sequence alignment of human AM(13–52), αCGRP(1–37), and CGRPmut. Numbers above the sequences correspond to AM amino acid position, whereas numbers below the sequences correspond to CGRP amino acid position. Dark-blue line represents the N-terminal disulfide linkage.
In our crystal structures of peptide-bound tethered MBP-RAMP1/2 ECD-CLR ECD fusion proteins, we used C-terminal antagonist fragments of a high-affinity CGRP N31D/S34P/K35F variant, which we termed CGRPmut for simplicity, and AM (Fig. 1, B and C) (Booe et al., 2015). The three CGRPmut substitutions were originally identified through systematic mutagenesis (Rist et al., 1998). The structures revealed that CGRPmut S34P promotes a key β-turn structure, and K35F contacts the CLR “turret loop” to enhance affinity (Fig. 1C). Given the shared peptide β-turn structures, we reasoned that Phe substitution at AM Q50, which is equivalent to CGRP K35 (Fig. 1G), should enhance its affinity. Using in silico mutagenesis in PyMOL, we modeled Q50F and substitutions at several positions in AM. We predicted that a Trp at position 50 could adopt either of two conformations to enhance affinity by contacting the turret loop, and that an additional Trp at position 45 in the helical turn would permit stacking of the Trp residues to stabilize the β-turn (Fig. 1D). In the context of AM(37–52), we measured the effects of the Q50F, Q50W, and S45W/Q50W substitutions on binding to the purified tethered fusion proteins using an AlphaLISA luminescent proximity peptide-binding assay similar to those previously described (Moad and Pioszak, 2013; Lee et al., 2016). In a competition binding format, the Q50F substitution enhanced affinity ∼10-fold for both tethered RAMP1-CLR ECD and RAMP2-CLR ECD (data not shown), whereas Q50W enhanced affinity ∼100-fold for both proteins (Fig. 1, E and F). The inclusion of S45W with Q50W further enhanced affinity such that the double mutant bound both tethered ECD complexes with IC50 values in the mid nanomolar range (Fig. 1, E and F).
The competition AlphaLISA assay is ill suited for determining accurate equilibrium constants, so we developed a fluorescence polarization/anisotropy peptide-binding assay similar to one we recently reported for the calcitonin receptor ECD (Lee et al., 2017) to enable a more thorough quantitative analysis of the binding of AM and CGRP variants to the purified tethered RAMP1/2-CLR ECD complexes (Supplemental Fig. S1A). As a probe, we used AM(37–52) S45W/Q50W labeled at the N terminus with FITC. The label had no effect on ECD-binding affinity (Supplemental Fig. S1, B and C). A robust, saturable anisotropy signal was observed upon titration of a fixed probe concentration with MBP-RAMP1-CLR ECD, but not with an unrelated MBP-RSPO fusion protein, and nonspecific binding was negligible (Supplemental Fig. S1, D and E). The FITC-AM(37–52) S45W/Q50W probe exhibited a pKD of 7.72 ± 0.12 for MBP-RAMP1-CLR ECD and a pKD of 6.83 ± 0.04 for MBP-RAMP2-CLR ECD (Fig. 2, A and B). In a competition binding assay format, unlabeled AM(37–52) S45W/Q50W exhibited a pKI of 7.40 ± 0.05 for MBP-RAMP1-CLR ECD and a pKI of 7.19 ± 0.06 for MBP-RAMP2-CLR ECD (Fig. 2, C and D; Table 1). In the competition assay, the affinity of AM(37–52) for both the tethered CGRP and AM1 receptor ECD complexes was increased by the rationally designed substitutions with a rank order of affinity-enhancement of S45W/Q50W > Q50W > Q50F (Table 1). Notably, the double mutant dramatically increased affinity from the mid- to high-micromolar range into the mid-nanomolar range.
FP peptide-binding assay utilizing an AM(37–52) S45W/Q50W probe. (A and B) Saturation binding assays with 7 nM FITC-AM(37–52) S45W/Q50W and indicated concentrations of purified MBP-RAMP1-(GS)5-CLR ECD (A) or MBP-RAMP2[L106R]-(GS)5-CLR ECD (B) tethered fusion proteins. (C and D) Competition binding assays using 7 nM FITC-AM(37–52) S45W/Q50W and 15 nM MBP-RAMP1-(GS)5-CLR ECD (C) or 110 nM MBP-RAMP2[L106R]-(GS)5-CLR ECD complexes (D) and indicated concentration of unlabeled competitor AM(37–52) S45W/Q50W. Each panel is representative of at least three independent experiments performed in duplicate. Individual data points for each technical replicate are shown.
Modeling predicted that a Trp at CGRP position 35 could adopt either of two conformations to contact the turret loop and possibly confer higher affinity than the K35F substitution in CGRPmut (Fig. 3A). We further predicted that Ser at position 36 could adopt either of two conformations in which its side chain hydrogen bonds with either the main chain carbonyl of D31 or the main chain G33 carbonyl and the C-terminal amide to stabilize the β-turn (Fig. 3B). In the FP assay, CGRPmut bound RAMP1-CLR ECD with a pKI of 7.48 ± 0.12 (KI ∼ 33 nM) in good agreement with the ∼30 nM IC50 value reported for its binding to the intact CGRP receptor in SK-N-MC cell membranes (Rist et al., 1998) (Fig. 3C; Supplemental Fig. S2B; Table 1). Replacing K35F in CGRPmut with K35W increased affinity for RAMP1-CLR ECD by ∼4.5-fold (pKI = 8.12 ± 0.11), and adding A36S into CGRPmut increased affinity ∼6-fold (pKI = 8.24 ± 0.04) (Supplemental Fig. S2B; Supplemental Table S2; Table 1). Combining these two substitutions yielded the CGRP(27–37) N31D/S34P/K35W/A36S variant with ∼28-fold improved affinity for RAMP1-CLR ECD (pKI = 8.93 ± 0.14) over that of CGRPmut and an approximately 1700-fold weaker affinity for RAMP2-CLR ECD (pKI = 5.64 ± 0.13) (Fig. 3, C and D; Supplemental Table S2; Table 1).
Enhancing CGRP peptide affinity through rational structure-based design. (A) Modeled CGRP K35W showing two possible conformations. (B) Modeled CGRP A36S showing hydrogen bonds in two possible conformations. (C and D) Representative competition binding FP assays using 7 nM FITC-AM(37–52) S45W/Q50W and 15 nM MBP-RAMP1-(GS)5-CLR ECD (C) or 110 nM MBP-RAMP2[L106R]-(GS)5-CLR ECD (D) with competitor peptides CGRPmut or CGRP(27–37) N31D/S34P/K35W/A36S. Each panel is representative of at least three independent experiments performed in duplicate. Individual data points for each technical replicate are shown.
Probing the Role of RAMP Allosteric Effects and RAMP-Peptide Contacts in Ligand Selectivity through AM and CGRP Variants.
AM Y52 and CGRP F37 are critical for receptor binding (Moad and Pioszak, 2013; Watkins et al., 2013), but their role in selectivity is less clear. CGRPmut F37 contacts RAMP1 W84 and AM Y52 hydrogen bonds with RAMP2 E101, whereas RAMP2 F111 cannot make the CGRP F37 contact and RAMP1 W74 cannot hydrogen bond with AM Y52, suggesting that distinct RAMP-peptide contacts contribute to selectivity (Fig. 4A). CGRPmut F37 is the only residue that contacts RAMP1, but AM binding is more complex because, in addition to Y52, K46 also contacts RAMP2 (Fig. 4B). Despite the apparent hydrogen bond network involving AM K46 and Y52 and RAMP2 R97, E101, and E105, only mutation of RAMP2 E101 impaired AM1 function in a signaling assay; the RAMP2 R97A and E105A mutants retained wild-type response to AM (Booe et al., 2015). The primary role of AM K46 appears to be to pack against Y52 and the CLR Trp shelf (Fig. 4B), and thus, only the AM C-terminal residue makes a contact of importance to the RAMP subunit. Given the paucity of direct RAMP-peptide contacts, we reasoned that RAMP alteration of CLR conformation might also contribute to selectivity.
Designing peptide substitutions to exploit RAMP-peptide contacts and putative RAMP allosteric modulation of CLR ECD conformation. (A) Detailed view of AM and CGRPmut C-terminal residue contacts with the RAMP subunits (Protein Data Bank 4RWF and 4RWG). Dotted line represents a hydrogen bond. (B) Hydrogen bond network involving AM K46, Y52, and RAMP2 R97, E101, and E105 (Protein Data Bank 4RWF). Hydrogen bond distances in angstroms are in black text. (C) Differences in the CLR β1-β2 loop position in the RAMP1- and RAMP2-CLR ECD complexes (Protein Data Bank 4RWF and 4RWG). The red arrow indicates the putative allosteric pathway propagating changes from the RAMP:CLR interface to the CLR β1-β2 loop (see also Supplemental Movie S1). (D) Modeling AM K46L suggests that it would push Y52 to a position similar to that of CGRPmut F37, allowing contact with G71 in RAMP1-CLR while clashing with G71 in RAMP2-CLR. Semitransparent space-filling spheres are shown for the modeled Leu at AM position 46, CGRPmut F37, and G71 in RAMP1-CLR ECD.
RAMP-mediated conformational changes in the CLR ECD were visualized by comparing the CGRPmut-bound RAMP1-CLR ECD and AM-bound RAMP2-CLR ECD structures in a molecular morph (Supplemental Movie S1). The morph highlights three areas of CLR change that appear to be RAMP-dependent. The first involves different conformations of CLR R119 on the turret loop adjacent to the RAMP α2-α3 loop. The second involves a subtle shift in the β1-β2 loop that appears to be a change propagated from the RAMP:CLR interface to the loop (Fig. 4C; Supplemental Movie S1). The third involves the N terminus of CLR α1, which is “unwound” in the RAMP2-CLR structure. We reasoned the CLR R119 change is not a significant selectivity determinant because the apparently RAMP1-induced “down” position of R119 would seem to sterically disfavor AM binding, but our aforementioned results indicated that AM variants can bind RAMP1-CLR ECD with high affinity. We also disregarded the change in the CLR α1 N terminus as it may be an artifact of tethering it to the RAMP ECD C terminus. We therefore focused on the subtle shift in the CLR β1-β2 loop, with the idea that it would alter the shape of the pocket occupied by the peptide C-terminal residue. Notably, this β1-β2 loop shift was also evident in ligand-free RAMP1:CLR and RAMP2:CLR ECD crystal structures (ter Haar et al., 2010; Kusano et al., 2012) (Supplemental Fig. S3).
Fortunately, the unique helical turn in AM appeared to provide an avenue to probe the putative allosteric effect. Using in silico mutagenesis, we predicted that a Leu at AM position 46 would push the AM Y52 side chain to adopt a position similar to that of CGRPmut F37 to favor RAMP1-CLR ECD binding due to better complementarity with the pocket while sterically disfavoring RAMP2-CLR ECD binding, because the altered Y52 position would clash with G71 in the β1-β2 loop (Fig. 4D). AM(37–52) peptides bearing K46 substitution with Leu or other hydrophobic amino acids were tested for binding to MBP-RAMP1- and RAMP2-CLR ECD proteins in a single-point competition FP assay (Fig. 5, A and B). We also tested K46R to see if introduction of an Arg would enable ionic interaction with RAMP2 E101 and/or E105. The AM K46L and K46M variants appeared to have enhanced affinity for RAMP1-CLR ECD (Fig. 5A), whereas all of the K46 variants had decreased binding to RAMP2-CLR ECD (Fig. 5B). Combining K46L with the affinity-enhancing Q50W or S45W/Q50W yielded variants that were 26- to 29-fold selective for RAMP1-CLR ECD over RAMP2-CLR ECD, with the triple mutant exhibiting affinity for RAMP1-CLR ECD in the low-nanomolar range (pKI = 8.14 ± 0.10) (Supplemental Fig. S2A; Table 1). The FP peptide-binding data for all peptide variants at purified ECD complexes are summarized in pKI correlation plots (Fig. 6), and representative competition binding curves are shown in Supplemental Fig. S2. Statistical analysis is presented in table (Supplemental Table S2; Table 1) and scatter plot (Supplemental Fig. S4) formats.
Binding of AM variants substituted at position 46 and the AMmut quadruple mutant to purified MBP-RAMP1/2-CLR ECD complexes. (A and B) Single-point competition FP assays using 7 nM FITC-AM S45W/Q50W and 15 nM MBP-RAMP1-(GS)5-CLR ECD (A) or 110 nM MBP-RAMP2[L106R]-(GS)5-CLR ECD (B) with indicated concentration of competitor peptides. Each of the variants is in the AM(37–52) scaffold. (C and D) Competition binding FP assays using 7 nM FITC-AM(37–52) S45W/Q50W and 15 nM MBP-RAMP1-(GS)5-CLR ECD (C) or 110 nM MBP-RAMP2[L106R]-(GS)5-CLR ECD (D) with competitor peptides AM(37–52) and the S45W/K46L/Q50W/Y52F variant (AMmut). (A) and (B) are representative of at least two independent experiments performed in duplicate with error bars shown as S.D., and (C) and (D) are representative of three independent experiments performed in duplicate with individual data points for each technical replicate shown.
Summary of AM and CGRP pKI values obtained in the FP assay with purified ECD complexes. Plot of the mean pKI values for AM variants (A) or CGRP variants (B) at either MBP-RAMP1-CLR ECD or MBP-RAMP2-CLR ECD purified fusion proteins. AM variants are in the context of AM(37–52), whereas CGRP variants are in the context of CGRP(27–37). Open circles indicate that pKI values were only obtainable at one of the purified RAMP1/2-CLR ECD complexes, whereas binding the other receptor had an estimated pKI of <4.3. Error bars indicate 95% confidence intervals. See Supplemental Fig. S4, Supplemental Table S2, and Table 1 for statistical analysis.
Next, we examined C-terminal residue swaps in the high-affinity AM and CGRP variants to probe the role of direct RAMP-peptide contacts in selectivity. Inclusion of Y52F with K46L/Q50W, S45W/K46L/Q50W, or S45W/K46M/Q50W in AM dramatically decreased affinity to RAMP2-CLR ECD and thereby enhanced selectivity for RAMP1-CLR ECD (Fig. 6A; Supplemental Figs. S2A and S4A; Table 1). The AM(37–52) S45W/K46L/Q50W/Y52F quadruple mutant, hereafter designated AMmut, exhibited strong nanomolar affinity for RAMP1-CLR ECD (pKI = 8.47 ± 0.21) and >2000-fold weaker micromolar affinity for RAMP2-CLR ECD (pKI = 5.13 ± 0.03) (Fig. 5, C and D; Table 1). Addition of F37Y into CGRP N31D/S34P/K35W/A36S had no significant effect on its affinity for RAMP1-CLR ECD, but enhanced affinity ∼5-fold for RAMP2-CLR ECD (pKI = 6.33 ± 0.11) (Fig. 6B; Supplemental Fig. S2B and S4B; Table 1).
Antagonism of AM and CGRP Variants at Full-Length Receptors in Cells.
The abilities of selected AM and CGRP variants to antagonize cAMP signaling at the three RAMP:CLR complexes transiently expressed in COS-7 cells was assessed to extend our results to the full-length receptors in a natural cell membrane environment. Importantly, this cell-based assay also allowed us to determine the potency of the variants at the AM2 receptor, for which we lack the corresponding purified RAMP3-CLR ECD complex due to difficulties with its expression and purification. Antagonist apparent pKB values were determined from the dextral displacement of the agonist concentration-response curves resulting from a fixed concentration of the variants. Representative data for AMmut and CGRP(27–37) N31D/S34P/K35W/A36S antagonism at each receptor are shown in Fig. 7, A and B, and representative data for the other AM and CGRP variants are shown in Supplemental Figs. S5 and S6, respectively. The apparent pKB values and the receptor selectivity of each AM and CGRP variant are listed in Table 2, and these data are summarized in apparent pKB heat maps (Fig. 8). Statistical analysis is shown in table (Table 2; Supplemental Table S3) and scatter plot (Supplemental Fig. S7) formats.
Antagonism of AM and CGRP variants at intact RAMP:CLR complexes in COS-7 cells. Concentration-response curves for CGRP, AM1, and AM2 receptors transiently expressed in COS-7 cells stimulated with αCGRP (CGRP receptor) or AM (AM receptors) agonists in the presence or absence of the indicated concentrations of AMmut (A) or CGRP(27–37) N31D/S34P/K35W/A36S (B). V, represents vehicle control with no agonist. Each panel is representative of three independent experiments performed in duplicate. Individual data points for each technical replicate are shown.
Heat map summary of AM and CGRP variant apparent pKB values obtained in the cell-based signaling assay in COS-7 cells. Apparent pKB values obtained for transiently expressed CGRP receptor (RAMP1), AM1 receptor (RAMP2), or AM2 receptor (RAMP3) stimulated with αCGRP (CGRP receptor) or AM (AM receptors) in the presence or absence of the indicated AM antagonist variants (A) or CGRP antagonist variants (B). AM variants are in the context of AM(37–52), whereas CGRP variants are in the context of CGRP(27–37). The mean apparent pKB values for AM(37–52) and AM Q50F at the CGRP receptor were estimated to be less than 5.3 (indicated by X symbol). See Table 2 for apparent pKB values with associated error and Supplemental Fig. S7 and Supplemental Table S3 for statistical analysis.
The Q50F and Q50W substitutions enhanced the antagonist potency of AM(37–52); however, these variants displayed selectivity for the AM receptors over the CGRP receptor (∼15-fold for Q50W) that was not observed in the ECD-binding assay (Fig. 8A; Table 2). Combining K46L with Q50W increased potency for the CGRP receptor and decreased potency at the AM1 receptor such that the K46L/Q50W double mutant was relatively nonselective with high-nanomolar range potency. Including the C-terminal Y52F swap in K46L/Q50W/Y52F appeared to further decrease potency at the AM1 receptor and increase potency at the CGRP receptor. Although these potency changes due to Y52F did not reach statistical significance (Supplemental Fig. S7A; Supplemental Table S3), the trends were consistent with a C-terminal Phe favoring RAMP1 and disfavoring RAMP2. The S45W/Q50W double mutant had increased potency at all three receptors, but again, in contrast to the ECD-binding data, this variant was selective for the AM receptors over the CGRP receptor. Adding K46L to S45W/Q50W increased potency for the CGRP receptor such that the S45W/K46L/Q50W triple mutant was nonselective with mid-nanomolar range potency. Last, addition of Y52F to yield AMmut significantly reduced potency at the AM1 receptor. AMmut thus had a dramatically altered receptor-binding profile as compared with the AM(37–52) parent such that it exhibited strong nanomolar potency at the CGRP (pKBapp = 7.60 ± 0.08) and AM2 receptors (pKBapp = 7.12 ± 0.21) with 56-fold selectivity for CGRP over AM1 and ∼19-fold selectivity for AM2 over AM1 (Fig. 8A; Supplemental Fig. S7A; Table 2).
Replacing K35F in CGRPmut with K35W increased potency at all three receptors ∼3- to 8-fold while not substantially altering selectivity (Fig. 8B; Supplemental Fig. S7B; Table 2). Inclusion of A36S increased potency for each receptor an additional ∼4- to 6-fold such that the quadruple N31D/S34P/K35W/A36S variant had strong low-nanomolar range potency for the CGRP receptor (pKBapp = 8.62 ± 0.05) with a selectivity profile similar to that of the CGRPmut parent (54-fold selectivity for CGRP over AM1 and ∼9-fold selectivity for AM2 over AM1). Notably, CGRP(27–37) N31D/S34P/K35W/A36S was comparable in activity to the much longer traditional CGRP(8–37) antagonist. Last, adding the F37Y C-terminal residue swap had no effect on potency at the CGRP receptor, but increased potency at the AM1 receptor. The N31D/S34P/K35W/A36S/F37Y variant thus exhibited strong nanomolar potency for all three receptors with a reduced ability to discriminate CGRP and AM1, although it still preferred the CGRP receptor over AM1 (Fig. 8B; Supplemental Fig. S7B; Table 2).
Structural Basis of Enhanced Affinity and Altered Selectivity of AMmut for the CGRP Receptor.
To see if the AM substitutions functioned as designed, we crystallized MBP-RAMP1-CLR ECD with bound AMmut (Supplemental Fig. S8A) and determined the structure at a resolution of 2.80 Å (Supplemental Table S4). Electron density was observed for the peptide in each of the four molecules present in the asymmetric unit, but the density for the ECD complex and peptide in Mol 4 (protein chain D:peptide chain H) was particularly poor, so we ignored this complex in our analysis. The best density and lowest B-factors were observed for Mol 2 (protein chain B:peptide chain F) and Mol 3 (protein chain C:peptide chain G) (Supplemental Fig. S8B; Supplemental Table S4). We largely focused on Mol 2 because, although AMmut in Mol 3 was the most well ordered, the MBP in this molecule is shifted closer to the peptide as compared with the other molecules, such that it packs against AMmut S45W, and we therefore cannot rule out an effect of MBP on the conformation of this residue (Supplemental Fig. S8, C and D).
AMmut occupies RAMP1-CLR ECD with the same mainchain conformation observed for AM bound to RAMP2-CLR ECD (Fig. 9, A and B). The Trp at position 50 adopts the rotamer in which the benzene ring of the indole contacts AM P49, thus presumably stabilizing the peptide β-turn in addition to making contact with the CLR turret loop as designed. The conformation of the Trp at position 45 was not unambiguously defined except in Mol 3, where MBP may have constrained its position. In Mol 3 (peptide chain G), Trp45 stacked against Trp50 (Supplemental Fig. S8D). In Mol 1 and Mol 2, the density for the Trp45 side chain suggested alternative conformations, and we modeled two such conformations in Mol 2 (peptide chain F) that involve the indole contacting either Trp50 or the introduced Leu at position 46 (Fig. 9A; Supplemental Fig. S8D). Leu46 packed against the phenyl ring of the Phe introduced at position 52 and pushed it to a position similar to that of CGRPmut F37 as designed (Fig. 9, C and D). RAMP1-CLR ECD with bound AMmut exhibited the same β1-β2 loop position observed with bound CGRPmut. CLR R119 adopted the “down” conformation presumably induced by RAMP1, but its position was slightly shifted to accommodate the difference in the AM β-turn position as compared with that of CGRPmut (Fig. 9C). Last, the CLR Trp shelf and β3-β4 loop adopted conformations similar to those observed in AM-bound RAMP2-CLR ECD, presumably because these changes are peptide-dependent rather than RAMP-dependent (Fig. 9, B and C).
Structural basis for AMmut binding to the RAMP1-CLR ECD complex. (A) Cartoon representation of the 2.8-Å resolution crystal structure of AMmut-bound MBP-RAMP1-CLR ECD with AMmut in gray, CLR ECD in blue, and RAMP1 ECD in red. Two alternate conformations were modeled for S45W. The (GSA)3 linker sequence connecting the two ECDs was disordered. MBP is omitted for clarity. (B) Superimposition of AMmut-bound RAMP1-CLR ECD [colored as in (A)] and AM (gold) bound to RAMP2 ECD (orange)-CLR ECD (cyan) (Protein Data Bank 4RWF). (C) Superimposition of AMmut-bound RAMP1-CLR ECD [colored as in (A)] and CGRPmut (violet) bound to RAMP1 ECD (yellow)-CLR ECD (green) (Protein Data Bank 4RWG). (D) Detailed view comparing peptide-binding interactions for the three peptide-bound RAMP ECD-CLR ECD heterodimer crystal structures (current structure and Protein Data Bank 4RWF and 4RWG). Key peptide and RAMP residues are labeled and colored as in (A)–(C).
Discussion
AM and CGRP signaling is complex because they activate three similar heterodimeric RAMP:CLR complexes. Crystal structures of the CGRPmut- and AM-bound tethered RAMP1-CLR and RAMP2-CLR ECD complexes, respectively, revealed that the peptides share a common binding site on CLR and make only a single contact of importance with the RAMP subunits (Booe et al., 2015). We posited that the RAMP subunits determine peptide selectivity through a combination of their distinct contacts with the AM/CGRP C-terminal residue and their modulation of CLR conformation, i.e., allosteric effects. Here, we explored this hypothesis through the use of rationally designed AM and CGRP antagonist variants. Our purpose was 3-fold: 1) develop short high-affinity AM and CGRP variants, 2) probe the contributions of RAMP allosteric effects and direct peptide contacts to ligand selectivity, and 3) inform efforts to develop novel potent and selective AM/CGRP-based pharmacological tools and therapeutics.
Through structure-guided design targeting AM positions 45 and 50 and CGRP positions 35 and 36, we increased the affinity of the minimal C-terminal antagonist fragment of AM from the micromolar into the nanomolar range and further increased the affinity of CGRPmut (Rist et al., 1998) approximately 20-fold. Our results are consistent with the idea that stabilizing the AM/CGRP β-turn in the receptor-bound conformation increases affinity. Indeed, exploiting turn structures is an established strategy in engineering peptide ligands for G protein–coupled receptors (Gruber et al., 2010). The AM Q50W substitution enhanced affinity more than Q50F, which presumably contacts the CLR turret loop similar to CGRPmut F35. The crystal structure of AMmut-bound RAMP1-CLR ECD revealed that the Trp50 indole contacts the CLR turret loop, but it also packs against AM P49 (Fig. 9). We predict that this intramolecular interaction stabilizes the β-turn and is largely responsible for the potent effect of this substitution. Addition of a Trp at position 45 presumably further stabilizes the receptor-bound conformation by stacking against Trp50 and/or Leu46 in AMmut. Trp at position 35 in CGRP enhanced affinity more than K35F in CGRPmut, and we suspect that the Trp indole packs against Pro at CGRP position 34. The CGRP A36S substitution also appears to enhance affinity by stabilizing the β-turn, probably via hydrogen bonds with the G33 carbonyl and the C-terminal amide (Fig. 3B). It is unclear if these putative β-turn–stabilizing effects occur in the free peptides or only when receptor-bound. Interestingly, evolution has yielded “substitutions” in other species similar to some of those found here at AM position 50 and CGRP position 36 (Supplemental Fig. S9, A and B). Zebrafish AM has a Tyr at position 50, which may enhance affinity by contacting the CLR turret loop, similar to a Phe at this position. A frog CGRP has Ser at position 36, and its (27–37) fragment exhibited nanomolar affinity (KI = 95 nM) for the human CGRP receptor (Ladram et al., 2008).
The high-affinity AM and CGRP variants provided ideal backgrounds in which to add substitutions designed to probe the mechanism of RAMP-altered CLR ligand selectivity because their enhanced affinity made it easier to accurately measure pKI and apparent pKB values. We focused on substitutions at AM position 46 to probe RAMP-dependent changes in the CLR β1-β2 loop position, i.e., allostery, and on swapping the C-terminal peptide residues to probe direct RAMP-peptide contacts. In the AM Q50W background, K46L increased affinity at the CGRP receptor and decreased affinity at the AM1 receptor (Figs. 6A and 8A), and the crystal structure indicated that K46L shifted the position of the peptide C-terminal residue phenyl ring as designed (Fig. 9D). We cannot rule out the possibility that the K46L effect is due to the loss of a hydrogen bond between the K46 side chain and RAMP2 E101 (Fig. 4B), but we think this is not the case because K46L not only decreased binding to RAMP2:CLR, but also increased binding to RAMP1:CLR. In the S45W/Q50W background, K46L only increased affinity for RAMP1:CLR (Figs. 6A and 8A). These data are consistent with K46L causing AM to be more complementary to the shape of the RAMP1:CLR pocket occupied by the C-terminal residue. In AMmut, the C-terminal residue swap Y52F further decreased antagonism at the AM1 receptor (Fig. 8A) because loss of the Y52 side chain hydroxyl removes the ability to hydrogen bond with RAMP2 E101, and RAMP2 F111 cannot contact the remaining phenyl ring (Fig. 4A). The CGRP conformation did not lend itself to probing the RAMP allosteric effects as we did with AM, but we probed the RAMP-CGRP contacts, and these data also provide insight into allostery. Giving CGRP N31D/S34P/K35W/A36S the ability to hydrogen bond with RAMP2 E101 via the F37Y substitution enhanced its ability to bind AM1, but it clearly still preferred the CGRP receptor (Figs. 6B and 8B). This result is probably due to RAMP1/2-dependent conformational differences in CLR because only the CGRP C-terminal residue contacts the RAMP. Altogether, our data are consistent with the hypothesis that RAMPs determine CLR ligand selectivity through a dual mechanism comprising distinct contacts with the peptide C-terminal residue and subtle allosteric effects on CLR conformation. There are also data suggesting an allosteric component to RAMP function at the related calcitonin receptor (Gingell et al., 2016; Lee et al., 2016).
It is informative to compare the ECD complex-binding data (Fig. 6) and the cell-based antagonism data (Fig. 8). For the CGRP variants, there was good agreement between the pKI and pKBapp values for RAMP1:CLR, suggesting that the binding determinants for these peptides reside mostly within the ECD. In contrast, the CGRP variants bound intact RAMP2:CLR in cells better than the tethered RAMP2-CLR ECD, and as a consequence, they were less selective in the cell-based assay than in the FP assay (Tables 1 and 2). For the AM variants, the relationship between the ECD binding and cell-based antagonism results was more complex. Some variants did not differentiate RAMP1-CLR ECD and RAMP2-CLR ECD, but preferred intact RAMP2:CLR over RAMP1:CLR (e.g., Q50W and S45W/Q50W). Consequently, whereas adding K46L into these variants caused them to prefer RAMP1-CLR ECD over RAMP2-CLR ECD, they became nonselective for intact RAMP1:CLR and RAMP2:CLR. AM(37–52)–based variants may be at a disadvantage at full-length RAMP1:CLR because their N-terminal region occupies the ECD complex differently than CGRP(27–37) (Fig. 9C). Similar to the CGRP variants, the AM variants were less selective in the cell-based assay than in the FP assay. The divergence in the degree of receptor selectivity between the ECD-binding and signaling data may be due to complexity in the intact receptors. The orientation of the ECD complex with respect to the remainder of the receptor within the full-length complexes is unknown, and portions of the membrane-embedded domains may also contact the C-terminal “ECD-binding” half of the peptides (Zhang et al., 2017). In addition, the ECD complexes were expressed in E. coli so they are not N-glycosylated, which might affect peptide binding as for the related calcitonin receptor (Lee et al., 2017; Liang et al., 2017). Despite some disagreement in the degree of receptor selectivity, the relative activities of the AM/CGRP variants were reasonably consistent across the two assays.
The cell-based data provide insight into the molecular determinants of ligand binding at the AM2 receptor, for which we lack the purified RAMP3-CLR ECD complex and a crystal structure. The two AM receptors are very similar, but there is evidence for distinct AM1/AM2 agonist-binding site properties (Watkins et al., 2016). Our data are consistent with these findings even for short antagonist fragments that are expected to primarily bind the ECD complex. Notably, AMmut discriminated the AM receptors with ∼19-fold selectivity, which is better than existing antagonists, although still not ideal for use as a tool. The AMmut C-terminal Phe probably did not decrease binding to AM2 because RAMP3 appears to be a RAMP1/2 hybrid in terms of how it augments the CLR pocket. RAMP3 W84 and E74 enable contact with peptides having either a C-terminal Phe or Tyr (Supplemental Fig. S10). The AM1/AM2-discriminating ability of AMmut appears to result in part from the ability of its C-terminal Phe to contact RAMP3 W84, but not RAMP2 F111. The conserved RAMP1/3 residue W84 suggests that peptide design strategies that seek to discriminate the RAMP1 and RAMP3 complexes may benefit from exploiting the putative allosteric effects of RAMPs rather than RAMP-peptide contacts.
This study highlights the challenge of developing peptides that strongly discriminate the CGRP, AM1, and AM2 receptors when starting from peptides that primarily bind CLR. Nonetheless, we developed short AM and CGRP variants with significantly enhanced affinities, and we were able to swap the selectivity of AM for the CGRP and AM1 receptors with the AMmut variant. The dramatically altered receptor binding preferences of AMmut resulted from substitutions designed to exploit the differences in how RAMP1 and RAMP2 both augment the CLR pocket and alter its shape. These results support a dual mechanism for RAMP function in which direct contact with the peptide C-terminal residue and allosteric modulation of CLR cooperate to determine ligand selectivity. The novel peptide variants developed herein and the new crystal structure should help guide future efforts to develop potent and selective AM/CGRP-based pharmacological tools and therapeutics.
Acknowledgments
The authors thank Joe Brunzelle for assistance with remote data collection at beamline 21-ID-G (LS-CAT) of the Advanced Photon Source.
Authorship Contributions
Participated in research design: Booe, Hay, Pioszak.
Conducted experiments: Booe, Warner, Roehrkasse.
Performed data analysis: Booe, Warner, Hay, Pioszak.
Wrote or contributed to the writing of the manuscript: Booe, Hay, Pioszak.
Footnotes
- Received October 20, 2017.
- Accepted January 19, 2018.
This work was supported by grants from the National Institutes of Health [R01GM104251] and the Oklahoma Center for the Advancement of Science and Technology [HR16-005] (A.A.P.). The microscope used for crystal imaging was supported by an Institutional Development Award (IDeA) from the National Institute of General Medical Sciences of the NIH [P20GM103640]. The Advanced Photon Source is a US Department of Energy Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357. Use of LS-CAT sector 21 was supported by the Michigan Economic Development Corporation and the Michigan Technology Tri-Corridor [Grant 085P1000817]. D.L.H. is supported by a James Cook Research Fellowship from the Royal Society of New Zealand.
The coordinates and structure factors for the AMmut-bound MBP-RAMP1 ECD-CLR ECD fusion protein were deposited with the worldwide protein data bank with ID code 5V6Y.
Peptide variants described in this manuscript are the subject of a patent application (AAP).
↵
This article has supplemental material available at molpharm.aspetjournals.org.
Abbreviations
- Ahx
- aminohexanoic acid
- AM
- adrenomedullin
- AMmut
- AM(37-52)NH2 S45W/K46L/Q50W/Y52F
- BSA
- bovine serum albumin
- CGRP
- calcitonin gene–related peptide
- CGRPmut
- CGRP(27-37)NH2 N31D/S34P/K35F
- CLR
- calcitonin receptor–like receptor
- DMEM
- Dulbecco’s modified Eagle’s medium
- ECD
- extracellular domain
- FITC
- fluorescein isothiocyanate
- FP
- fluorescence polarization
- MBP
- maltose-binding protein
- MR
- molecular replacement
- RAMP
- receptor activity–modifying protein.
- Copyright © 2018 by The American Society for Pharmacology and Experimental Therapeutics