HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

J. Biol. Chem., Vol. 275, Issue 31, 24003-24012, August 4, 2000

This Article Abstract Full Text (PDF) All Versions of this Article: 275/31/24003    most recent M002325200v1 Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Nicole, P. Articles by Laburthe, M. Search for Related Content PubMed PubMed Citation Articles by Nicole, P. Articles by Laburthe, M. Social Bookmarking                      What's this?

Identification of Key Residues for Interaction of Vasoactive Intestinal Peptide with Human VPAC₁ and VPAC₂ Receptors and Development of a Highly Selective VPAC₁ Receptor Agonist

ALANINE SCANNING AND MOLECULAR MODELING OF THE PEPTIDE^*

Pascal Nicole, Laurence Lins, Christiane Rouyer-Fessard, Cyrille Drouot^§, Pierre Fulcrand^§, Annick Thomas, Alain Couvineau, Jean Martinez^§, Robert Brasseur^¶, and Marc Laburthe

From  Unité INSERM U410 de Neuroendocrinologie et Biologie Cellulaire Digestives, Faculté de Médecine Xavier Bichat, Paris, 75018, France, the ^§ UMR CNRS 5810, Universités Montpellier I and II, Faculté de Pharmacie, Montpellier, 34060, France, and the ^¶ Centre de Biophysique Moléculaire Numérique, Faculté des Sciences Agronomiques, B-5030 Gembloux, Belgium

Received for publication, March 15, 2000, and in revised form, May 3, 2000

	ABSTRACT

TOP ABSTRACT INTRODUCTION EXPERIMENTAL PROCEDURES RESULTS DISCUSSION REFERENCES

The widespread neuropeptide vasoactive intestinal peptide (VIP) has two receptors VPAC₁ and VPAC₂. Solid-phase syntheses of VIP analogs in which each amino acid has been changed to alanine (Ala scan) or glycine was achieved and each analog was tested for: (i) three-dimensional structure by ab initio molecular modeling; (ii) ability to inhibit ¹²⁵I-VIP binding (K_i) and to stimulate adenylyl cyclase activity (EC₅₀) in membranes from cell clones stably expressing human recombinant VPAC₁ or VPAC₂ receptor. The data show that substituting residues at 14 positions out of 28 in VIP resulted in a >10-fold increase of K_i or EC₅₀ at the VPAC₁ receptor. Modeling of the three-dimensional structure of native VIP (central -helice from Val⁵ to Asn²⁴ with random coiled N and C terminus) and analogs shows that substitutions of His¹, Val⁵, Arg¹⁴, Lys¹⁵, Lys²¹, Leu²³, and Ile²⁶ decreased biological activity without altering the predicted structure, supporting that those residues directly interact with VPAC₁ receptor. The interaction of the analogs with human VPAC₂ receptor is similar to that observed with VPAC₁ receptor, with three remarkable exceptions: substitution of Thr¹¹ and Asn²⁸ by alanine increased K_i for binding to VPAC₂ receptor; substitution of Tyr²² by alanine increased EC₅₀ for stimulating adenylyl cyclase activity through interaction with the VPAC₂ receptor. By combining 3 mutations at positions 11, 22, and 28, we developed the [Ala^11,22,28]VIP analog which constitutes the first highly selective (>1,000-fold) human VPAC₁ receptor agonist derived from VIP ever described.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION EXPERIMENTAL PROCEDURES RESULTS DISCUSSION REFERENCES

The vasoactive intestinal peptide (VIP)¹ is a prominent neuropeptide with wide distribution in both peripheral and central nervous systems and a large spectrum of biological actions in mammals (1, 2). VIP-containing nerves and VIP effects have been described in digestive tract, cardiovascular system, airways, reproductive system, immune system, endocrine glands, and brain (1). Besides its short-term actions on exocrine secretions, hormone release, muscle relaxation, and metabolism (1, 2), VIP has been also characterized as a growth regulator for fetuses and tumor cells and during embryonic brain development (3). There are recent evidences for an important role of VIP in the perception of pain (4) and suppression of inflammation (5). Finally, VIP has been involved in diseases such as the watery diarrhea syndrome and clinical applications of VIP have been already suggested in impotence, asthma, lung injury, a variety of tumors and neurodegenerative diseases (1-3).

VIP belongs to a large family of structurally related peptides (2, 6, 7) that comprises VIP, pituitary adenylate cyclase-activating peptide PACAP-27, and its C-terminal extended form PACAP-38, secretin, glucagon, and glucagon-like peptides-1 and -2, gastric inhibitory polypeptide, peptide histidine methionine amide, growth hormone-releasing factor (GRF), and peptides isolated from the venom of the Gila Monster. VIP and PACAP are the most closely related peptides in terms of structure and function (2, 6). They share two common receptors, VPAC₁ and VPAC₂, which display high affinity for both VIP and PACAP (2, 8). These receptors together with receptors for VIP-related peptides (see above) clearly constitute an original subfamily within the superfamily of G protein-coupled receptors (2, 9, 10). This subfamily referred to as class II (2) also comprises receptors for parathyroid hormone, calcitonin, corticotropin-releasing factor, and the so called EGF-TM7 receptors (11). Class II family of receptors for peptides display several common properties including large N-terminal extracellular domains containing highly conserved cystein residues, N-terminal leader sequences, and complex gene organization with many introns (2).

