|
|
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Molecular Pharmacology and Experimental Therapeutics, Mayo Clinic Scottsdale, Scottsdale, Arizona (M.D., F.G., K.H., D.I.P., L.J.M.); Department of Molecular Biology, Scripps Research Institute and Molsoft LLC, La Jolla, California (P.C.-H.L., R.A.); and Department of Pharmacology, Monash University, Clayton, Victoria, Australia (P.M.S.)
Received February 22, 2007; accepted May 2, 2007
| Abstract |
|---|
|
|
|---|
Insights into mechanisms of agonist ligand binding and activation of the secretin receptor have come from natural ligand structure-activity relationships and from receptor mutagenesis and photoaffinity labeling studies (Gourlet et al., 1996
; Dong et al., 2003
). These have identified the critical importance of the long and structurally complex amino terminus of this receptor for agonist ligand binding. This theme has been consistent for multiple other class II GPCR family members (Juppner et al., 1994
; Cao et al., 1995
; Gourlet et al., 1996
). Alignment of the receptor amino terminal sequences reveals conserved features, which include six cysteine residues with three intradomain disulfide bonds (Grauschopf et al., 2000
; Lisenbee et al., 2005
). Insight into the structure of this region was substantially advanced with the solution of an NMR structure of the isolated amino terminus of the corticotropin-releasing factor (CRF2β) receptor (Grace et al., 2004
). However, it remains unclear whether the amino terminus of the intact CRF receptor will have the same or related conformation, how similar the amino terminal region of other receptors in this family will be to that structure, or how ligands dock to these receptors.
Natural ligands for class II GPCRs are moderately long peptides (greater than 25 residues) having diffuse pharmacophoric domains (Ulrich et al., 1998
). Like natural ligands for other members of this family, the amino-terminal region of secretin contains key determinants for receptor selectivity and activation, whereas its carboxyl-terminal region contains determinants for high-affinity binding. Considerable experimental evidence indicates that the carboxyl-terminal regions of ligands for class II receptors interact with the amino-terminal regions of their receptors (Juppner et al., 1994
; Gourlet et al., 1996
; Pham et al., 2004
; Tan et al., 2006
), whereas the amino terminus of the peptides is critical for agonist activity and interacts with parts of the receptor core and loop regions (Bisello et al., 1998
; Dong et al., 2004
). Critical to understanding receptor activation is knowledge of how these ligands dock to their receptors. The most useful insights have come from intrinsic photoaffinity labeling studies that provide physical constraints on the relative position of peptide and receptor amino acids. The class II GPCR that has been most extensively studied this way is the secretin receptor, where spatial approximations have been established between residues scattered throughout the pharmacophore of secretin, in positions 6, 12, 13, 14, 18, 22, and 26, and residues within the receptor amino terminus (Dong et al., 1999a
, b
, 2000
, 2002
, 2003
; Zang et al., 2003
).
To date, the only photoaffinity labeling probes shown to covalently interact with receptor residues in regions other than the amino terminus are those having photolabile residues substituted into the amino-terminal end of the peptide ligand (Dong et al., 2004
). This type of observation, also reported for the parathyroid hormone receptor, has resulted in the hypothesis that natural peptide ligands of class II GPCRs can act as tethers interacting with both amino terminus and body of their receptors, exerting tension to induce the conformational change associated with activation and G protein association (Bisello et al., 1998
). We have proposed an alternative hypothesis (Dong et al., 2006
) that describes a ligand-induced conformational change in the receptor amino terminus that exposes a previously "constrained" epitope (WDN in the secretin receptor) that can act as an endogenous agonist that interacts with and activates the body of the receptor.
Evaluation of these two apparently divergent activation mechanisms requires a more detailed understanding of the mechanism of binding and of the conformation of the complex between secretin and the amino terminus of the secretin receptor. We have therefore performed additional photoaffinity labeling experiments with new probes that incorporate photolabile residues into the important carboxyl-terminal region, in positions 21 or 23 of secretin. Both probes bound specifically and saturably to the secretin receptor and were full agonists. They were also able to covalently label single and distinct residues within the amino terminus of the receptor, at residues Arg15 (for position 21 probe) and Arg21 (for position 23 probe). We also have now used the NMR data for the amino-terminal domain of the structurally related CRF2β receptor (Grace et al., 2004
) to refine that structure and to build a homology model of the analogous domain of the secretin receptor. We have oriented this relative to a homology model of the bovine rhodopsin helical bundle (Palczewski et al., 2000
) to yield a preliminary full secretin receptor model. The previously identified spatial approximation constraints from photoaffinity labeling studies were used to provide initial ligand docking. This model was found to also nicely accommodate the two new spatial approximation constraints identified in the current work, with all of the currently available constraints contributing toward the most meaningful model of ligand docking yet proposed.
| Materials and Methods |
|---|
|
|
|---|
Synthetic Peptides. Photolabile secretin probes used in this study were [Tyr10,Bpa21]rat secretin (Bpa21 probe) and [Tyr10,Bpa23]rat secretin (Bpa23 probe). They incorporated a photolabile p-benzoyl-L-phenylalanine (Bpa) to replace Arg21 and Leu23 located at the carboxyl-terminal half of the ligand, respectively. Like the well characterized radioligand [Tyr10]rat secretin, they incorporated a Tyr to replace Leu10, which has been shown to be well tolerated. Peptides were synthesized as described previously (Powers et al., 1988
), radioiodinated oxidatively using the solid-phase oxidant Iodobeads, and purified to homogeneity using reversed-phase high-performance liquid chromatography to yield specific radioactivities of 2000 Ci/mmol.
Receptor-Expressing Cell Lines. Several cell lines that had been established and well characterized previously were used as sources of receptor for this study. These included Chinese hamster ovary (CHO) cell lines stably expressing the wild-type secretin receptor (SecR) (Ulrich et al., 1993
), the HA-tagged wild-type secretin receptor (SecR-HA37) (Dong et al., 1999b
), and HA-tagged mutant secretin receptors (SecR-V16M-HA37 and SecR-V13M-HA37) (Dong et al., 2000
). Numbering of the residues in the secretin receptor reflects the cleavage of the 22-residue signal peptide sequence (Dong et al., 1999a
). These cell lines were cultured at 37°C in a 5% CO2 environment on Falcon tissue culture plasticware in Ham's F-12 medium supplemented with 5% Fetal Clone-2 (HyClone Laboratories, Logan, UT). Cells were passaged twice a week and lifted mechanically before use. Plasma membranes were prepared using discontinuous sucrose gradient centrifugation (Hadac et al., 1996
). Membranes were then resuspended in Krebs-Ringers/HEPES (KRH) buffer containing 25 mM HEPES, pH 7.4, 104 mM NaCl, 5 mM KCl, 2 mM CaCl2, 1 mM KH2PO4, 1.2 mM MgSO4, 0.01% soybean trypsin inhibitor, and 1 mM phenylmethylsulfonyl fluoride for storage at -80°C until they were to be used.
Binding Studies. Binding activities of the Bpa21 and Bpa23 probes were performed with enriched plasma membranes from the CHO-SecR cells in standard competition-binding assays, using conditions previously established (Hadac et al., 1996
). Membranes (containing
10 µg of protein) were incubated with a constant amount of radioligand [125I-[Tyr10]rat secretin (5-10 pM)] and increasing concentrations of nonradioactive ligand (0-1 µM Bpa21 or Bpa23 probe) in KRH buffer containing 0.2% bovine serum albumin for 1 h at room temperature (final volume, 500 µl). Separation of bound from free radioligand was performed with a Skatron cell harvester (Molecular Devices, Sunnyvale, CA), using receptor-binding filter mats that had been presoaked in 0.3% Polybrene. Bound radioactivity was quantified using a gamma-spectrometer. Nonspecific binding was determined in the presence of 1 µM unlabeled secretin and represented less than 15% of total binding. Data were graphed using Prism software (GraphPad Software, San Diego, CA) and were analyzed using the nonlinear least-squares curve-fitting routine in LIGAND (Munson and Rodbard, 1980
).
Biological Activity Assay. Biological activities of the Bpa21 and Bpa23 probes were studied using assays for stimulation of intracellular cAMP in the receptor-bearing CHO-SecR cell line, as described previously (Ganguli et al., 1998
). CHO-SecR cells were incubated with increasing concentrations of the probe (0-1 µM) for 30 min at 37°C in KRH buffer containing 0.2% bovine serum albumin, 0.01% soybean trypsin inhibitor, 0.1% bacitracin, and 1 mM 3-isobutyl-1-methylxanthine. Reactions were stopped by adding 6% (w/v) perchloric acid, and the pH was adjusted to 6.0 with 30% KHCO3. The cAMP concentrations in cell lysates were measured using a competition-binding assay (Diagnostic Products Corporation, Los Angeles, CA) (Ganguli et al., 1998
). Radioactivity was quantified by scintillation counting in a Beckman LS6000 (Beckman Coulter, Fullerton, CA).
Photoaffinity Labeling. The detailed procedure has been described previously (Dong et al., 1999b
). Plasma membranes (
50 µg) were incubated with the 125I-labeled Bpa21 or Bpa23 probes (
0.1 nM) in KRH buffer for 1 h in the dark at room temperature in the presence of increasing concentrations of unlabeled secretin. The reaction was then exposed to photolysis for 30 min at 4°C using a Rayonet Photochemical Reactor (Southern New England Ultraviolet Co., Bradford, CT) equipped with 3500 Å lamps. Membranes were then washed with KRH buffer, solubilized in SDS sample buffer, and separated by gel electrophoresis on 10% SDS-polyacrylamide gels (Laemmli, 1970
). The apparent molecular masses of radiolabeled receptors, visualized by autoradiography, were determined by interpolation on a plot of the mobility of ProSieve protein standards (Lonza Rockland, Inc., Rockland, ME) versus the log values of their apparent masses.
Peptide Mapping. Photoaffinity-labeled secretin receptor was prepared in larger scale using
200-µg aliquots of plasma membranes and 0.5 nM 125I-labeled Bpa21 or Bpa23 probes. Labeled bands were cut from gels, eluted, lyophilized, and precipitated with ethanol before subsequent cleavage. Purified labeled receptor was deglycosylated with endoglycosidase F and digested with CNBr and endoproteinase Lys-C under conditions described previously (Dong et al., 1999b
). Cleavage products were separated on 10% NuPAGE gels (Invitrogen, Carlsbad, CA) using MES running buffer. The apparent molecular masses of radiolabeled receptor fragments were determined by interpolation on a plot of the mobility of Multimark protein standards (Invitrogen) versus the log values of their apparent masses.
The identity of the affinity-labeled secretin receptor fragment resulting from CNBr cleavage was further established by immunoprecipitation with the 12CA5 anti-HA monoclonal antibody using the affinity-labeled HA-tagged receptor construct (SecR-HA37), previously shown to bind and signal identically to wild-type receptor (Dong et al., 1999b
). The immunoprecipitated material was eluted into sample buffer and resolved by NuPAGE gel electrophoresis.
Radiochemical Edman degradation sequencing (Dong et al., 1999b
) was used to determine the specific sites of attachment after achieving definitive identification of the receptor fragments labeled with each probe. This involved the cross-linking of the receptor fragment of interest through cysteine residues to N-(2-aminoethyl-1)-3-aminopropyl glass beads. We used V13M-HA37 and V16M-HA37 mutant receptor constructs expressed in CHO cell lines for the Bpa21 and Bpa23 probes, respectively. Labeled receptor was purified and cleaved by CNBr and the resultant fragments were gel-purified to radioactive homogeneity before being covalently coupled through cysteine residues to maleimidobenzoyl succinimide-activated N-(2-aminoethyl-1)-3-aminopropyl glass beads. After this, repetitive cycles of manual Edman degradation were performed with quantitation of radioactivity released in each cycle. This procedure was repeated three times in independent experiments.
Molecular Modeling. All molecular modeling activities were conducted using a stochastic global energy optimization procedure implemented in Internal Coordinate Mechanics (ICM) (Abagyan et al., 1994
). This procedure consisted of the following iterative steps: 1) random conformational change of a dihedral angle according to the biased-probability Monte Carlo method (Abagyan and Totrov, 1994
); 2) local minimization of all free dihedral angles; and 3) acceptance or rejection of the new conformation based on the Metropolis criterion at the simulation temperature, usually at 600 K (Metropolis et al., 1953
). This procedure can generate and search through diverse sets of conformations by actively sampling a selected set of dihedral angles. All calculations were carried out on computers equipped with 3.4-GHz Intel XEON-EMT processors.
Refinement of the NMR Models of the Mouse CRF2β Receptor Amino-Terminal Domain. The 20 NMR models of the mouse CRF2β receptor amino-terminal domain and the associated NMR distance restraints that were deposited into Protein Data Bank (PDB code 1U34) were obtained. The highly flexible regions at both ends of this domain were removed, retaining only residues 37 to 122. The models were then completed by attaching missing hydrogens, assigning ECEPP/3 atom types and charges (Nemethy et al., 1992
), followed by regularization to impose ideal covalent geometries on all bond angles and bond lengths. Eight hundred and three NMR restraints were converted to ICM-readable format and were imposed on the models using a quadratic restraint energy with a strength of 10 kcal/(mol Å2) (Bordner and Abagyan, 2006
). Each of the 20 models was then subjected to local minimization, followed by 10 independent cycles of Monte Carlo sampling of both the backbone and side-chain dihedral angles at 600 K in ICM, each simulation lasting for more than 144 h (Abagyan and Totrov, 1999
). The lowest energy conformation for each of the 20 models was retained.
Homology Modeling of the Amino-Terminal Domain of the Rat Secretin Receptor. The multiple sequence alignment of class II GPCRs from the Swissprot database was used as a guide to align the sequence of the rat secretin receptor to the mouse CRF2β receptor amino-terminal domain (Boeckmann et al., 2003
). The 20 refined models of this domain of the CRF2β receptor were then used to generate 20 homology models of the amino terminus of the secretin receptor (Cardozo et al., 1995
). A subset of 115 NMR restraints of the CRF2β receptor amino terminus was imposed on the residues that are conserved between the secretin receptor and the CRF2β receptor. Each of the 20 models was then subjected to local minimization, followed by 10 independent cycles of Monte Carlo sampling of backbone and side-chain conformations at 600 K, each simulation lasting for more than 216 h. The lowest energy conformation for each of the 20 models was retained and further refined by Monte Carlo simulation at 100 K for 36 h.
Homology Modeling of the Transmembrane Domain of the Rat Secretin Receptor. The sequence of the transmembrane helices of the rat secretin receptor were aligned to those of bovine rhodopsin (PDB code 1U18) using the "cold spot" method (Frimurer and Bywater, 1999
). A homology model of the transmembrane helical bundle of the rat secretin receptor was then built in ICM. The packing of the model was further optimized by side-chain sampling and backbone minimization.
Alignment of the Amino-Terminal Domain of the Secretin Receptor with Its Transmembrane Helical Bundle. A distance restraint was imposed between the WDN epitope and the top of the sixth transmembrane segment. A second distance restraint was imposed between the carboxyl end of the amino-terminal domain and the amino-terminal end of helix 1. The amino-terminal domain was then docked onto the helical bundle by Monte Carlo sampling of the six positional variables of the amino-terminal domain. The contact between the two domains was further optimized by side chain sampling to generate a preliminary model of the secretin receptor incorporating both the amino terminus and receptor core (Fernández-Recio et al., 2003
).
|
| Results |
|---|
|
|
|---|
The probes were then tested for their ability to stimulate intracellular cAMP accumulation in CHO-SecR cells. As shown in Fig. 1, both probes were fully efficacious agonists, stimulating full intracellular cAMP responses that were not different from those achieved in response to natural secretin. They were, however, less potent than natural secretin (EC50 values: secretin, 0.06 ± 0.01 nM; Bpa21 probe, 0.37 ± 0.02 nM; Bpa23 probe, 5.22 ± 0.97 nM).
Photoaffinity Labeling of the Secretin Receptor. We explored the ability of the Bpa21 and Bpa23 probes to covalently label the secretin receptor. As shown in Fig. 2, both probes specifically and saturably labeled the secretin receptor, with the labeled protein bands migrating at approximate Mr = 70,000 and shifting to approximate Mr = 42,000 after deglycosylation with endoglycosidase F. This labeling was competed by secretin in a concentration-dependent manner (Bpa21 probe, IC50 = 18.1 ± 2.8 nM; Bpa23 probe, IC50 = 3.8 ± 0.5 nM). No radioactive band was observed in the affinity-labeled nonreceptor-bearing CHO cell membranes.
|
|
|
|
To further localize the region of labeling by each of the probes, two previously established secretin receptor mutant constructs, V13M-HA37 and V16M-HA37 (Dong et al., 2000
), were used. Both receptor mutants were specifically and saturably labeled with each probe (Fig. 6). CNBr cleavage of the V13M-HA37 mutant labeled with the Bpa21 probe yielded a fragment migrating on a 10% NuPAGE gel at approximate Mr = 19,000 that shifted to Mr = 9000 after deglycosylation, representing the fragment between Arg14 and Met51. CNBr cleavage of the V16M-HA37 mutant labeled with the Bpa21 probe yielded a nonglycosylated fragment migrating at approximate Mr = 5000, representing the fragment between Ala1 and Met16 (Fig. 6). Taken together, these data indicate that the site of labeling with the Bpa21 probe was within the small segment between Arg14 and Val16.
|
Radiochemical Edman degradation sequencing of the purified labeled fragments resulting from the CNBr cleavage of the labeled V13M-HA37 (for the Bpa21 probe) and V16M-HA37 (for the Bpa23 probe) receptor mutants was performed. Shown in Fig. 7 are the profiles of eluted radioactivity in which a peak was found in cycle 2 for the Bpa21 probe and in cycle 5 for the Bpa23 probe. This identified Arg15 and Arg21 of the secretin receptor as the sites of labeling with the Bpa21 probe and for the Bpa23 probe, respectively.
|
1/
2 rotamer normality, and backbone conformation after simulation.
|
Nederveen et al. (2005
) published the RECOORD database, in which more than 500 NMR models from the Protein Data Bank were recalculated. The largest improvement in the model quality was achieved by restrained molecular dynamics in a hydrated environment. Comparing their study with ours shows that refinement in ICM improved the model quality similar to that of Nederveen's study in all categories in PROCHECK and WHAT_CHECK, bringing the overall quality of the amino terminus of the CRF2β receptor to that of an average NMR model. Note that although the Z-score of the backbone conformation after refinement was still -8.4, it represented an improvement over the Z-score of the original models of -11.4. Consistent with this, Nederveen's refinement also did not improve the backbone conformations significantly.
To see how the refinement affected the overall structure of the models, we calculated the pairwise C
root-mean-square deviation (rmsd) by superimposing the C
of the stable central core residues 59 to 64, 70 to 81, and 100 to 113. For each of the 20 models, we superimposed the model before and after refinement; the resultant C
rmsd for the central core and the remaining disordered region were 1.22 ± 0.16 Å and 7.8 ± 1.6 Å, respectively. These values reflected the degree of movement in the central core and the disordered region during refinement. As a reference, the pairwise C
rmsd for the central core and the disordered region within the ensemble of the original 20 NMR models were 0.84 ± 0.17 Å and 7.7 ± 1.9 Å, respectively, which indicate the intrinsic flexibility/uncertainty in the two regions in the original NMR models. These results demonstrate that our simulation enables large movement in the disordered region while limiting the movement of the stable central core within experimental uncertainty, as shown in Fig. 8.
|
The ensemble of 20 rat secretin receptor amino terminus models after homology modeling and refinement is shown in Fig. 9. We calculated the pairwise C
rmsd within the ensemble by superimposing the C
of the stable central core residues 43 to 48, 53 to 64, and 82 to 93. The C
rmsd for the central core and the remaining disordered region were 2.11 ± 0.60 Å and 11.3 ± 2.0 Å, respectively. A comparison of these numbers with those of the original NMR models of the amino terminus of the CRF2β receptor (0.84 ± 0.17 and 7.7 ± 1.9 Å) shows a wider spread of the supposed stable central core, probably reflecting the smaller number of distance restraints used in the modeling of this region of the secretin receptor.
|
Only seven of the 20 models were found to contain an intact salt bridge between Asp49 and Arg83, corresponding to the salt bridge between Asp65 and Arg101 proposed to be present in the amino terminus of the mouse CRF2β receptor. We inspected the distance restraints deposited in Protein Data Bank and found that there was no direct distance restraint connecting Asp65 and Arg101 in the original NMR models, even though this salt bridge was reported as a prominent feature of those conformations.
Molecular Model of the Intact Secretin Receptor. To study how secretin interacts with its receptor, we first constructed an initial model of the receptor including the receptor amino terminus and the seven-transmembrane segment bundle. A preliminary model of the transmembrane helical bundle domain was generated using a previously reported "cold-spot" method (Frimurer and Bywater, 1999
). We attached the seven secretin receptor amino-terminal domain models having an intact Asp49 to Arg83 salt bridge (that were most similar to the reported conformation of the CRF2β receptor) to this receptor transmembrane domain. The relative orientation of the two domains was restrained by the connection of the two domains at the amino-terminal end of helix 1 and by the proposed interaction between the WDN sequence and the transmembrane domain (Dong et al., 2006
). Seven models of the receptor that were consistent with these restraints were generated.
Secretin was first docked to the grid representation of the seven full models to maximize sampling efficiency. Because of the rigid nature of the receptor in grid docking, we subsequently carried out flexible backbone refinement. The C
distances between cross-linked residues in secretin and the secretin receptor from photoaffinity labeling studies are shown for three of the best models in Table 2. Note that the two photoaffinity labeling groups, p-benzoyl-L-phenylalanine and p-benzoylbenzoyl-L-lysine, can span distances of 10 and 15 Å, respectively. Therefore, the seven distance restraints between secretin and the amino-terminal domain of the secretin receptor coming from previous photoaffinity labeling studies were all satisfied.
|
Most importantly, the last two entries of Table 2 show that the two pairs of cross-linked residues being demonstrated for the first time in the current study (Arg21 of secretin to receptor residue Arg15 and Leu23 of secretin to receptor residue Arg21) were found to be within cross-linking distances in our models, even though these two sets of distance restraints had not been used during the peptide docking and refinement. Figure 10 shows the best working model for secretin docking to the receptor amino terminus and highlights the potential accessibility of the WDN sequence within the amino terminus of the receptor previously postulated as playing an endogenous agonist role by interacting with the body of the secretin receptor (Dong et al., 2006
).
|
| Discussion |
|---|
|
|
|---|
Useful structural insights have come from the solution of the NMR structure of the amino-terminal domain of another class II GPCR, the CRF2β receptor (Grace et al., 2004
). In the current work, we have carefully examined those reported structures and the distance constraints that had contributed to the reported models. In the current report, the CRF2β receptor model was refined, and the insights coming from it were applied to the analogous region of the structurally related secretin receptor. This probably represents a more meaningful template for a homology model for this region of the secretin receptor than the amino-terminal region of ribonuclease Mc1, which was previously used as template in the absence of the CRF2β receptor structure (Dong et al., 2003
). Indeed, the core regions of the amino-terminal structures of the secretin receptor and of the CRF2β receptor have now been found to be very similar.
Most of the experimentally defined sites of residue-residue approximation that have been determined for the secretin receptor are within the first 38 residues of the distal amino terminus in poorly constrained regions of this structure. Of particular interest is that the regions of the secretin receptor that are homologous to the relatively stable core of the amino terminus of the CRF2β receptor have not been covalently labeled by any of positions of photolabile residues spread throughout the pharmacophoric domain of secretin receptor probes that have been used to date. This raises the possibility that this stable core provides a lattice deep in the structure with the portions of the receptor that interact with the natural peptide ligands situated above it.
Although placement of the peptide ligands within the amino-terminal domain is crucial to understanding high-affinity binding for class II receptors, an additional key factor to resolving the mechanism of receptor activation by agonists is the positioning of both the peptide ligand and receptor amino-terminal domain relative to the heptahelical transmembrane core.
There are currently few empirical data to direct the organization of the receptor amino-terminal domain and the transmembrane receptor core, and this has led to disparate models of the relative positioning of these two domains. For the CRF2β receptor, the orientation used was based on a proposed interaction between a negatively charged surface of the receptor amino-terminal domain and positively charged residues present in the extracellular loops of the receptor core (Grace et al., 2004
), whereas the VPAC1 receptor amino-terminal domain was arbitrarily placed above the receptor core in the recently proposed model of the intact VPAC1 receptor (Tan et al., 2006
). In both cases, the docking of peptides to their receptors relied significantly on general consideration of the functional roles of the amino- and carboxyl-terminal ends of the peptide ligands across class II receptors that are dependent upon the orientation of the receptor amino-terminal domain and receptor core. These models are further complicated by lack of direct knowledge of regions of the amino terminus that are disordered in the NMR structure that form the basis for receptor amino-terminal domain modeling.
In the current report, we have used two restraints in determining how to best position the amino-terminal domain of the secretin receptor with its helical bundle domain. The first restraint is the integrity of the peptide backbone involved in connecting the two domains, bringing the amino-terminal end of the first transmembrane helix to the carboxyl-terminal end of the receptor amino terminus. This is a common feature of each of the receptor models; however, the innate flexibility of the distal end of the receptor amino-terminal domain allows multiple surfaces to potentially abut the receptor core. The second restraint is specific to the secretin receptor and assumes that the WDN sequence within the amino terminus of the receptor interacts with the helical bundle above transmembrane segment six (Dong et al., 2006
). This was based on a series of studies demonstrating the ability of this receptor sequence to act as an endogenous agonist ligand, with its site of interaction demonstrated by covalent labeling of the receptor body using two distinct photoaffinity labeling probes derived from the critical WDN tripeptide (Dong et al., 2006
). It is noteworthy that similar activation of VPAC1 and calcitonin receptors by homologous endogenous peptide sequences has also been shown to occur, suggesting that such activation may be a common mechanism within the class II GPCR family (Dong et al., 2006
).
Secretin was docked to this preliminary model of the secretin receptor using the seven published photoaffinity labeling spatial approximation constraints. Of particular interest is that, with the peptide docked this way and given the proposed orientation of the two domains of the secretin receptor, the amino terminus of the peptide was found to be adjacent to the helical bundle in the region it has been independently shown to photoaffinity label (Dong et al., 2004
). In addition, the two spatial approximation constraints newly identified in the current report were also compatible with this docking and orientation. It is encouraging that these independently derived observations are fully compatible with the working model.
Critical comparison of the currently proposed model of the mode of peptide docking to the amino-terminal domain of the secretin receptor with those of the astressin/CRF-occupied CRF2β receptor (Grace et al., 2004
) and of the vasoactive intestinal polypeptide-occupied VPAC1 receptor (Tan et al., 2006
), shows clear differences. As discussed above, the orientations of the amino-terminal domains with the helical bundle domains of these receptors are all different. The secretin receptor orientation was guided by the WDN motif and was, therefore, fully consistent with its interaction with the relevant region of the receptor body. However, in the CRF2β receptor model, the amino-terminal domain interacts with the transmembrane domain through the "back side" of the model displayed in Figs. 8 and 9. As a consequence, Leu64 to Gln66 of the CRF2β receptor, corresponding to the WDN motif, is facing upward and away from the analogous region of that receptor. In the VPAC1 receptor model, the amino-terminal domain makes contact with the transmembrane domain near the end of the second β-sheet containing residues Ile61 to Cys63 and Pro77 to Gly79. In this orientation, the region analogous to the WDN motif of the secretin receptor is facing sideward. However, there was no clear rationalization for this orientation with respect to the helical bundle.
Comparison of the peptide-binding regions of these receptors also reflects differences. In the current model of the secretin receptor, the peptide is situated between the distal amino terminus, the region containing the majority of the sites of cross-linking, and part of the stable core. Glu9 to Ala17 of secretin interacts with the tip of the first β-sheet, including residues Asp49 to Ser52, and the disordered loop comprising residues Lys77 to Ser80. Gln20 to Val27 of secretin interacts with the disordered loop, including residues Glu31 to Gly37 (Fig. 10). In the proposed CRF2β receptor structure, regions proposed to represent the peptide-binding site based on chemical-shift perturbations upon antagonist (astressin) binding include Ile67 to Thr69, Gly90 to Asn93, Glu102 to Cys103, and Arg112 to Ser116. The peptide-binding region in both the secretin and CRF2β receptor models includes the tip of the first β-sheet and the disordered loop near the palm of the second β-sheet. In the proposed VPAC1 receptor structure, the peptide-binding interface was inferred from three pairs of cross-linked residues, at positions 6, 22, and 24 of the peptide ligand. Here, this region includes receptor residues Gln80 to Leu84, Val101 to Cys105, and Trp110 to Cys122. It is noteworthy that this interface would be at the "bottom side" of the model of the secretin receptor displayed in Figs. 8 and 9, opposite the proposed peptide-binding sites in both the currently proposed model of the secretin receptor and that of the CRF2β receptor. At present, we do not have sufficient empirical data to confirm whether the alternate models proposed for different receptors reflect divergence in mechanism of action or potential problems in generation of working models.
Another structural motif of substantial interest is the proposed intradomain salt bridge between Asp65 and Arg101 of the CRF2β receptor. Of particular interest, the region of the secretin receptor amino terminus that has been shown to have endogenous agonist activity is centered on the residue analogous to CRF2β receptor residue Asp65. If such a bond were present and stable in the secretin receptor (it would involve the secretin receptor residues Asp49 and Arg83), it might make the endogenous agonist hypothesis less likely as a mechanism of activation of these receptors. In the secretin receptor amino terminus models that were built based on homology with the CRF2β receptor structure, the salt bridge between Asp49 and Arg83 actually did not exist in the majority of conformations. Whether the absence of this salt bridge in the secretin receptor models is an artifact caused by an inadequate number of distance restraints in modeling or a genuine feature of the secretin receptor amino terminus remains to be explored. It should be noted, however, that there is no explicit salt-bridge present from the CRF2β NMR distance restraints, whereas mutation of Asp49 to glycine has no negative effect on the biological activity of the secretin receptor (Dong et al., 2006
). It is also recognized that solvent-accessible salt bridges are easily broken. Such a breakable salt-bridge between Asp49 and Arg83 of the secretin receptor is still consistent with the endogenous agonist hypothesis (Dong et al., 2006
). However, in the best conformation of the secretin-occupied secretin receptor currently being proposed, the WDN motif is readily accessible to its proposed target of action, even with the salt-bridge incorporated.
It is encouraging that, in the current work, multiple experimentally derived residue-residue distance constraints were easily accommodated within the molecular model of the secretin receptor. With initial peptide docking accomplished based on the previously described seven constraints, it was even more encouraging that the currently derived two additional constraints were also then well accommodated. In the best working model, all existing experimentally derived constraints were fully accommodated, and the WDN motif was available for interaction with the body of the receptor.
| Acknowledgements |
|---|
| Footnotes |
|---|
M.D. and P.C.-H.L. are coprimary authors.
ABBREVIATIONS: GPCR, G protein-coupled receptor; CRF, corticotropin-releasing factor; CHO, Chinese hamster ovary; CRF2β, corticotropin-releasing factor receptor; HA, hemagglutinin; Bpa, p-benzoyl-L-phenylalanine; KRH, Krebs-Ringers/HEPES; MES, 2-(N-morpholino)ethanesulfonic acid; PDB, Protein Data Bank; ICM, internal coordinate mechanics; rmsd, root-mean-square deviation; SecR, secretin receptor; VPAC1, vasoactive intestinal polypeptide receptor.
Address correspondence to: Laurence J. Miller, M.D., Mayo Clinic, 13400 East Shea Blvd., Scottsdale, AZ 85259. E-mail: miller{at}mayo.edu
| References |
|---|
|
|
|---|
Abagyan R, Totrov M, and Kuznetsov D (1994) ICM—a new method for protein modeling and design: applications to docking and structure prediction from the distorted native confirmation. J Comput Chem 15: 488-506.[CrossRef]
Abagyan RA and Totrov M (1999) Ab initio folding of peptides by the optimal-bias Monte Carlo minimization procedure. J Comput Phys 151: 402-421.[CrossRef]
Bisello A, Adams AE, Mierke DF, Pellegrini M, Rosenblatt M, Suva LJ, and Chorev M (1998) Parathyroid hormone-receptor interactions identified directly by photocross-linking and molecular modeling studies. J Biol Chem 273: 22498-22505.
Boeckmann B, Bairoch A, Apweiler R, Blatter MC, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I, et al. (2003) The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res 31: 365-370.
Bordner AJ and Abagyan R (2006) Ab initio prediction of peptide-MHC binding geometry for diverse class I MHC allotypes. Proteins 63: 512-526.[CrossRef][Medline]
Cao YJ, Gimpl G, and Fahrenholz F (1995) The amino-terminal fragment of the adenylate cyclase activating polypeptide (PACAP) receptor functions as a high affinity PACAP binding domain. Biochem Biophys Res Commun 212: 673-680.[CrossRef][Medline]
Cardozo T, Totrov M, and Abagyan R (1995) Homology Modeling by the Icm Method. Proteins 23: 403-414.[CrossRef][Medline]
Clore GM, Nilges M, Brunger A, and Gronenborn AM (1988) Determination of the backbone conformation of secretin by restrained molecular dynamics on the basis of interproton distance data. Eur J Biochem 171: 479-484.[Medline]
Dong M, Asmann YW, Zang M, Pinon DI, and Miller LJ (2000) Identification of two pairs of spatially approximated residues within the carboxyl terminus of secretin and its receptor. J Biol Chem 275: 26032-26039.
Dong M, Li Z, Pinon DI, Lybrand TP, and Miller LJ (2004) Spatial approximation between the amino terminus of a peptide agonist and the top of the sixth transmembrane segment of the secretin receptor. J Biol Chem 279: 2894-2903.
Dong M, Li Z, Zang M, Pinon DI, Lybrand TP, and Miller LJ (2003) Spatial approximation between two residues in the mid-region of secretin and the amino terminus of its receptor. Incorporation of seven sets of such constraints into a three-dimensional model of the agonist-bound secretin receptor. J Biol Chem 278: 48300-48312.
Dong M, Pinon DI, Asmann YW, and Miller LJ (2006) Possible endogenous agonist mechanism for the activation of secretin family G protein-coupled receptors. Mol Pharmacol 70: 206-213.
Dong M, Wang Y, Hadac EM, Pinon DI, Holicky E, and Miller LJ (1999a) Identification of an interaction between residue 6 of the natural peptide ligand and a distinct residue within the amino-terminal tail of the secretin receptor. J Biol Chem 274: 19161-19167.
Dong M, Wang Y, Pinon DI, Hadac EM, and Miller LJ (1999b) Demonstration of a direct interaction between residue 22 in the carboxyl-terminal half of secretin and the amino-terminal tail of the secretin receptor using photoaffinity labeling. J Biol Chem 274: 903-909.
Dong M, Zang M, Pinon DI, Li Z, Lybrand TP, and Miller LJ (2002) Interaction among four residues distributed through the secretin pharmacophore and a focused region of the secretin receptor amino terminus. Mol Endocrinol 16: 2490-2501.
Fernández-Recio J, Totrov M, and Abagyan R (2003) ICM-DISCO docking by global energy optimization with fully flexible side-chains. Proteins 52: 113-117.[CrossRef][Medline]
Frimurer TM and Bywater RP (1999) Structure of the integral membrane domain of the GLP1 receptor. Proteins 35: 375-386.[CrossRef][Medline]
Ganguli SC, Park CG, Holtmann MH, Hadac EM, Kenakin TP, and Miller LJ (1998) Protean effects of a natural peptide agonist of the G protein-coupled secretin receptor demonstrated by receptor mutagenesis. J Pharmacol Exp Ther 286: 593-598.
Gourlet P, Vilardaga JP, De Neef P, Waelbroeck M, Vandermeers A, and Robberecht P (1996) The C-terminus ends of secretin and VIP interact with the N-terminal domains of their receptors. Peptides 17: 825-829.[CrossRef][Medline]
Grace CR, Perrin MH, DiGruccio MR, Miller CL, Rivier JE, Vale WW, and Riek R (2004) NMR structure and peptide hormone binding site of the first extracellular domain of a type B1 G protein-coupled receptor. Proc Natl Acad Sci U S A 101: 12836-12841.
Grauschopf U, Lilie H, Honold K, Wozny M, Reusch D, Esswein A, Schafer W, Rucknagel KP, and Rudolph R (2000) The N-terminal fragment of human parathyroid hormone receptor 1 constitutes a hormone binding domain and reveals a distinct disulfide pattern. Biochemistry 39: 8878-8887.[CrossRef][Medline]
Hadac EM, Ghanekar DV, Holicky EL, Pinon DI, Dougherty RW, and Miller LJ (1996) Relationship between native and recombinant cholecystokinin receptors: role of differential glycosylation. Pancreas 13: 130-139.[Medline]
Hooft RWW, Vriend G, Sander C, and Abola EE (1996) Errors in protein structures. Nature 381: 272.[Medline]
Juppner H, Schipani E, Bringhurst FR, McClure I, Keutmann HT, Potts JT Jr, Kronenberg HM, Abou-Samra AB, Segre GV, and Gardella TJ (1994) The extracellular amino-terminal region of the parathyroid hormone (PTH)/PTH-related peptide receptor determines the binding affinity for carboxyl-terminal fragments of PTH-(1-34). Endocrinology 134: 879-884.
Laemmli UK (1970) Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 227: 680-685.[CrossRef][Medline]
Laskowski RA, Macarthur MW, Moss DS, and Thornton JM (1993) Procheck—a program to check the stereochemical quality of protein structures. J Appl Crystallogr 26: 283-291.[CrossRef]
Lisenbee CS, Dong M, and Miller LJ (2005) Paired cysteine mutagenesis to establish the pattern of disulfide bonds in the functional intact secretin receptor. J Biol Chem 280: 12330-12338.
Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH, and Teller E (1953) Equation of state calculations by fast computing machines. J Chem Phys 21: 1087-1092.[CrossRef]
Munson PJ and Rodbard D (1980) Ligand: a versatile computerized approach for characterization of ligand-binding systems. Anal Biochem 107: 220-239.[CrossRef][Medline]
Nederveen AJ, Doreleijers JF, Vranken W, Miller Z, Spronk C, Nabuurs SB, Guntert P, Livny M, Markley JL, Nilges M, et al. (2005) RECOORD: A recalculated coordinate database of 500+ proteins from the PDB using restraints from the BioMagResBank. Proteins 59: 662-672.[CrossRef][Medline]
Nemethy G, Gibson KD, Palmer KA, Yoon CN, Paterlini G, Zagari A, Rumsey S, and Scheraga HA (1992) Energy parameters in polypeptides. 10. Improved geometrical parameters and nonbonded interactions for use in the ECEPP/3 algorithm, with application to proline-containing peptides. J Phys Chem 96: 6472-6484.[CrossRef]
Palczewski K, Kumasaka T, Hori T, Behnke CA, Motoshima H, Fox BA, Le Trong I, Teller DC, Okada T, Stenkamp RE, et al. (2000) Crystal structure of rhodopsin: a G protein-coupled receptor. Science 289: 739-745.
Pearson RK, Miller LJ, Hadac EM, and Powers SP (1987) Analysis of the carbohydrate composition of the pancreatic plasmalemmal glycoprotein affinity labeled by short probes for the cholecystokinin receptor. J Biol Chem 262: 13850-13856.
Pham V, Wade JD, Purdue BW, and Sexton PM (2004) Spatial proximity between a photolabile residue in position 19 of salmon calcitonin and the amino terminus of the human calcitonin receptor. J Biol Chem 279: 6720-6729.
Powers SP, Pinon DI, and Miller LJ (1988) Use of N,O-bis-Fmoc-D-Tyr-ONSu for introduction of an oxidative iodination site into cholecystokinin family peptides. Int J Pept Protein Res 31: 429-434.[Medline]
Tan YV, Couvineau A, Murail S, Ceraudo E, Neumann JM, Lacapere JJ, and Laburthe M (2006) Peptide agonist docking in the N-terminal ectodomain of a class II G protein-coupled receptor, the VPAC1 receptor. Photoaffinity, NMR, and molecular modeling. J Biol Chem 281: 12792-12798.
Ulrich CD 2nd, Holtmann M, and Miller LJ (1998) Secretin and vasoactive intestinal peptide receptors: members of a unique family of G protein-coupled receptors. Gastroenterology 114: 382-397.[CrossRef][Medline]
Ulrich CD 2nd, Pinon DI, Hadac EM, Holicky EL, Chang-Miller A, Gates LK, and Miller LJ (1993) Intrinsic photoaffinity labeling of native and recombinant rat pancreatic secretin receptors. Gastroenterology 105: 1534-1543.[Medline]
Zang M, Dong M, Pinon DI, Ding XQ, Hadac EM, Li Z, Lybrand TP, and Miller LJ (2003) Spatial approximation between a photolabile residue in position 13 of secretin and the amino terminus of the secretin receptor. Mol Pharmacol 63: 993-1001.
This article has been cited by other articles:
![]() |
M. Dong, P. C.-H. Lam, D. I. Pinon, P. M. Sexton, R. Abagyan, and L. J. Miller Spatial Approximation between Secretin Residue Five and the Third Extracellular Loop of Its Receptor Provides New Insight into the Molecular Basis of Natural Agonist Binding Mol. Pharmacol., August 1, 2008; 74(2): 413 - 422. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Dong, F. Gao, D. I. Pinon, and L. J. Miller Insights into the Structural Basis of Endogenous Agonist Activation of Family B G Protein-Coupled Receptors Mol. Endocrinol., June 1, 2008; 22(6): 1489 - 1499. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. G. Harikumar, P. C.-H. Lam, M. Dong, P. M. Sexton, R. Abagyan, and L. J. Miller Fluorescence Resonance Energy Transfer Analysis of Secretin Docking to Its Receptor: MAPPING DISTANCES BETWEEN RESIDUES DISTRIBUTED THROUGHOUT THE LIGAND PHARMACOPHORE AND DISTINCT RECEPTOR RESIDUES J. Biol. Chem., November 9, 2007; 282(45): 32834 - 32843. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||