Cécile Echalierab,
Aleksandra Kalistratovaac,
Jérémie Ciccioneab,
Aurélien Lebruna,
Baptiste Legranda,
Emilia Naydenovac,
Didier Gagnea,
Jean-Alain Fehrentza,
Jacky Mariea,
Muriel Amblarda,
Ahmad Mehdib,
Jean Martineza and
Gilles Subra*a
aInstitut des Biomolécules Max Mousseron (IBMM), UMR5247 Université de Montpellier, CNRS, ENSCM, Faculté de Pharmacie, 15 avenue Charles Flahault, 34093 Montpellier, France. E-mail: gilles.subra@univ-montp1.fr; Tel: +33 4 11 75 96 06
bInstitut Charles Gerhardt (ICGM), UMR5253 Université de Montpellier, CNRS, ENSCM, Place Eugène Bataillon, 34095 Montpellier, France. Tel: +33 4 67 14 38 32
cUniversity of Sofia, bulevard “Sveti Kliment Ohridski” 8, Sofia, Bulgaria
First published on 24th March 2016
We developed a simple and straightforward way to dimerize unprotected peptide sequences that relies on a chemoselective condensation of hybrid peptides bearing a hydroxydimethylsilyl group at a chosen position (either C-ter, N-ter or side-chain linked) to generate siloxane bonds upon freeze-drying. Interestingly, the siloxane bond sensitivity to hydrolysis is strongly pH-dependent. Thus, we investigated the stability of siloxane dimers in different experimental conditions. For that purpose, 29Si, 13C and 1H NMR spectra were recorded to accurately quantify the ratio of dimer/monomer. More interestingly, we showed that 1H resonances of the methylene and methyl groups connected to the Si can be used as sensitive probes to monitor siloxane hydrolysis and to determine the half-lives of the dimers. Importantly, we showed that the dimers were rather stable at pH 7.4 (t1/2 ≈ 400 h) and we applied the dimerization strategy to bioactive sequences. Once optimized, three dimers of the growth hormone releasing hexapeptide (GHRP-6) were prepared. Interestingly, their pharmacological evaluation revealed that the activity of the dimeric ligands could be switched from agonist to inverse agonist depending on the position of dimerization.
In this study, we focused our efforts on an alternate dimerization process based on hydroxydimethylsilane peptides. To our knowledge, most of the existing strategies require orthogonal protecting groups and multistep procedures to obtain homodimeric peptides at the exception of dimerization of unprotected peptide sequence containing cysteine residue through the formation of a disulfide bond.11 Advantageously, the use of dimethylhydroxysilane peptides allows the dimerization of any peptides, even those containing a Cys residue and greatly simplifies and speeds up the process. Indeed, dimerization occurs chemoselectively in the presence of unprotected side chains.
Moreover, it is of particular interest that hybrid dimethylhydroxysilane monomers are prepared by classical peptide synthesis and can be purified before dimerization. We recently applied this strategy to the synthesis of linear and branched peptide-polymers and tripeptide dimers.12 In this work, we aimed at demonstrating that the strategy is not only applicable to the synthesis of N-ter dimers but also to C-ter and internal bioactive peptide sequence dimerization. Interestingly we proved here that siloxane formation simply occurred during lyophilization after preparative RP-HPLC of the hybrid monomer. The 29Si NMR signals were different in siloxane bonds and silanols but we demonstrated that, 1H NMR could be used advantageously to measure the SiOH/SiOSi ratio. Thanks to quick and sensitive 1H NMR analyses, we monitored the stability of siloxane bonds in different solvents but more importantly in water at various pHs. At last, we synthesized three types of dimers based on the sequence of the growth hormone releasing hexapeptide (GHRP-6). The dimerization strategy was applied at three different positions of the peptide sequence: C-terminal, N-terminal and intra-sequence through a modified lysine residue (Scheme 1). Then, the binding and signaling properties of these dimers were evaluated on intact HEK293T cells expressing the growth hormone secretagogue receptor 1a (GHS-R1a) involved in important physiological pathways.13
![]() | ||
| Scheme 1 Synthesis of GHRP-6 peptide homodimers based on dimethylhydroxysilyl peptides. 5: N-terminal dimerization; 7: C-terminal dimerization; 9: internal residue dimerization. | ||
The following abbreviations were used: ICPDMCS – isocyanatopropyldimethylchlorosilane. Other abbreviations used are recommended by the IUPAC-IUB Commission.
120 Hz). For determination of silicon chemical shifts, a DEPT sequence was used to benefit from the polarization transfer of hydrogens near the silicon atom (256 scans, spectral width of 11
900 Hz and recovery delay of 2 s). 1H NMR spectra were typically recorded using 16 scans with a spectral width of 6000 Hz and a recovery delay of 2 s. Spectra were processed and visualized either with Topspin 3.2 (Bruker Biospin) on a Linux station or with MestReNova 6.0.2 (Mestrelab research).
000 cells per well. IP1 production was measured using the IP-One HTRF kit (Cisbio Bioassays) as previously described.16 Briefly, cells were stimulated for 30 min at 37 °C with the ligand to be tested in 70 μl of IP1 stimulation buffer. An anti-IP1 antibody labelled with Lumi4-Tb (15 μl) and an IP1-d2 derivative (15 μl), were added on cells. The medium was incubated for 1 h at room temperature. Signals at 665 nm and 620 nm were detected using a PHERAstar (BMG LABTECH) fluorescence reader. Values are expressed as ΔF. ΔF corresponded to: (ratio 665 nm/620 nm of the assay − ratio 665 nm/620 nm of the negative control)/ratio 665 nm/620 nm of the negative control. The negative control corresponded to the Lumi4-Tb blank and was used as an internal assay control. Inositol phosphate accumulation was expressed as the percentage of the maximal ghrelin response using the formula: (ΔF mock cells − ΔF receptor transfected cells)/(ΔF mock cells − ΔF maximal ghrelin stimulation for receptor transfected cells). EC50 values were determined from dose–response curves using GraphPad Prism software (GraphPad Software, Inc., San Diego, CA).
![]() | ||
| Fig. 1 Model compounds for NMR determination of 1H, 13C and 29Si shifts. Representative atoms of interest for the identification of both dimer and monomer are marked in bold. | ||
Compound 1a was obtained from commercially available 1,3-bis(3-aminopropyl)tetramethyl disiloxane poured into hydrogen chloride ether solution. On the contrary of its neutral counterpart, 1a was soluble in water enabling the NMR analysis in aqueous solutions. Compound 1a was solubilized either in DMSO-d6 or D2O and analyzed by 1H, 13C and 29Si NMR. Results are reported in Table 1 (first column, dimer 1a). As expected, when compound 1a was left in D2O for a prolonged time, monomer 2a appeared upon hydrolysis. This latter was analyzed by NMR (Table 1). Hydrochloric acid was added to dimer 1a in DMSO-d6 to force hydrolysis and subsequently identify the chemical shifts of monomer 2a (Table 1). Since all the peptides were purified in trifluoroacetic acid-containing solvents and therefore recovered as trifluoroacetate salts, dimer 1b was also prepared. It was prepared from 1,3-bis(3-aminopropyl)tetramethyl disiloxane poured into trifluoroacetic acid ether solution and solubilized in DMSO-d6 or D2O. Partial hydrolysis occurred with time and the mixture of monomer 2b and dimer 1b was analyzed by 1H, 13C and 29Si NMR. Chemical shifts of monomer 2b and dimer 1b are reported in Table 1.
| Solvent | Nucleus | Dimer 1a | Monomer 2a | Dimer 1b | Monomer 2b | ||||
|---|---|---|---|---|---|---|---|---|---|
| –CH2Si | CH3Si | –CH2Si | CH3Si | –CH2Si | CH3Si | –CH2Si | CH3Si | ||
| DMSO-d6 | 1H | 0.53 | 0.06 | 0.48 | 0.00 | 0.52 | 0.06 | 0.48 | 0.01 |
| 13C | 14.6 | 0.2 | 14.5 | 0.0 | 14.5 | 0.1 | 14.4 | −0.1 | |
| 29Si | 8.0 | 10.9 | 8.0 | 10.8 | |||||
| D2O | 1H | 0.71 | 0.22 | 0.74 | 0.23 | 0.70 | 0.21 | 0.72 | 0.22 |
| 13C | 13.6 | −1.2 | 12.9 | −2.3 | 13.6 | −1.2 | 12.9 | −2.3 | |
| 29Si | 11.4 | 18.2 | 11.4 | 18.1 | |||||
Analyses of monomer/dimer mixtures by 29Si NMR showed two signals at 8.0 and 10.8 (±0.1) ppm in DMSO-d6 and 11.4 and 18.1 (±0.1) ppm in D2O that corresponded to the silicon signal involved in siloxane Si–O–Si and hydroxysilane Si–OH bonds respectively (Fig. 2). Whatever the conditions of analysis, 29Si NMR unambiguously allowed the identification of hybrid monomeric and dimeric species. The 13C NMR spectra also allowed us to readily discriminate the monomer from the dimer signals looking at the methyl and the methylene resonances connected to the Si. 29Si and 13C NMR have the advantage to display a large spectral width and only few 29Si and 13C signals are present in our systems. Nevertheless, due to the low natural abundance (4.7% and 1.1%, respectively) and low NMR sensitivity of 29Si and 13C, long acquisition times are required even at 0.1 M concentration. In this context, we investigated the characterization of all derivatives using 1H NMR spectra which are far more sensitive than 29Si and 13C NMR despite a narrowest spectral width and the larger number of 1H resonances. In this context, we switched our attention to the NMR signals of the hydrogen on carbon atoms in alpha position to the silicon atoms (Fig. 1). Interestingly, 1H chemical shifts of –CH2Si and CH3Si of the monomer and dimer were localized in a rather clear region of the 1H NMR spectra of the hybrid peptides (≤0.7 ppm) (Table 1). These resonances should not overlapped with those of the methyl groups of the Leu, Ile and Val expected around 0.7–0.8 ppm. In addition, we showed with the model compounds 1a–b and 2a–b that in DMSO-d6 1H NMR chemical shifts of the monomer and dimer were clearly different whatever the amine counter-ion, thereby hybrid dimers and monomers could be unambiguously identified by a simple 1H NMR. In D2O, 1H NMR methyl signals of monomer and dimer as hydrochloride salts were partially overlayed. However, the counter-ion switch from hydrochloride to trifluoroacetate induced a splitting of monomer and dimer signals (ESI Fig. S22†). As a consequence, TFA salts should be preferred for a better signal separation in D2O. In contrast to the 29Si and the 13C resonances, the 1H chemical shifts values of the –CH2Si (∼0.7 ppm) and CH3Si (∼0.2 ppm) were rather close in both species and sensitive to the pH which could lead to confusion between monomer and dimer during assignment. When both monomer and dimer were present, signal assignment was easy as NMR peaks always came out in the same order: in D2O, dimer protons were typically more shielded than monomer hydrogens. In case of any doubt about a sample containing only one species, 29Si NMR analyses had to be undertaken. In order to significantly decrease the experiment time, we used a DEPT24 sequence instead of the typical 29Si quantitative sequence (ESI Fig. S23†).
![]() | ||
| Fig. 2 29Si, 13C and 1H NMR spectra of the dimer 1b/monomer 2b mixture recorded in DMSO-d6. Chemical shifts of relevant signals (protons, carbons, and silicon atoms of –CH2Si and CH3Si) are reported in Table 1. Integrations suggested that the sample contained 83% of dimer 1b and 17% of monomer 2b. | ||
When both monomer and dimer are present in a sample, determination of their proportion is valuable. For that, we used 13C and 29Si quantitative sequences which require long acquisition times. Interestingly, the 29Si DEPT sequence was not expected to be quantitative; however we noticed that when applied to our hybrid systems, signal integrations reflected the proportion of monomer and dimer in the sample (Fig. 2, 29Si spectra). To go further, signal integrations of 1H spectra were compared with 29Si NMR spectra. Interestingly, integration of NMR signals related to dimer or monomer showed that the monomer/dimer ratio was unchanged (Fig. 2). Thus, 1H could be used to determine the proportion of monomer and dimer in a sample. This was particularly advantageous for kinetic studies where high rate sampling was required to monitor the hydrolysis of siloxane bond. Indeed, long acquisition times provided an average ratio that did not make sense when sample evolution was fast. Furthermore, due to the high sensitivity of 1H NMR, the spectra always presented a higher ratio signal to noise by comparison to those of 29Si, which led to a lower detection limit and more accurate integrations.
We showed that simple 1H NMR can be routinely used to analyze hybrid peptides and quantify hydroxydimethylsilane/dimethylsiloxane ratio. The presence of TFA salts seemed to be preferable in D2O for a better signal separation. However, when possible, DMSO-d6 should be preferred as solvent since chemical shifts in DMSO-d6 were less sensitive to salts and pH, resulting in an easier signal assignment. The analysis of hybrid peptide ligands (4–9) used for biological assays was performed by 1H NMR in DMSO-d6.
As expected, the closer the pH was to the neutrality, the slower was the dimer hydrolysis. When the media is moved to acidic and basic pH (4 and 10, respectively), the hydrolysis rate increased, particularly in basic medium. The half-life of dimer 1a was 203 min at pH 4 and 54 min at pH 10. Most of the time, bioactive peptides are administered intravenously. As a consequence, a good stability of hybrid peptide dimers at blood pH was mandatory. Importantly, at pH 7.4, the half-life of dimer 1a was 439 h, with less than 10% of hydrolysis within 2 days. This stability allows envisioning biological applications for hybrid dimers.
In order to mimic the RP-preparative HPLC experimental conditions, the behavior of dimer 1a in more acidic conditions (pH 2) was studied. As expected a few minutes after solubilization of dimer 1a in D2O/0.1% TFA, only monomer was detected (Fig. S33†). Thus, when RP-preparative HPLC is required, hybrid compounds need to be purified as monomer species prior to dimerization.
We evaluated the impact of GHRP-6 dimerization in terms of binding affinity and more importantly, in terms of bioactivity. Indeed, we hypothesized that, depending on the position of the bridge between monomers, the agonism behavior of GHRP-6 sequence could be impacted. To answer these questions, three hybrid derivatives of GHRP-6 were prepared (Fig. 4).
Peptide sequences of the three hybrid derivatives of GHRP-6 were prepared on solid support by conventional Fmoc/tBu strategy. For N-terminus dimerization, a β-alanine spacer was added after the last amino acid and the peptidyl resin was reacted with the 3-isocyanatopropyldimethylchlorosilane (ICPDMCS) to yield the hybrid monomer 4 after cleavage and side-chain deprotection. In contrast, the C-terminus modification was performed through introduction of an ethylene diamine. For that purpose, ethane-1,2-diamine was anchored on 2-chloro chlorotrityl PS resin and the sequence was elongated to the N-terminus via conventional Fmoc/tBu SPPS. The N-Boc, fully protected peptide was cleaved from the solid support by mild acidic treatment (HFIP/DCM 1/4 v/v, 30 min) to maintain the side chain protecting groups. The protected peptide was then reacted with ICPDMCS in solution, precipitated in diethyl ether and treated by TFA to remove the side-chain protecting groups to yield the hybrid peptide 6. To conserve the N and C-termini unchanged, dimerization within the peptide sequence was performed, which might induce a non-negligible steric hindrance in the center of the sequence. Alanine residue in position 3 of GHRP-6 was swapped with a lysine, which served as anchoring point for ICPDMS. Selective removal of Alloc protecting group on the Lys3 side-chain was performed with phenylsilane and Pd0(PPh3)4 on the sequence still anchored on solid support. The ε-amino group was reacted with ICPDMCS before cleavage of the hybrid peptide from Rink amide resin, yielding hybrid peptide 8.
In all cases, the last step of the synthesis of dimethylhydroxysilyl hybrid peptides was a TFA treatment, which removed the acid sensitive protecting groups of the side chains and cleaved the linker when the synthesis was entirely performed on solid support (compounds 4 and 8). TFA was removed under reduced pressure and the hybrid peptide was recovered as its monomeric form by precipitation in diethyl ether. Hybrid monomers 4, 6 and 8 were purified by preparative RP HPLC. We then applied the dimerization procedure by solubilization of the hybrid peptide monomers in aqueous solution at pH 7.4 at high concentration (>0.2 M) as previously described for polymers obtained from bis-silylated precursors.12 In most cases, a precipitate appeared at this stage and condensation products were recovered by centrifugation. This procedure proved to be efficient for dimer 5, which precipitated and was recovered by centrifugation with good yield (90%). Only dimer 5 was observed by 13C NMR. However, a more sensitive 1H NMR analysis indicated the presence of 3% of monomer 4. Unfortunately, no precipitation occurred for dimers 7 and 9, probably because of the presence of the free N-terminal amine, which enhanced water solubility.
First of all, all dimers were able to bind to ghrelin receptor in the 100 nM range. However, compared to GHRP-6, a 5 to 10-fold decrease of binding was observed for compounds 5 and 7. Internal dimerization (compound 9) induced a 15-fold loss of affinity. N-Terminus dimerization (compound 5) had the less impact on the binding affinity. These results could be related to the steric hindrance of the dimeric molecule to approach GHS-R1a binding sites, the transmembrane region 3 of the GHS-R1a in which peptide ligands and small molecules can bind.23 However, binding behavior cannot be associated to bioactivity by a direct relationship. We thus investigated the biological activity of the ligands by measuring their ability to promote inositol phosphate (IP) production in HEK293T cells expressing the GHS-R1a. The production of IP, reported as % of ghrelin-induced maximal responses, are reported in Fig. 6 and summarized in Table 3. The peptide [D-Arg1, D-Phe5, D-Trp7,9, Leu11]-substance P (noted SPA for substance P analogue)24 is a non selective ligand of GHS-R1a, which behaves as an inverse agonist. It was also included in the same experiments.
| Compound | N-ter dimer 5 JMV 6187 | C-ter dimer 7 JMV 6186 | Internal dimer 9 JMV 618 | SPA | Compound 3 GHRP-6 | Ghrelin |
|---|---|---|---|---|---|---|
| a EC50 and Emax were obtained from dose–response curve analysis with GraphPad Prism v5 software. For agonist compounds, Emax is expressed as the % of maximal ghrelin. For inverse agonists, Emax is expressed as the maximal basal inhibition (100% inhibition corresponding to the basal of mock-transfected HEK293T cells). Values are mean ± S.E.M of three independent experiments. | ||||||
| EC50a (nM) | 20.4 ± 6.8 | 2.57 ± 0.35 | 160 ± 44 | 27 ± 2 | 1.18 ± 0.82 | 0.85 ± 0.29 |
| Emax (%) | 86.5 ± 2.7 | 97.8 ± 1.7 | 38 ± 3 | 77.5 ± 7.5 | 99.5 ± 2.2 | 100 |
| Behaviour | Partial agonist | Full agonist | Inverse agonist | Inverse agonist | Full agonist | Full agonist |
Interestingly very different behaviors were obtained for the three dimers. N-Terminal dimer 5 which was the best binder within the dimer series, proved to be a partial agonist with an EC50 of 20 nM. It induced 86.5 ± 2.7% of the maximal effect of the ghrelin. In contrast, the C-terminal dimer 7, which had a three-fold lower affinity for GHS-R1a, was a full agonist with an EC50 of 2.57 nM, close to that of GHRP-6 (1.2 nM). The most striking result was obtained with dimer 9, which behaved as an inverse agonist with an efficacy of 38 ± 3% (Table 3). Interestingly, it was reported that the peptide D-Lys3-GHRP-6 displayed some inverse agonist activity in the sea bream25 but behaved as an antagonist in mammals.26,27 In the case of dimer 9, a L-lysine was introduced as an anchoring point of the dimethylhydroxysilane moiety. Therefore, it could not be excluded that the inverse agonist effect was not induced by dimerization but simply by the Lys/Ala substitution, even though the stereochemistry was not the same (D-Lys in the reported monomeric inverse agonist and L-Lys in the dimer 9). Thus, to know the significance of either the dimerization or the modification of the linear sequence itself in inducing the inverse agonist activity, two other compounds were prepared (Fig. 7). Compound 10 was obtained in the same way as monomer 8 except isocyanopropyltrimethylsilane was used instead of ICPDMCS for the formation of the urea bond. Lys3-GHRP-6 (compound 11) was also prepared. Bioactivity of 10, 11 and of Lys3 N-ε acetylated analogue 12 was also evaluated on HEK293T cells expressing the GHS-R1.
Interestingly none of them displayed any inverse agonist activity and all behaved as neutral ligands (ESI Fig. S54†). This last result showed that the dimerization in itself induced a complete switch of activity and the position of the dimerization within the active sequence is accordingly of great impact. Further pharmacological studies will be performed to investigate the structure–activity relationships of GHRP-6 analogues.
Footnote |
| † Electronic supplementary information (ESI) available: NMR spectra, kinetic studies, compound characterization, monomeric ligand binding assays and bioactivity. See DOI: 10.1039/c6ra06075g |
| This journal is © The Royal Society of Chemistry 2016 |