Although the structure-function relationship of VIP receptors, including VPAC₁ and VPAC₂, has been recently documented (2, 9, 12-20), the structure-function relationship of VIP itself is still poorly understood. Some old studies carried out before the characterization and cloning of VIP receptor subtypes (21-23) indicated that: (i) the entire sequence of VIP is required for full biological function (21, 22); (ii) VIP-related peptides having significant sequence homologies with VIP such as peptide histidine methionine amide, secretin, GRF, and helodermin behave as low potency VIP agonists (23); (iii) there are important differences between species in the pharmacology of VIP receptors, especially between rodents and humans (2, 15). In this context, the present study explores the structure-function relationship of VIP for interacting with the two human VIP receptor subtypes VPAC₁ and VPAC₂. The contribution of each side chain of VIP was investigated by the systematic single exchange of each residue of VIP 1-28 by alanine or glycine. All modified VIP peptides were then analyzed for their binding affinities and abilities to stimulate adenylyl cyclase in CHO cell clones stably expressing human recombinant VPAC₁ or VPAC₂ receptor. Moreover, the contribution of each residue was also analyzed after molecular modeling of all VIP analogs. These studies allowed us to determine key residues for interaction with human VPAC₁ and VPAC₂ receptors and also resulted in the development of the most highly selective VPAC₁ receptor agonist ever described.

	EXPERIMENTAL PROCEDURES

TOP ABSTRACT INTRODUCTION EXPERIMENTAL PROCEDURES RESULTS DISCUSSION REFERENCES

Materials-- Enzymes and vectors for cloning were obtained from Promega (Charbonnière, France). The recombinant human VPAC₁ receptor was stably expressed in CHO cells as described (24). The human VPAC₂ cDNA (25) was a gift from Dr. M. Svoboda and Dr. P. Robberecht (Brussels, Belgium). Culture medium was obtained from Life Technologies, Inc. (Cergy-Pontoise, France). In addition to VIP analogs synthesized by us (see below), [Gly⁴]VIP, [Ala¹²]VIP, [Ala²³]VIP, and [Ala^11,22,28]VIP analogs were obtained by custom peptide synthesis from Neosystem (Strasbourg, France). ¹²⁵I-VIP was prepared and purified as described (26). Other highly purified chemicals used were from Sigma (Saint-Quentin-Fallavier, France).

Synthesis, Purification, and Analysis of VIP Analogs-- The peptides were synthesized by solid phase on a PerSeptive Biosystems Pioneer Instrument using a Fmoc strategy. We have used a Fmoc-PAL-PEG-PS resin (27) (PerSeptive Biosystems GmbH, Hamburg, Germany) loaded at 0.16 mmol/g. N-protected Fmoc amino acids were used. For functionalized amino acids, the following derivatives were used: Fmoc-Asp(OBut)-OH, Fmoc-Glu(OBut)-OH, Fmoc-Arg(Pmc)-OH, Fmoc-His(Trt)-OH, Fmoc-Ser(But)-OH, Fmoc-Thr(But)-OH, Fmoc-Tyr(But)-OH, Fmoc-Gln(Xan)-OH, Fmoc-Asn(Trt)-OH, and Fmoc-Lys(Boc)-OH and were from Advanced ChemTech (Machelen, Belgium). The coupling reagent was 2-(1H-benzotriazol-1-yl)-1,1,3,3-tetramethyluronium hexafluorophosphate (28) for all peptides except for [Ala⁵]VIP, [Ala¹⁶]VIP, [Ala¹⁷]VIP, and [Ala²⁷]VIP for which HATU (29) was used. Syntheses of the different peptides were performed from 0.625 g of Fmoc-PAL-PEG-PS resin. An excess of 4 eq of each amino acid was used.

A double coupling was performed when Fmoc-Asn(Trt)-OH, Fmoc-Val-OH, Fmoc-Gln(Xan)-OH, Fmoc-Arg(Pmc)-OH, Fmoc-Thr(But)-OH, and Fmoc-Tyr(But)-OH were involved in coupling. Fmoc deprotection was performed with a solution of piperidine in dimethylformamide in a 2:8 (v/v) ratio. Final deprotection of the peptides from the resin was performed in a mixture containing trifluoroacetic acid, water, thioanisole, ethylenedithiol, phenol in a 82.5/5/5/2.5/5 (v/v) ratio (30) for 2 h. Peptides were then precipitated by addition of cold diethyl ether and dissolved in a mixture of 0.1% trifluoroacetic acid in water/acetonitrile and lyophilized. They were purified by high performance liquid chromatography on a Waters instrument equipped with a Waters Delta Packs C18 column (340 × 100 mm) using a gradient of solvents composed of water, 0.1% trifluoroacetic acid, and acetonitrile containing 0.1% trifluoroacetic acid. Purity of all peptides was checked by analytical high performance liquid chromatography on a Waters instrument using a C18 column (Novarpack, 5 µm, 100 Å, 180 × 3.5 mm) and all were at least 97% pure (UV detection at 214 and 254 nm). They were characterized by electrospray mass spectrometry using a Jeol LX 2000 spectrometer.

Molecular Modeling of VIP and VIP Analogs-- The conformations of native VIP and VIP analogs were calculated by molecular modeling using the OSIRIS software which is an ab initio approach allowing the calculation of three-dimensional structure of soluble proteins of less than 150 residues (31). The first step of OSIRIS is used in the present study since it is sufficient to correctly determine most secondary structures of different -helix proteins (32). Briefly, the primary sequence of the peptide is introduced and 179 pairs of  and  angles are successively attributed to all amino acids of the sequence.

For each rotation axis, the energy is calculated using the following equation,

(Eq. 1)

where E(k) is the energetic component associated to the torsion axis k, between atoms w and w+1, with E(k)_tor as torsional energy, E_vdw^ij as Van der Waals energy, E_elec^ij as electrostatic energy, E_sol.in^ij as internal solvation energy between atoms i and j, and E_sol.out^ij as solvent solvation energy, and q as an energy quantum equal to the interaction energy between atoms i and j, divided by the number of torsion axes between atoms i and j. A structure is then obtained from the linear combination of each residue energy minimum. The energy of this construct is calculated and the residues having the lowest energy are frozen. Then the calculation is started again, keeping the structure of frozen residues constant and testing the 179 pairs of ,  angles on the others. At each step of calculation, µ residues are frozen, µ varying between 5 and 28. Increasing the µ value allows to test the variability of the calculated VIP or VIP analog structures. Calculations were performed on RAMSES (Rapid Analysis Master/Slaves Extensible System), a parallel hardware of 21 tracer Europa Pentium Pro PC cadenced at 180 Mhz connected by a 100 Mbytes Network and controlled by a HP Vectra VA Pentium Pro cadenced at 200 Mhz. Molecular views were drawn with WinMGM 1.0 (Ab Initio Technology, Obernai, France) as described (33).

