![]()
|
|
||||||||
J. Biol. Chem., Vol. 278, Issue 41, 40010-40019, October 10, 2003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

¶

||
From the
Laboratoire de Pharmacochimie de la Communication Cellulaire, UMR 7081 CNRS, 74 Route du Rhin, 67401 Illkirch, France and ¶INSERM U469, 141 Rue de la Cardonille, 34094 Montpellier Cedex 5, France
Received for publication, February 3, 2003 , and in revised form, July 16, 2003.
| ABSTRACT |
|---|
|
|
|---|
-chloroacetamide moieties. The F225C, F308C, and K128C mutants of the V1a receptor were expressed in COS-7 or Chinese hamster ovary cells, and their pharmacological properties toward SR49059 and its sulfydryl-reactive analogues were analyzed. We demonstrated that treatment of the F225C mutant with the isothiocyanate-derivative compound led to dose-dependent inhibition of the residual binding of the radio-labeled antagonist [125I]HO-LVA. This inhibition is probably the consequence of a covalent irreversible chemical modification, which is only possible when close contacts and optimal orientations exist between reactive groups created both on the ligand and the receptor. This result validated the three-dimensional model hypothesis. Thus, we propose that residue Phe225, located in transmembrane domain V, directly participates in the binding of the V1a-selective nonpeptide antagonist SR49059. This conclusion is in complete agreement with all our previous data on the definition of the agonist/antagonist binding to members of the oxytocin/vasopressin receptor family. | INTRODUCTION |
|---|
|
|
|---|
Up to now, identification of the binding site of the SR49059 into the V1aR has not been investigated at a molecular level. Several studies based on a combination of receptor three-dimensional modeling and site-directed mutagenesis experiments have suggested that an AVP-binding pocket is buried into a 1520-Å-deep central cavity of the V1a receptor, defined by the transmembrane helices and surrounded by the extracellular loops (17, 18). Extracellular residues also play a role in the binding of the hormone (19, 20). Binding domains for synthetic peptide antagonists overlap those for peptide agonists into the AVP/OT receptors; however, discrimination of agonist versus antagonist ligands is achieved by conserved aromatic residues (Trp304, Phe307, and Phe308) located at the bottom of the binding pocket in the transmembrane helix VI (2124). Although the binding mode of SR49059 into the V1aR is not yet defined, nonconserved residues Thr333 and Ala334 located in transmembrane region VII have been shown to control the V1aR/V2R receptor binding selectivity for SR49059 and for cyclic peptide antagonists as well (21). To precisely investigate the SR49059-V1aR binding interactions at a molecular level, we first constructed a three-dimensional model of the antagonist docked into the receptor based on the x-ray crystal structure of bovine rhodopsin (25). To validate or invalidate the hypotheses of potential ligand-receptor interactions, we developed a site-directed irreversible labeling strategy that combines mutagenesis of the receptor and binding of chemically reactive probes derived from SR49059. This approach generates a chemical bond between the nucleophilic moiety of a cysteine residue incorporated into the receptor and the electrophilic moiety of a sulfydryl-reactive group (isothiocyanate or an
-chloroacetamide) in the SR49059 analogue. Because such a chemical bond formation is only possible when a close contact exists between reactive elements created both on the ligand and on the protein, the covalent link can be taken as a proof of direct interaction versus a long range effect (26). This strategy combines the advantages of three-dimensional modeling and site-directed mutagenesis techniques with those of a direct unambiguous identification of contact regions between a receptor and a specific ligand, like in photoaffinity labeling studies. This is an elegant complement to receptor photolabeling; however, it is much easier to develop than photolabeling, because it does not require the use of a radiolabeled ligand and the overproduction of the target receptor. However, this strategy necessitates introduction of a mutation in the receptor, which could be by itself detrimental to the overall structure of the protein. In this study, we demonstrated that the use of cysteine mutants of V1aR as chemical sensors for sulfydryl-reactive ligands allows unambiguous localization of the SR49059 nonpeptide-binding site and offers a reliable molecular model of the V1a receptor.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Alignment of Amino Acid SequencesThe amino acid sequences of the three human AVP receptor subtypes and that of the OT receptor were retrieved from the Swiss-Prot data base (accession numbers were as follows: V1a receptor, P37288 [GenBank] ; V1b receptor, P47901 [GenBank] ; V2 receptor, P30518 [GenBank] ; and OTR, P30559 [GenBank] ) and aligned to the sequence of bovine rhodopsin (accession number P02699 [GenBank] ) using the ClustalW multiple alignment program (28). A slow pairwise alignment using BLOSUM matrix series (29) and a gap opening penalty of 15.0 were chosen for aligning the amino acid sequences to the sequence of bovine rhodopsin. Because the disulfide bridge occurring between the third transmembrane segment (TM III) and the second extracellular loop (E II) in the structure of bovine rhodopsin is conserved in all AVP/OT receptors, we manually adjusted the alignment of E II to align the corresponding cysteines.
Preparation of Starting Protein CoordinatesThe three-dimensional model of the human V1a receptor was constructed by mutating the side chains of the amino acids in rhodopsin. Standard geometries for the mutated side chains were given by the BIOPOLYMER module of SYBYL (SYBYL 6.8, Tripos Inc., St. Louis, MO). Whenever possible, the side chain torsional angles were kept to the values occurring in bovine rhodopsin. Otherwise, a short scanning of side chain angles was performed to remove steric clashes between the mutated side chain and the other amino acids. The third intracellular loop between helices V and VI, which shows a high degree of variability, was not included in any of the three models. This loop is not involved in direct interactions with the ligand (31). We therefore assume that omitting this loop should not influence our docking results. Insertions/deletions occurred only in loops but not in secondary structure elements (
-helix and
-sheet). The insertions/deletions in the loops were achieved through a simple knowledge-based loop search procedure using the LOOPSEARCH module of the SYBYL package. In this procedure, a set of 1478 high resolution x-ray structures was searched for loops of similar length and similar distance between the C
atoms of the residues delimiting the loop window. The loop showing the highest sequence identity and the lowest root mean square deviation was then selected for insertion in the model. Special caution had to be given with E II, which has been described in bovine rhodopsin to fold back over transmembrane helices (25) and thereby limit the size of the active site. Hence, amino acids of this loop could be involved in direct interactions with the ligands. A driving force to this peculiar fold of the E II loop might be the presence of a disulfide bridge between cysteines in TM III and E II. Because this covalent link is conserved in all receptors modeled in the current study, the E II loop was modeled using a rhodopsin-like constrained geometry around the E II-TM III disulfide bridge. After the heavy atoms were modeled, all of the hydrogen atoms were added, and the protein coordinates were then minimized with AMBER (32) using the AMBER95 force field (33). Topology and charge parameters for the nonpeptide V1a antagonist were computed using a previously described procedure (34). The minimizations were carried out by 1,000 steps of steepest descent followed by conjugate gradient minimization until the root mean square gradient of the potential energy was less than 0.05 kcal/mol·Å. A twin cut-off (10.0 and 15.0 Å) was used to calculate nonbonded electrostatic interactions at every minimization step, and the nonbonded pair list was updated every 25 steps. A distance-dependent (
= 4r) dielectric function was used.
Modeling the Antagonist-bound V1a Receptor ModelTo obtain an "antagonist-bound" model of the V1a receptor, the above-described coordinates were refined by the following procedure. Conivaptan, a known V1a antagonist (35) was first manually docked into the putative active site (21) so that its shape optimally fit in that of the binding pocket. The resulting protein-ligand complex was then refined by minimization using the above-described AMBER parameters. Removing the ligand atoms from the minimized complex finally yielded one set of coordinates for an antagonist-bound receptor model.
Automated Docking of a Selective V1a Antagonist (SR49059)The Gold 1.2 docking program (36) was used to automatically dock the selective V1a antagonist SR49059. For each of the 10 independent genetic algorithm (GA) runs, a maximum number of 1,000 GA operations was performed on a single population of 50 individuals. Operator weights for cross-over, mutation, and migration were set to 100, 100, and 0, respectively. To allow poor nonbonded contacts at the start of each GA run, the maximum distance between hydrogen donors and fitting points was set to 5 Å, and nonbonded Van der Waals energies were cut off at a value equal to the kij well depth of the van der Waals energy for the atom pair i and j). To further speed up the calculation, the GA docking was stopped when the top three solutions were within 1.5 Å root mean square deviation. If this criterion is met, we can assume that these top solutions represent a reproducible pose for the ligand.
Site-directed Mutagenesis of the Human V1aRPoint mutations F225C, Y300C, F308C, and K128C introduced in the cDNA sequence of the human V1aR vasopressin receptor were directly achieved on the eukaryotic expression vector pCMV (37) using the QuikChangeTM site-directed mutagenesis kit (Stratagene). All of the mutations were verified by direct dideoxynucleotide sequencing (T7TM sequencing kit; Amersham Biosciences).
Cell Culture, Receptor Expression, and Membrane Preparations The wild-type human V1aR and its mutants were transiently expressed in COS-7 cells or Chinese hamster ovary cells by electroporation as previously described (17, 22). Briefly, the cells were suspended in electroporation buffer (107 cells/0.3 ml) and incubated with plasmid DNA (20 µg of carrier DNA and 14 µg of expression vector containing the mutant cDNA insert) for 10 min at room temperature before being pulsed (280 V, 950 microfarads; GeneZapper system, Kodak Scientific Imaging). After electroporation, the cells were plated in Petri dishes and grown in Dulbecco's modified Eagle's medium (Invitrogen) supplemented with 10% fetal calf serum, 4 mM L-glutamine, 500 units/ml penicillin and streptomycin each, and 0.25 µg/ml amphotericin B in an environment containing 95% air and 5% CO2 at 37 °C. To increase the level of expression of some mutants, the cells were treated overnight with 5 mM sodium butyrate (38, 39). As already published, this treatment does not modify the pharmacological properties of the receptors (22, 23). Transfected cells were harvested 4872 h after electroporation, and the membranes were prepared as already described (17). In short, the cells were washed twice in phosphate-buffered saline without Ca2+ and Mg2+, Polytron-homogenized in lysis buffer (15 mM Tris-HCl, pH 7.4, 2 mM MgCl2, 0.3 mM EDTA), and centrifuged at 800 x g for 5 min at 4 °C. The supernatants were recovered and centrifuged at 44,000 x g for 20 min at 4 °C. The pellets were washed in Buffer A (50 mM Tris-HCl, pH 7.4, 5 mM MgCl2) and centrifuged at 44,000 x g for 20 min at 4 °C. The membranes were suspended in a small volume of Buffer A, and the protein content was determined by the Bradford method (Bio-Rad) using bovine serum albumin as the standard. Aliquots of membranes were used immediately for binding assays or stored at 80 °C.
Radioligand Binding AssaysThe binding assays were performed at 30 °C using [125I]HO-LVA or [3H]AVP as the radioligands and 13 µg (for assays with [125I]HO-LVA) or 1015 µg (for assays with [3H]AVP) of membrane proteins in standard radioligand saturation and competition binding assays as previously described (17). Briefly, the membranes were incubated in Buffer A supplemented with 1 mg/ml bovine serum albumin (binding buffer) and with radiolabeled and displacing ligands for 30 min (with [3H]AVP) or 60 min (with [125I]HO-LVA). Affinities (Kd) for [125I]HO-LVA (concentrations from 50 to 800 pM) as well as for [3H]AVP (concentrations from 0.1 to 20 nM) were directly determined in saturation experiments. Affinities (Ki) or apparent affinities (apparent Ki) for the unlabeled ligands were determined by competition experiments using [125I]HO-LVA (100500 pM) or [3H]AVP (12 nM) as the radioligands. The concentration of the unlabeled ligands varied from 1 pM to 10 µM. In saturation and competition experiments, depending on the radiolabeled peptide, nonspecific binding was determined by adding unlabeled AVP (10 µM) or HO-LVA (400 nM). Bound and free radioactivity were separated by filtration over Whatman GF/C filters presoaked in a 10 mg/ml bovine serum albumin solution (with [3H]AVP) or in a 0.5% polyethylenimine solution (with [125I]HO-LVA) for 34 h. The ligand binding data were analyzed by nonlinear least squares regression using the computer program Ligand (40). All of the assays were performed in triplicate on at least three separate batches of electroporated cells.
Chemical Modification of the Mutant Receptors and Measure of Residual Radioligand BindingTo covalently and irreversibly label the cysteine mutants of the human V1aR, membrane preparations of COS-7 cells were first incubated with sulfydryl-modifying antagonist ligands (100 nM-10 µM), their nonreactive analogues (110 µM), or vehicle (Me2SO) in binding buffer for 30 min at 37 °C. The protein content was controlled before this treatment. Next, the membranes were pelleted at 21,000 x g for 20 min at 4 °C, washed, and centrifuged four times with binding buffer without the ligands. The protein content of the pelleted and washed membranes was checked again before the radioligand binding assays. Residual binding of V1aR was assessed by incubating the homogenate with different concentrations of [125I]HO-LVA (from 50 to 800 pM) or [3H]AVP (from 0.1 to 20 nM) depending on the mutant (saturation experiments). Nonspecific binding was determined in parallel assays using saturating concentrations of HO-LVA (400 nM) or AVP (10 µM). Protein concentrations used in the receptor irreversible labeling assays are not equivalent to the ones used in classical binding assays (13 µg for assays with [125I]HO-LVA and 1015 µg for assays with [3H]AVP). Indeed, because several washes of the membranes are necessary for eliminating the irreversible ligand after preincubation of the receptors, more than half of the protein content is lost during this protocol. Thus, 2540 µg of proteins are preincubated with the ligand when using the F225C mutant. Only 510 µg of proteins are preincubated with the irreversible ligand when using the wild-type V1a receptor because of a receptor expression level much higher than that of the mutant. After preincubating the membranes with the irreversible ligand and then extensive washing, the protein content used for measuring [125I]HO-LVA residual binding is much lower, approximately 1015 µg/assay for the F225C mutant and 24 µg/assay for the wild-type V1a receptor. For instance, in the residual binding experiments illustrated in Fig. 6, protein contents are as follows: 3.7, 2.3, and 2.3 µg/assay for wild-type V1a-expressing membranes that have been pretreated with vehicle, compound 2, and compound 4, respectively; 12.3, 14.6, and 13.4 µg/assay for mutant F225C-expressing membranes that have been pretreated with vehicle, compound 2, and compound 4, respectively.
|
|
| RESULTS |
|---|
|
|
|---|
|
The sequence of the three human AVP receptors and that of the OT receptor were first aligned to that of bovine rhodopsin in a consistent manner. As seen in Fig. 2, class I GPCR-specific fingerprints (41) are all aligned at key positions of each of the seven transmembrane domains. Because the two cysteine residues forming a disulfide bridge between TM III and E II in the rhodopsin x-ray crystal structure (25) are conserved in the V1aR, we postulated the existence of the same disulfide linkage in the target receptor, although no experimental data yet support this hypothesis. Taking into account the alignment described on Fig. 2 and the x-ray crystal structure of bovine rhodopsin, a three-dimensional model of the V1aR was constructed as detailed under "Experimental Procedures." The model is suitable for an automated docking of the selective antagonist SR49059 and is compatible with most experimental data derived from studies undertaken to map the V1a antagonist-binding sites (2124, 42, 43).
|
As seen in Fig. 3, the dimethoxyphenyl group of SR49059 could be extensively buried among hydrophobic and aromatic amino acids (Phe225, Tyr300, Trp304, Phe307, and Phe308) of transmembrane segments V and VI. Mutation of some of these residues is detrimental to antagonist binding, particularly those of Phe307 and Phe308 (21). The nearby indoline moiety could also develop hydrophobic contacts to TM III, V and VI, but the corresponding subsite is partly hydrophilic. Two glutamine residues (Gln185 and Gln311) known to map the SR49059 binding site (17, 21) are in close proximity to the antagonist. Interestingly, the Lys128 side chain that has previously been shown to play a key role in SR49059 binding to the V1a receptor does not contribute by the polar nitrogen atom but by its apolar carbon atoms (Fig. 3A). Similar interactions have also been evidenced by x-ray diffraction analysis of major histocompatibility complex-peptide complexes (44). The present model also accounts for the selectivity of this compound toward the human V1a receptor. Hence, the two residues (Thr333 and Ala334) located on top of TM VII and previously reported to control V1a/V2 binding selectivity (21) directly interact with SR49059 in the present model
|
To validate the three-dimensional model described above, different analogues of SR49059 derivatized with sulfydryl-reactive moieties were synthesized. These compounds containing an isothiocyanate moiety or an
-chloroacetamide as chemical reactive entities as well as their precursors are shown in Fig. 4. All of the chemical substitutions were introduced on the dimethoxyphenyl ring.2 Assuming a conserved binding mode for close analogues of SR49059, the isothiocyanate moiety of compound 4 (or the
-chloroacetamide of compound 5) should be located close enough to three residues of the V1a receptor (Phe225, Tyr300, and Phe308) to afford covalent coupling with the respective cysteine mutants (F225C, Y300C, or F308C). These three residues were chosen to maximize chances of covalent linkage for three main reasons: 1) they are close enough to the dimethoxyphenyl moiety of SR49059, 2) their side chains are almost coplanar to the dimethoxyphenyl ring of the reference antagonist, and 3) their cysteine mutants could react with synthetically accessible reactive analogues of SR49059. We thus decided to construct these receptor mutants and to verify this hypothesis. Lys128, which is not supposed to directly interact with the isothiocyanate moiety of compound 4 (or the
-chloroacetamide of compound 5), was substituted with a cysteine residue as well as a negative control.
|
Affinities of the nonpeptide antagonist compounds were first determined for the V1a receptor and compared with those of the parent ligand SR49059. As seen in Table I, all of the chemical substitutions led to decreases in affinity, with the SR49059 being the most potent. Interestingly, the two sulfydryl-modifying compounds 4 and 5 still displayed a relative high affinity for the V1a receptor (6.3 and 16.3 nM, respectively), which is a good prerequisite for investigating their capacity to establish a covalent binding to the cysteine mutants incorporated into the receptor.
|
In a second step, the affinities of the SR49059 and its nonpeptide analogues to the cysteine mutants of the V1a receptors were calculated. The results are also indicated in Table I. We could not express the Y300C V1a mutant. Thus, only binding properties of mutants F225C, F308C, and K128C are displayed in Table I. We first measured the affinity of the radiolabeled probes, [3H]AVP and [125I]HO-LVA. For the mutant F308C, [3H]AVP was used as the radioligand, and its Kd was 0.85 ± 0.18 nM (5), a value equivalent to the one determined for the wild-type receptor, 0.67 ± 0.17 nM (21, 22). We could not determine the affinity of [125I]HO-LVA for this mutant, suggesting a significant decrease in affinity for this compound. On the contrary, mutant receptors F225C and K128C cannot bind [3H]AVP, and thus [125I]HO-LVA was used as the radioligand. [125I]HO-LVA could bind to F225C and K128C mutants with Kd values equal to 160 ± 37 pM (n = 6) and 420 ± 140 pM (n = 3), respectively. These values are in the range of that measured for the wild-type V1a receptor, Kd = 50 ± 4 pM (21, 22). Because compounds 4 and 5 are supposed to potentially create a covalent bond with the mutant receptors and because affinity of an irreversible ligand cannot be determined reliably in a conventional competition binding assay, only apparent Ki values instead of true Ki values were measured for these two nonpeptide ligands. The apparent affinity determined in this way may actually represent both reversible and irreversible binding. Thus, simultaneous covalent and competitive binding should overestimate the apparent Ki for the mutant receptors.
In the case of each receptor mutant, affinities for the SR49059 and for its analogues were all diminished (Table I). For the F225C V1aR mutant, the Ki values for SR49059 and compounds 25 are decreased three to six times when compared with that obtained with the wild-type receptor, leading to concentrations still in the 10100 nM range. The moderate decrease in the binding affinity of SR49059 for the F225C V1a receptor mutant can be easily explained from the current molecular model. Hence, Phe225 is only one of the eight amino acids delimiting a strong hydrophobic pocket interacting with the dimethoxyphenyl moiety of SR49059 (Fig. 3B). The single point mutation does not significantly modify the polarity and the physicochemical property of that subsite. Docking SR49059 into the three-dimensional model of F225C mutant alters neither the binding mode (root mean square deviation of 0.25 Å from the wild-type-bound SR49059 coordinates) nor the GOLD docking score (63.52 for wild type versus 60.49 for the F225C mutant). For F308C V1aR and K128C V1aR mutants, the affinities are much lower, being in the 1001000 nM range. Like for photoactivatable ligands (22, 23), which can be covalently bound to their receptors upon UV irradiation, the affinity for sulfydryl-modifying affinity ligands is expected to be high to favor such covalent chemical reactions. Taking into account the previous results, using compound 4 (or compound 5) in combination with the F225C V1aR mutant receptor and at a much lesser extent with F308C and K128C V1aR mutants, thus represented a particularly favorable situation for applying the site-directed irreversible labeling approach.
In the next step, potential irreversible binding of the sulfydryl-reactive analogues of SR49059 was verified on the different cysteine mutants of the human V1a receptor. To facilitate the covalent chemical modification of the receptor mutants, compound 4 was first used at a very high concentration (105M) during a 30-min incubation with COS cell membranes. As described under "Experimental Procedures," residual binding of [125I]HO-LVA or [3H]AVP was measured after extensive washing of the COS membranes. A negative control (Me2SO as vehicle only) was included in the assays to ascertain the specificity of the effects. The residual binding gives an idea about the number of receptor sites still available for the radioligand and consequently the proportion of receptor sites irreversibly occupied by analogue 4. Different concentrations of [125I]HO-LVA or [3H]AVP were used to be able to construct Scatchard representations of the experiments and to measure Kd and Bmax parameters. The results are shown in Table II. As demonstrated, treatment of the membranes expressing mutants F308C or K128C by the sulfydryl-reactive analogue did not lead to any change in the binding properties of the radioligand ([125I]HO-LVA or [3H]AVP, depending on the mutant tested). Neither the affinity nor the number of sites (Bmax) were modified. These two parameters were equivalent to those measured in the corresponding control assays with the vehicle (Me2SO) only. On the other hand, treatment of the membranes expressing the F225C mutant with the analogue 4 had a drastic effect on the binding of the [125I]HO-LVA radioligand. Compared with the corresponding negative control (vehicle only) in which the Kd and Bmax values were 0.15 ± 0.02 nM and 0.15 ± 0.016 pmol/mg of membrane proteins, respectively, some residual binding was measured (1030%), but we were unable to calculate the Kd and Bmax parameters in this condition (Table II). This result suggested an irreversible and almost complete binding of analogue 4.
|
The results obtained with F225C mutant V1a receptor were very encouraging and led us to check whether the effect of the sulfydryl-reactive analogue was dose-dependent. The membranes expressing the F225C mutant were treated with different concentrations of analogue 4 and tested for [125I]HO-LVA residual binding. As illustrated in Fig. 5, a dose-dependent effect of the nonpeptide antagonist was unambiguously measured. At 100 nM, the SR49059 analogue was almost without effect because residual binding of [125I]HO-LVA represented 94.2 ± 2.9% of that of the control condition (membranes treated with vehicle only). At 1 µM, the residual binding was decreased, being 44.7 ± 6.9%, whereas at 10 µM, binding of the radioligand only represented 26.6 ± 3.4%. Concentrations of compound 4 required to irreversibly inactivate more than 50% of receptors are higher than measured Ki affinity values, suggesting that the covalent chemical reaction between the two partners is slower than their association/dissociation. One should also remember that, as explained above, because affinity of an irreversible ligand cannot be determined reliably in a conventional competition binding assay, only an apparent Ki value (25 nM here for the F225C mutant) instead of a true Ki value was measured for this nonpeptide ligand. The apparent affinity determined in this way may actually represent both reversible and irreversible binding. Thus, simultaneous covalent and competitive binding should overestimate the apparent Ki for the mutant receptors. Consequently, the true Ki value for this compound could be much higher, in the 100 nM range or even worse.
The specificity of the irreversible binding reaction using the wild-type V1a receptor and the nonreactive analogue 2 (which is also the precursor of analogue 4) as two negative controls was then investigated. To be able to quantify residual binding (Kd and Bmax values) for [125I]HO-LVA after treatment of the membranes, we used analogue 4 at 1 µM. This concentration has been demonstrated to inhibit approximately 50% of the radio-labeled HO-LVA ligand binding (Fig. 5). The membranes expressing either the wild type or the F225C mutant were treated with vehicle, analogue 2 (1 µM), or analogue 4 (1 µM). Representative Scatchard plots of the residual binding of [125I]HO-LVA are illustrated in Fig. 6. As expected for the wild-type receptor, analogues 2 and 4 were without effect on the binding of [125I]HO-LVA, and Scatchard representations were almost superimposable with that of the vehicle condition (Kd were 46.5, 77, and 77 pM for vehicle, analogue 2, and analogue 4, respectively; corresponding Bmax values were 1.05, 1.18, and 0.94 pmol/mg of membrane protein). On the contrary, treatment of the membranes expressing the F225C mutant with the sulfydryl-reactive analogue 4 significantly decreased the number of sites for [125I]HO-LVA (76 compared with 142 fmol/mg of membrane protein for the control condition), whereas analogue 2 was without effect (166 versus 142 fmol/mg of membrane protein). In the three conditions, Kd values for the radiolabeled antagonists were equivalent; values were 188, 196, and 228 pM for vehicle, analogue 2, and analogue 4, respectively. Pretreatment of the membranes expressing either the wild-type V1a or the mutant F225C with 10 µM analogue 2 was without effect on [125I]HO-LVA binding (data not shown). In conclusion, analogue 4 (R = NCS) is able to bind irreversibly and specifically to mutant F225C in a dose-dependent manner. Unfortunately, these results could not be reproduced with a high concentration (10 µM) of analogue 5, which carries an
-chloroacetamide as a sulfydryl-reactive group, whereas analogue 4 carries an isothiocyanate moiety (data not shown).
| DISCUSSION |
|---|
|
|
|---|
-chloroacetamide instead of an isothiocyanate moiety as a sulfydryl-reactive group could not irreversibly label the F225C human V1a receptor mutant, suggesting that an optimal orientation of the sulfydryl group of the cysteine residue and the electrophilic moiety of the reactive ligand is necessary to afford such a covalent bond. Moreover, the
-chloroacetamide group is much more bulky than an isothiocyanate moiety, a steric parameter that could be crucial as well. Taking into account the result of affinity-directed irreversible binding, the model is in agreement with previous experimental results. First, we have already demonstrated that the SR49059 could be in contact with residues Thr333 and Ala334 located at the top of transmembrane domain VII (21). Based on a gain-of-function mutagenesis strategy, these two amino acids were shown to control the V1a/V2 receptor subtype binding selectivity of the nonpeptide antagonist. This result and the one described in this study make the three-dimensional model of the SR49059-V1a receptor complex much more reliable. In addition, based on the three-dimensional model, Lys128 (TM III), Gln185 (TM IV), Trp304, Phe307, and Phe308 (TM VI) could also participate in the binding site of the nonpeptide ligand, an hypothesis compatible with previous site-directed mutagenesis results. Indeed, mutation of these residues was detrimental to the antagonist binding (1719, 2124). In conclusion, SR49059 seems to bind deeply in the transmembrane region of the human V1a receptor, surrounded by a hydrophobic/aromatic cage. The involvement of Phe225 confirmed the importance of aromatic residues in the binding of all classes of V1a receptor-specific antagonists, both cyclic and linear peptides as well as nonpeptides (2124).
The binding sites for SR49059 and AVP only partially overlap. This overlap is constituted with polar residues such as Lys128 and Gln185 (21, 24) and to a certain degree Gln131 and Gln311 (17). On the contrary, although AVP and agonists interact with extracellular residues of the receptor such as Tyr115 in the extracellular loop I, which is responsible for receptor subtype selectivity (19), or residue Arg46 in the N-terminal part (20), SR49059 does not. In addition, Phe225 does not play any role in the binding of agonists (21), but the homologous residue in the OTR (Tyr209) modulates the response of this receptor to the partial agonist AVP (45). AVP was demonstrated to become a full agonist on OTR when this aromatic residue was replaced by the residue present at the equivalent position in the V1a receptor (F225). Moreover, transmembrane V, in which Phe225 is located, seems to play an important role for other nonpeptide V1a antagonists. It has been experimentally demonstrated that Ile224, adjacent to Phe225, controls the rat/human V1a receptor selectivity toward the OPC21268 nonpeptide antagonist binding (42). Based on molecular modeling hypotheses, transmembrane domain V could also interact with the nonpeptide YM087, a potent V1a/V2 antagonist (43). Finally, it is now commonly accepted that Phe225, like other conserved hydrophobic/aromatic residues in the transmembrane domains of all GPCRs, participates into an intramolecular network of interactions that maintains the receptors as a constrained/inactive conformation (46). SR49059, while interacting with Phe225, could reinforce this network of intramolecular interactions, thus behaving as an antagonist.
Cysteine scanning has been widely used to allow the incorporation of chemical probes into membrane receptors (26). However, incorporating such a cysteine mutant to be labeled with a sulfydryl-reactive affinity ligand has been applied very few times, for example in the case of the
-aminobutyric acid, type A receptor with noncompetitive channel blockers (47) or in the case of the human
2A-adrenergic receptor (30). For this GPCR, chloroethylclonidine, an alkylating derivative of the
2A-adrenergic agonist clonidine, bound irreversibly to the receptor by forming a covalent bond with the sulfydryl side chain of a cysteine residue exposed in the binding cavity, leading to inactivation of the receptor. Using cysteine mutants as chemical sensors for interacting with sulfydryl-reactive affinity ligands is a powerful strategy to investigate ligand-receptor interactions. This approach is complementary to molecular modeling, site-directed mutagenesis and photoaffinity labeling. Like photoaffinity labeling, this approach results in unambiguous identification of residues directly involved in the binding of the ligands and in the precise definition of the binding pockets of GPCRs. Interestingly, this approach does not require the use of a radiolabeled ligand, the overproduction of the target protein, and its tedious sequencing to locate the labeled residue. However, this strategy necessitates a reliable three-dimensional model of the hypothetical ligand-receptor interactions and the introduction of a mutation in the receptor, which could itself be detrimental to the overall structure of the protein.
AVP and OT receptors represent a logical target for drug development, and the strategy we described in the present paper would be of high interest for defining the binding sites of nonpeptide antagonists and agonists selective to the vasopressin V2 and V1b receptors or to the OTR. Such characterization should facilitate the rational design of potent and selective therapeutic molecules useful in many physiopathological conditions.
| FOOTNOTES |
|---|
These authors contributed equally to this work. ![]()
|| To whom correspondence should be addressed: Laboratoire de Pharmacochimie de la Communication Cellulaire, UMR CNRS 7081, Faculté de Pharmacie de Strasbourg, 74 Route du Rhin, BP 24, 67401 Illkirch Cedex, France. Tel.: 33-3-90-24-42-32; Fax: 33-3-90-24-43-10; E-mail: mhibert{at}pharma.u-strasbg.fr.
1 The abbreviations used are: AVP, arginine vasopressin; V1aR, vascular V1a vasopressin receptor; V2R, renal V2 vasopressin receptor; V1bR, pituitary V1b vasopressin receptor; OT, oxytocin; OTR, OT receptor; HO-LVA, HO-phenylacetyl1-D-Tyr(Me)2-Phe3-Gln4-Asn5-Arg6-Pro7-Arg8NH2; TM, transmembrane segment; E, extracellular loop; GA, genetic algorithm; GPCR, G protein-coupled receptor. ![]()
2 C. Tahtaoui, M. N. Balestre, P. Klotz, D. Rognan, C. Barberis, B. Mouillac, and M. Hibert, manuscript in preparation. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|