Statistical analysis of the predicted structures of VIP and VIP analogs was made using the SICLA program (32). This procedure starts from n × n distance tables calculated from the local root mean square distances of all pairs of structures obtained using different E_mvi (µ) values. The root mean square distances was calculated by fitting the backbone of 5 amino acids of a structure to the 5 corresponding residues of the reference structure (here, the predicted structure obtained with µ = 5). The root mean square distances of the fragment is attributed to the central amino acid, the window is then moved by 1 residue along the sequence and calculations are repeated. The clustering procedure used in SICLA is an automatic nonhierarchical classification that compares pairs of structures (34, 35). This procedure is described in detail elsewhere (32). From the clustering of structures calculated with different µ values, a mean distance between the central structure and the others is obtained. This distance is representative of the peptide structural variability.

Stable Expression of cDNA Encoding Human VPAC₁ and VPAC₂ Receptors in CHO Cells-- Full-length VPAC₁ or VPAC₂ receptor cDNAs were ligated in the EcoRI site of pcDNA3 (Invitrogen). Both receptors were tagged at the C terminus with a marker dodecapeptide (Tag) as described (14). The recombinant plasmid encoding human VPAC₂ was transfected into CHO cells (ECACC 85050302) by electroporation. Briefly, four million CHO cells, cultured as described (14), were preincubated on ice for 5 min with 20 µg of plasmid DNA and 20 µg of salmon sperm DNA carrier in cold Ham's F-12 medium. After electroporation (330 V, 1,000 microfarads, infinite resistance), cells were kept on ice for 5 min and then transferred to Petri dishes containing 10 ml of culture medium (Ham's F-12, 10% (v/v) fetal calf serum, 100 IU/ml penicillin, 100 µg/ml streptomycin). After 48 h, cells were selected by addition of geneticin (G418) at the final concentration of 800 µg/ml for 2 weeks. Resistant cells were cloned by limiting dilution as described (24). A clone, referred to as clone 10, was isolated and selected on the basis of its binding parameters for VIP (see below). Clone 10 cells were grown in the above described culture medium containing 100 µg/ml G418 in a humidified atmosphere of 95% air and 5% CO₂ at 37 °C. A clone (referred to as clone 15) of CHO cells stably expressing the human VPAC₁ receptor was previously isolated and characterized (24). Both clone 10 and clone 15 cells were passaged every 7 days in 25-cm² plastic culture flasks and used between passages 8 and 30, under conditions for which the density of VPAC₂ and VPAC₁ receptors did not vary significantly (see below).

Membrane Preparation-- Cells, either clone 10 or clone 15, were grown to confluency for 3-4 days. After removing the culture medium, attached cells were washed three times with phosphate-buffered saline, then harvested with a rubber policeman and centrifuged for 10 min at 2,000 × g. The cell pellet was exposed for 30 min at 4 °C to hypoosmotic 5 mM Hepes buffer, pH 7.4, and a particulate fraction, referred to as membrane preparation, was obtained as described (17). Membranes were stored at 80 °C until use. Protein content was measured by the procedure of Bradford (36) with bovine serum albumin as standard.

Receptor Binding Assay-- Membranes (100 µg of protein/ml) were incubated at equilibrium for 60 min at 30 °C with 0.05 nM ¹²⁵I-VIP and increasing concentrations of native VIP or VIP analogs. The incubation medium was 20 mM Hepes, pH 7.4, containing 2% (w/v) bovine serum albumin and 0.1% (w/v) bacitracin. The reaction was stopped by transferring the incubation medium to an excess of ice-cold buffer. Bound and free peptides were immediately separated by centrifugation at 20,000 × g for 10 min, and membrane pellets were washed twice with 10% (w/v) sucrose in ice-cold 20 mM Hepes. The radioactivity was then counted with a -counter. Specific binding was calculated as the difference between the amount of ¹²⁵I-VIP bound in the absence (total binding) and presence (nonspecific binding) of 10 µM unlabeled native VIP. The nonspecific binding represented about 10 and 15% of total radioactivity bound to clone 15 and clone 10 cell membranes, respectively. All binding data were analyzed using the LIGAND computer program (38). The dissociation constant (K_d) and binding capacity (B_max) of VIP binding to clone 10 and clone 15 cell membranes were determined by Scatchard analysis. From the linear Scatchard plots, the K_d of VPAC₁ receptor in clone 10 cell membranes and VPAC₂ receptor in clone 15 cell membranes were 0.4 and 0.7 nM, respectively. The B_max were 1.6 and 1.3 pmol/mg of protein, respectively. The constants K_i for the inhibition of ¹²⁵I-VIP binding by unlabeled peptides were calculated using the Cheng-Prusoff equation as described (39). All competition curves for VIP and VIP analogs fit with a monophasic binding profile for both VPAC₁ and VPAC₂ binding assays.

Adenylyl Cyclase Assay-- Adenylyl cyclase activity in membranes from clone 10 or clone 15 cells was assayed in the presence of increasing concentrations of native VIP or VIP analogs as described (37). Dose-response curves were fitted and concentration of peptide giving half-maximal response (EC₅₀) were calculated using the Prism software suite (GraphPad Software, San Diego, CA).

	RESULTS

TOP ABSTRACT INTRODUCTION EXPERIMENTAL PROCEDURES RESULTS DISCUSSION REFERENCES

Alanine Scanning Analysis of VIP for Interacting with Human VPAC₁ Receptor-- Twenty-six analogs of VIP were synthesized in which each side chain of VIP 1-28 was systematically replaced with alanine. Two other analogs were also synthesized in which glycine was substituted for alanine in position 4 or 18 of native VIP. The ability of each analog to interact with the recombinant human VPAC₁ receptor expressed in CHO clone 15 cells was first analyzed by competitive inhibition of ¹²⁵I-VIP binding to cell membranes. For all VIP analogs, the dose-response curves for inhibiting ¹²⁵I-VIP binding were parallel to that of native VIP and complete inhibition could be observed at 10 µM peptide concentration (not shown). All competitor curves fitted with a monophasic binding profile consistent with a single binding site. Analysis of K_i (Table I) showed that substituting residues at 14 positions including positions 1, 3, 5, 6, 7, 8, 10, 12, 14, 15, 20, 21, 23, and 26 of VIP resulted in an important decrease in affinity for VPAC₁ receptors, i.e. more than one log. Fig. 1 shows the competition curves for ¹²⁵I-VIP binding to VPAC₁ receptor for some of these VIP analogs with decreased affinity e.g. H1A, F6A, K20A, and I26A VIP analogs. Substitution at other positions resulted in no change of affinity or a small decrease of affinity of less than one log. The most pronounced decreases of affinity occurred with H1A, F6A, R12A, R14A, and L23A analogs for which a 100-200-fold decrease was observed. Other significant decreases of affinity were observed with D3A, V5A, T7A, D8A, Y10A, K15A, K20A, K21A, and I26A analogs, e.g. >10-fold. All VIP analogs were also tested for their ability to stimulate adenylyl cyclase activity in membranes of CHO clone 15 cells expressing the human VPAC₁ receptor. The dose-response curves for all analogs paralleled that of VIP and 10 µM analogs achieved the same maximal stimulation of enzyme activity as native VIP (not shown). Fig. 1 shows the dose-response curves for stimulating enzyme activity for some of these analogs with decreased potency. In general, there is a good correlation between the EC₅₀ for stimulating enzyme activity and the K_i for inhibiting ¹²⁵I-VIP binding (see Table I) (Fig. 2). All analogs behaved as VPAC₁ receptor agonists with identical or lower potencies than native VIP.

View larger version (28K): [in this window] [in a new window]   Fig. 1.   Binding assay and adenylyl cyclase assay of VIP analogs in CHO cell clone stably expressing human VPAC1 receptor. A, dose effects of native VIP and some VIP analogs for inhibition of 125I-VIP binding to membranes from CHO cells expressing VPAC1 receptor. B, stimulation of adenylyl cyclase activity by native VIP and some VIP analogs in membranes from CHO cells expressing VPAC1 receptor. The data are expressed as the percentage of maximal stimulation above basal obtained with 10 µM native VIP. All data are mean ± S.E. of at least three experiments performed in duplicate. , native VIP; , H1A mutant; , F6A mutant; , K20A mutant; , I26A mutant.

View larger version (18K): [in this window] [in a new window]   Fig. 2.   Correlation between the EC50 values and Ki values determined for native VIP and all single-substituted analogs of VIP in the adenylyl cyclase assay and binding assay in CHO cell clone expressing human VPAC1 receptor. All analogs are numbered according to the position of the amino acid that was substituted. Wt, native VIP. See Table I and legend to Fig. 1 for details.

Altogether, the data were consistent with the importance of many residues along the VIP 1-28 sequence. One can notice the particular contributions of the N-terminal 1-8 sequence and also of the central part of the molecule with five crucial basic residues i.e. Arg¹², Arg¹⁴, Lys¹⁵, Lys²⁰, and Lys²¹. Two hydrophobic residues in the C-terminal part also play a role, i.e. Leu²³ and Ile²⁶. A recurrent issue for explaining alanine scanning data for peptide ligand is to determine whether a given substitution by alanine disrupts a specific interaction between peptide and receptor and/or alters the global structure of the peptide. To address this important question, modeling of three-dimensional structure of all VIP analogs was performed.

Molecular Modeling of VIP and VIP Analogs-- The conformation of native VIP as well as those of all VIP analogs were analyzed using OSIRIS (see "Experimental Procedures"). The first step of this ab initio procedure allows to assign a secondary structure to each residue. Iterative calculations were performed and at each iteration µ residues are frozen up to the entire sequence, µ varying between 5 and 28 (µ = 5, 7, 9, 11, 12, 16, 20, and 28 residues). This mimics nucleation centers observed during the folding of a protein and was used here to estimate the structural variability of native VIP and each analog. Fig. 3 shows the different conformers obtained for native VIP. The conformational variability is essentially located at the N- and C-terminal extremities, the central domain from Val⁵ to Asn²⁴ being clearly helical. A statistical analysis was then carried out to assess the structural dispersion shown in Fig. 3. A statistical mean distance (D) based on the local root mean square between all pairs of conformers was calculated using SICLA (see "Experimental Procedures"). The structural dispersion of several VIP analogs is clearly higher than that of native VIP. Indeed, D3A, F6A, T7A, D8A, Y10A, R12A, and K20A VIP analogs exhibit D values >1.5 Å versus 1.2 Å for native VIP. This indicates that the corresponding analogs have a significant change of conformation as compared with native VIP (see below). In contrast, other mutations did not affect this dispersion with D values <1.2 Å, e.g. H1A, S2A, A4G, K15A, Q16A, M17A, A18G, V19A, Y22A, L23A, S25A, I26A, and N28A VIP analogs. These latter analogs have predicted structures which are very similar to that of native VIP (for example, see the L23A analog in Fig. 3). An intermediate behavior was observed for the V5A, T11A, R14A, K21A, N24A, and L27A VIP analogs which exhibited D values between 1.2 and 1.5 Å. For these analogs, we noticed that the dispersion is mostly due to one of the eight conformers corresponding to one peculiar µ value as shown, for example, for the V5A VIP analog (Fig. 3). In this particular example, when µ = 28 the predicted structure of the V5A analog does not fit as well in the N-terminal domain as for the other µ values. We assumed that a single deviation is not significant and we concluded that an intermediate D value is not representative of a structural change in the analog as compared with native VIP. For analogs with D values >1.5 Å, the high dispersions are clearly due to differences of several, if not all, conformers. The D3A, F6A, T7A, D8A, and Y10A mutants exhibit changes in their N-terminal domain as compared with native VIP (Fig. 3). The R12A mutant exhibits a global alteration of the predicted structure due to partial disruption of the central helical domain (Fig. 3). Finally, the K20A mutant exhibits a change in its C-terminal domain as compared with native VIP (Fig. 3).

View larger version (21K): [in this window] [in a new window]   Fig. 3.   Molecular modeling of native VIP and some VIP analogs: structural modifications triggered by alanine substitutions. The figure shows superimposition of all subsets of possible conformations. Ribbon representation and superposition of the structures obtained with different µ values are shown. Green is for µ = 5; yellow for µ = 7; dark blue for µ = 9; pink for µ = 11; purple for µ = 12; green-yellow for µ = 16; orange for µ = 20; light blue for µ = 28. See "Experimental Procedures" for details. The lowest energy conformation is shown in purple. Native VIP is essentially a -helice from Val5 to Asn24. The N-terminal domain (His1-Ala4) is shown on the left, and C-terminal domain (Ser25-Asn28) on the right are less structured with several possible low energy conformations.

Structure-Activity Relationship of VIP for Interacting with Human VPAC₁ Receptor-- The 28 residues of native VIP can be classified in three categories with respect to the interaction of VIP with human VPAC₁ receptor: (i) residues which can be substituted into alanine or glycine without significant alteration of their binding affinity or biological potency. This category includes Ser², Ala⁴, Asn⁹, Thr¹¹, Leu¹³, Gln¹⁶, Met¹⁷, Ala¹⁸, Val¹⁹, Tyr²², Asn²⁴, Ser²⁵, Leu²⁷, and Asn²⁸. It is quite interesting to note that, compared with native VIP, mutants at these positions do not exhibit significant change in their predicted structure; (ii) residues whose substitution into alanine is associated with a significant alteration of the binding affinity or biological potency of the corresponding analog as well as a change of the predicted analog structure. This category includes Asp³, Phe⁶, Thr⁷, Asp⁸, Tyr¹⁰, Arg¹², and Lys²⁰. Although we can speculate that the decreased binding affinity of the corresponding mutant is due to its altered structure, the possible direct involvement of the category ii residues in the VIP-VPAC₁ receptor interaction cannot be ruled out; (iii) residues whose substitution into alanine results in a decreased binding affinity or biological potency of the corresponding analog without change of the predicted analog structure. This category includes His¹, Val⁵, Arg¹⁴, Lys¹⁵, Lys²¹, Leu²³, and Ile²⁶. These residues are likely to participate in the direct interaction between VIP and human VPAC₁ receptor. The VIP sequence with the three categories of residues is shown in Fig. 4.

View larger version (41K): [in this window] [in a new window]   Fig. 4.   Amino acid sequence of human VIP showing residues which are critical for binding affinity to VPAC1 (top) or VPAC2 receptor (bottom) and/or for maintaining three-dimensional structure. The figure highlights the data shown in Tables I and II. The sequence of human PACAP-27 is aligned with that of VIP to show the identities (solid bar) and homologies (double arrow).

It is worth noting that all analogs exhibiting alteration of the predicted peptide structure, as compared with native VIP, have a decreased binding affinity or biological potency. In this context, it was interesting to align the amino acid sequences of human VIP and human PACAP-27 (Fig. 4) since both peptides have the same high affinity for human VPAC₁ receptors (38). It clearly appears that residues belonging to categories ii and iii are identical in VIP and PACAP or highly conserved, e.g. exchanges between valine and isoleucine at positions 5 and 26.

Alanine Scanning Analysis of VIP for Interacting with Human VPAC₂ Receptor: Selectivity of VIP Analogs-- The ability of each VIP analog to interact with the recombinant human VPAC₂ receptor expressed in CHO clone 10 cells was analyzed by competitive inhibition of ¹²⁵I-VIP binding to cell membranes. All competitor curves fitted with a monophasic binding profile, consistent with a single binding site (not shown). Fig. 5 shows the competition curves for ¹²⁵I-VIP binding to VPAC₂ receptor for some VIP analogs, e.g. T11A, Y22A, and N28A. Analysis of K_i (Table II) showed that substituting residues at 16 positions of VIP resulted in a >1 log decrease of affinity for VPAC₂ receptors. This includes the same 14 positions as for VPAC₁ receptor and two additional spots at positions 11 and 28. Since the corresponding analogs T11A and N28A VIP do not exhibit changes in their predicted structure as compared with native VIP, it can be suggested that Thr¹¹ and Asn²⁸ participate in the direct interaction between VIP and human VPAC₂ receptor. Further differences between VPAC₁ and VPAC₂ receptors can be noticed since some VIP analogs exhibit a much higher decrease in binding affinity for the VPAC₂ receptor than for the VPAC₁ receptor, i.e. D8A, Y10A, and L23A. All VIP analogs were also tested for their ability to stimulate adenylyl cyclase activity in membranes of CHO clone 10 cells expressing the human VPAC₂ receptor. The dose-response curves for all analogs paralleled that of VIP and 10 µM analogs achieved the same maximal stimulation of enzyme activity as native VIP (not shown). The dose-response curves for T11A, Y22A, and N28A analogs are shown in Fig. 5. In general, there is a good correlation between the EC₅₀ for stimulating enzyme activity and the K_i for inhibiting ¹²⁵I-VIP binding (Fig. 6 and Table II). However, the EC₅₀ of Y22A VIP for stimulating adenylyl cyclase activity is much higher than its K_i for binding to VPAC₂ receptor (Fig. 5 and Table II). This property is specific for VPAC₂ receptor since the EC₅₀ and K_i of Y22A VIP at the VPAC₁ receptor are very similar (Table I).

View larger version (23K): [in this window] [in a new window]   Fig. 5.   Binding assay and adenylyl cyclase assay of VIP analogs in CHO cell clone stably expressing human VPAC2 receptor. A, dose effects of native VIP and some VIP analogs for inhibition of 125I-VIP binding to membranes from CHO cells expressing VPAC2 receptor. B, stimulation of adenylyl cyclase activity by native VIP and some VIP analogs in membranes from CHO cells expressing VPAC2 receptor. The data are expressed as the percentage of maximal stimulation above basal obtained with 10 µM native VIP. All data are mean ± S.E. of at least three experiments performed in duplicate. , native VIP; , T11A mutant; , Y22A mutant; , N28A mutant. Note the shift for the Y22A mutant when comparing the binding assay (A) and the adenylyl cyclase assay (B).

View larger version (19K): [in this window] [in a new window]   Fig. 6.   Correlation between the EC50 values and Ki values determined for native VIP and all single-substituted analogs of VIP in the adenylyl cyclase assay and binding assay in CHO cell clone expressing human VPAC2 receptor. All analogs are numbered according to the position of the amino acid that was substituted. Wt, native VIP. See Table II and the legend to Fig. 5 for details. Note the Y22A mutant for which the EC50 for stimulating adenylyl cyclase activity is much higher than its Ki for binding.

From the above described data, it appears that several analogs discriminate between VPAC₁ and VPAC₂ receptors. They exhibit a much higher binding affinity and/or biological potency for VPAC₁ than for VPAC₂ receptors. The reverse is not true since none of the analogs had a higher affinity for VPAC₂ than for VPAC₁. In this context, we synthesized new analogs with the aim to develop more selective human VPAC₁ receptor agonists. For that purpose, we first combined mutations at positions 11 and 28 which individually resulted in a decreased affinity for VPAC₂ receptor without any change in the affinity for VPAC₁ receptor. The [Ala^11,28]VIP analog clearly discriminated between VPAC₁ and VPAC₂ receptors (Tables I and II). Indeed, this VIP analog had the same affinity as native VIP for VPAC₁ receptor whereas it displayed a 44-fold lower affinity than native VIP for VPAC₂ receptor (Fig. 7). Similar data were obtained in the adenylyl cyclase assay (Fig. 7). The [Ala^11,28]VIP analog was as potent as native VIP for stimulating enzyme activity via VPAC₁ receptor whereas it was 21-fold less potent than VIP for stimulating enzyme activity via VPAC₂ receptor (Fig. 7). Since we also noted that position 22 in VIP was important for discriminating between VPAC₁ and VPAC₂ receptors with respect to adenylyl cyclase activation (see above), we further synthesized a VIP analog which combines 3 mutations at positions 11, 22, and 28. Quite interestingly, the [Ala^11,22,28]VIP analog was highly selective for VPAC₁ receptor. Indeed, its binding affinity for VPAC₂ receptor and subsequent biological potency in stimulating adenylyl cyclase were decreased by >1,000-fold as compared with native VIP (Fig. 7). In sharp contrast, this analog with a triple mutation retained a binding affinity for VPAC₁ receptor closed to that of native VIP (Fig. 7). The same held true for biological response, the triple mutant and native VIP having identical potencies for stimulating adenylyl cyclase activity in VPAC₁ receptor-transfected cells (Fig. 7). It is interesting to note that 10⁸ M of this analog which triggers maximal response at VPAC₁ receptor is inactive at VPAC₂ receptor (Fig. 7). This makes [Ala^11,22,28]VIP the most highly selective human VPAC₁ receptor agonist derived from VIP ever described.

View larger version (24K): [in this window] [in a new window]   Fig. 7.   VIP analogs discriminating between human VPAC1 and VPAC2 receptors. Top, dose effects of native VIP and two VIP analogs for inhibition of [125I]VIP binding to membranes from CHO cells expressing VPAC1 (A) or VPAC2 (B) receptor. Bottom, stimulation of adenylyl cyclase activity by native VIP and two VIP analogs in membranes from CHO cells expressing VPAC1 (C) or VPAC2 (D) receptor. The data are expressed as the percentage of maximal stimulation above basal obtained with 10 µM native VIP. All data are mean ± S.E. of at least three experiments performed in duplicate. , native VIP; , [Ala11,28]VIP analog; , [Ala11,22,28]VIP analog.

	DISCUSSION

TOP ABSTRACT INTRODUCTION EXPERIMENTAL PROCEDURES RESULTS DISCUSSION REFERENCES

Although several studies have long been reported regarding the pharmacology of analogs of VIP and related peptides (40-43), no systematic evaluation of the functional and structural role of every residues of VIP has been yet performed since the cloning of the two VIP receptor subtypes VPAC₁ and VPAC₂ (44). Moreover most previous studies were performed using VIP receptors from rodents although it is clear that there are important differences between species in the pharmacology of VIP receptors (8, 15). In the present study, we have synthesized a total of 28 single alanine or glycine mutants of VIP and have analyzed both their predicted three-dimensional structure by molecular modeling and their biological properties by binding assay and adenylyl cyclase assay in CHO cell clones expressing recombinant human VPAC₁ or VPAC₂ receptor. Our data provides new information: (i) they delineate key amino acid residues of VIP which play a role in biological activity for interacting with human VIP receptor subtypes and/or contribute to the three-dimensional structure of the peptide; (ii) they point to the role of specific amino acid residues in VIP for discriminating VPAC₁ and VPAC₂ receptors; (iii) they provide a rationale for developing a new highly specific agonist of the human VPAC₁ receptor by combining three single mutations at key positions of VIP. This new information represents a significant advance over the current knowledge on this widespread neuropeptide.

Analysis of the interaction of the 28 single alanine or glycine mutants of VIP indicates that 14 residues out of 28 in native VIP could not be changed without a significant decrease in affinity for binding to human VPAC₁ receptor. These important residues are distributed along the peptide chain from position 1 to 26. This distribution is in good agreement with previous observations indicating that VIP fragments obtained by deletion of the N-terminal, central, or C-terminal domains of VIP retained very low affinity, if any, for VIP receptors (21). In a previous Ala scan of the VIP derivative Ro 23-7059 the binding activity of analogs had been tested in guinea pig or human lung membranes and emphasized the importance of only 6 positions including Asp³, Phe⁶, Thr⁷, Tyr¹⁰, Tyr²², and Leu²³ (22). These data are somewhat different from ours since we found many other important positions and identified Tyr²² as an unimportant position, at least for the interaction of VIP with VPAC₁ receptor. The reasons for these discrepancies probably include several factors: (i) a VIP derivative was used, i.e. Ac-[Lys¹²,Nle¹⁷,Val²⁶,Thr²⁸]VIP (22); (ii) the authors used lung membranes as a source of VIP receptors whereas we used recombinant receptor subtypes in this study. In this context, it has been shown since the initial study of O'Donnell et al. (22) that lung is endowed with both VPAC₁ and VPAC₂ receptors (45); (iii) the authors utilized guinea pig tissues in their studies (22). It is known that guinea pig peptide hormones or neuropeptides have divergent amino acid sequences as compared with most mammals (46). In the case of VIP, the sequence is unchanged in all mammals with the notable exception of guinea pig whose VIP has four changes (47); (iv) the bioassay used in their study (22) consisting of relaxation of guinea pig trachea appears to be less sensitive than the adenylyl cyclase assay used in our study.

In the present study, we identified His¹, Asp³, Val⁵, Phe⁶, Thr⁷, Asp⁸, Tyr¹⁰, Arg¹², Arg¹⁴, Lys¹⁵, Lys²⁰, Lys²¹, Leu²³, and Ile²⁶ as being important for binding to human VPAC₁ receptor and subsequent stimulation of adenylyl cyclase activity (see Table I). It is quite interesting to note that all these amino acids are conserved in PACAP (see Fig. 4) which has the same affinity as VIP for human VPAC₁ receptor (39). Indeed, the amino acids are either identical or are highly homologous with Val⁵ and Ile²⁶ in VIP being replaced by Ile⁵ and Val²⁶ in PACAP, respectively (see Fig. 4). This conservation strongly validates the present data. Futher analysis of the 14 important amino acid residues in VIP may be done by comparing the amino acid sequences of VIP and VIP-related peptides which do (PACAP, secretin, peptide histidine methionine amide, and GRF) or do not (glucagon, glucagon-like peptides and GIP) behave as VIP agonists at the VPAC₁ receptor (12, 39, 44). This comparison supports that His¹, Asp³, Phe⁶, Arg¹², Lys²⁰, Lys²¹, Leu²³, and Ile²⁶ which emerged from the present Ala scan as important residues in VIP for interacting with human VPAC₁ receptor, are strictly conserved or replaced by homologous residues in VIP-related peptides which behave as VIP agonists. This suggests that these residues may play a similar role in allowing the binding of VIP, PACAP, peptide histidine methionine amide, GRF, and secretin to human VPAC₁ receptor but altogether are not sufficient to ensure high affinity binding. To further substantiate this issue, it is worth pointing out that (i) we could not find well conserved amino acid residues in VIP-related peptides which do not play a role in the activity of VIP; (ii) conversely, we did identify nonconserved residues which are important for VIP activity at the human VPAC₁ receptor.

At this stage of the Ala scan study of VIP and more generally of other similar studies with any bioactive peptide, the question arises of whether one given residue is important because it directly interacts with the receptor and/or it plays a role in maintaining the three-dimensional structure of the peptide, an issue which is not often addressed. In the present study, the conformations of native VIP and VIP analogs were calculated by molecular modeling using an ab initio approach (31). When this approach was applied to native VIP, the predicted three-dimensional structure of VIP was consistent with a central -helical domain from Val⁵ to Asn²⁴ and conformational variability essentially located at the N- and C-terminal domains (see Fig. 3). This is in good agreement with previous experimental studies of the VIP structure in methanol/water solutions using nuclear magnetic resonance and circular dichroism which also indicated that the structure of VIP is mostly helical (48) with the existence of a central well defined -helix, the remaining residues being not ordered (49). Similarly, PACAP38 was shown to have a disordered N-terminal domain followed by a central -helical structure (50). Interestingly, helodermin, a peptide isolated from the lizzard Gila Monster (51), which is highly homologous to VIP (51) and behaves as a VIP agonist (52), has also been shown to exhibit a central -helice with random coiled N and C terminus (53). With respect to three-dimensional structure, VIP analogs in which amino acid substitution into alanine was associated with alteration of binding affinity and biological potency could be classified into two categories: (i) analogs which exhibited significant change in predicted peptide structure as compared with native VIP. This includes D3A, F6A, T7A, D8A, Y10A, R12A, and K20A mutants. For these analogs, it can be suggested that altered activity is related to altered structure, although we cannot strictly rule out the possibility that the corresponding residues may directly participate in the interaction of VIP with the VPAC₁ receptor; (ii) analogs which exhibited no change in predicted peptide structure as compared with native VIP. This includes H1A, V5A, R14A, K15A, K21A, L23A, and I26A mutants. It is tempting to speculate that the corresponding amino acid residues in native VIP (see Fig. 8) participate in the interaction of the neuropeptide with VPAC₁ receptor. This is in line with other observations: (i) His¹ in VIP (22) as well as His¹ in some VIP-related peptides (54-56) is known to play an important role in biological activity; (ii) the fact that half of the residues falling into the category of candidates for direct interaction with VPAC₁ receptors are basic residues including Arg¹⁴, Lys¹⁵, and Lys²¹, is probably significant in view of the fact that site-directed mutagenesis of the human VPAC₁ receptor previously identified several important acidic residues in the receptor for VIP binding, including Glu³⁶ (16), Asp⁶⁸ (12), and Asp¹⁹⁶ (57). This supports the idea that the very basic VIP with an isoelectric point >11 (6) may use some of its basic residues for an electrostatic interaction with acidic residues of the receptor; (iii) finally, the remaining important residues of VIP including Val⁵, Leu²³, and Ile²⁶ are highly hydrophobic in consonance with the old idea that VIP, or at least a binding region of VIP, has an hydrophobic environment within the receptor (58). Recent observations support the idea that important residues in the VPAC₁ receptor for binding VIP are indeed hydrophobic including Trp⁶⁷, Trp⁷³, and Trp¹¹⁰ (12).² Finally, it is worth pointing out that all analogs which have the same affinity as native VIP for VPAC₁ receptor exhibit no alteration of predicted three-dimensional structure, further validating our molecular modeling approach.

View larger version (63K): [in this window] [in a new window]   Fig. 8.   Ribbon representation of the VIP molecule highlighting residues which are likely to participate in the direct interaction of VIP with the human VPAC1 (A) and VPAC2 (B) receptors. These residues shown in red are His1, Val5, Arg14, Lys15, Lys21, Leu23, and Ile26 for interaction with VPAC1 and VPAC2 receptors. Three additional residues shown in blue are important only for interaction with VPAC2 receptor, e.g. Thr11, Tyr22, and Asn28.

The data obtained from the analysis of the interaction of the 28 single alanine or glycine mutants of VIP with human VPAC₂ receptor are very similar to those obtained with VPAC₁ receptors with some specific differences. The great similarities indicate that VIP binding to VPAC₁ and VPAC₂ receptor is mostly ensured by a common set of crucial amino acids along the VIP peptide sequence. However, two important differences can be noted: (i) some VIP analogs exhibit a much lower binding affinity for the VPAC₂ than for the VPAC₁ including the D8A, Y10A, and L23A mutants. Two other analogs have no alteration of binding to VPAC₁ receptor whereas they exhibit a significant decrease of affinity for VPAC₂ receptor, i.e. T11A and N28A analogs (see Fig. 5 and Table II). When these two modifications are combined in one single analog, the resulting peptide behaves as a VIP agonist which clearly discriminates between human VPAC₁ and VPAC₂ receptors (see Fig. 7). (ii) The Y22A VIP analog exhibits an increased EC₅₀ for stimulating adenylyl cyclase activity through interaction with the human VPAC₂ receptor whereas no change was observed at the VPAC₁ receptor (see Tables I and II). Altogether, these data indicate that single amino acid substitutions in VIP either alter similarly peptide activity at the two VIP receptor subtypes or evoke a higher decrease of activity at the VPAC₂ receptor than at the VPAC₁ receptor. In other words, no selective VPAC₂ receptor ligand was observed in the series of analogs obtained by the Ala scan. Conversely, we took advantage of the selectivity observed toward the VPAC₁ receptor (see above) to develop a VIP analog which combined 3 mutations at positions 11, 22, and 28 which were shown to be individually important for discriminating human VPAC₁ and VPAC₂ receptors. The resulting [Ala^11,22,28]VIP analog is an agonist which exhibits an exquisite high selectivity for VPAC₁ receptor with EC₅₀ <1 nM and >1,000 nM in stimulating adenylyl cyclase activity through interaction with human VPAC₁ and VPAC₂ receptors, respectively (see Fig. 7). Two VPAC₁ receptor selective agonists derived from secretin and GRF were previously described (59, 60). Chicken [Arg¹⁶]secretin, although it discriminates between rat VPAC₁ and VPAC₂ receptors, has a high affinity for rat secretin receptor (59). Moreover, it has a rather weak affinity for the human VPAC₁ receptor (60), in line with the important differences between species in the pharmacology of VPAC₁ receptors (15, 44). The chimeric, substituted peptide [Lys¹⁵,Arg¹⁶,Leu²⁷]VIP (1-7)/GRF (8-27) has an important selectivity for VPAC₁ receptor but its possible interaction with GRF receptors was not directly evaluated (61). In this context, the highly selective human VPAC₁ receptor agonist developed in this study by rationale combination of three mutations in VIP itself constitutes henceforth the most operative pharmacological tool derived from VIP for characterizing VPAC₁ receptor-mediated events.

In conclusion, the present activity data on the series of VIP analogs in combination with molecular modeling provide the first characterization of the functional properties of the VIP molecule for interaction with the two human VIP receptor subtypes. Key residues for interaction of VIP with human VPAC₁ or VPAC₂ receptor are highlighted in the ribbon representation of the VIP molecule (Fig. 8). In view of considerable interest of this widespread neuropeptide in health and diseases, this work already provides useful pharmacological tools and paves the way for further development in the design of new analogs and mimetic versions of VIP.

	FOOTNOTES

^* The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

To whom correspondence should be addressed: INSERM U410, Faculté de Médecine Xavier Bichat, 75018 Paris, France. Tel.: 33-01-44-85-61-35; Fax: 33-01-44-85-61-24; E-mail: laburthe@bichat.inserm.fr.

Published, JBC Papers in Press, May 5, 2000, DOI 10.1074/jbc.M002325200

² P. Nicole, L. Lins, C. Rouyer-Fessard, C. Drouot, P. Fulcrand, A. Thomas, A. Couvineau, J. Martinez, R. Brasseur, and M. Laburthe, unpublished data.

	ABBREVIATIONS

The abbreviations used are: VIP, vasoactive intestinal peptide; PACAP, pituitary adenylate cyclase-activating peptide; GRF, growth hormone-releasing factor; VPAC receptor (for official nomenclature see Ref. 8), CHO, Chinese hamster ovary cells; Fmoc, N-(9-fluorenyl)methoxycarbonyl.

	REFERENCES

TOP ABSTRACT INTRODUCTION EXPERIMENTAL PROCEDURES RESULTS DISCUSSION REFERENCES