Efficacy of inverso isomer of CendR peptide on tumor tissue penetration.

The dense extracellular matrix and high interstitial fluid pressure of tumor tissues prevent the ability of anti-tumor agents to penetrate deep into the tumor parenchyma for treatment effects. C-end rule (CendR) peptides can enhance the permeability of tumor blood vessels and tumor tissues via binding to neuropilin-1 (NRP-1), thus aiding in drug delivery. In this study, we selected one of the CendR peptides (sequence RGERPPR) as the parent l-peptide and substituted d-amino acids for the l-amino acids to synthesize its inverso peptide D(RGERPPR). We investigated the NRP-1 binding activity and tumor-penetrating ability of D(RGERPPR). We found that the binding affinity of D(RGERPPR) with NRP-1 and the cellular uptake was significantly higher than that of RGERPPR. Evans Blue tests revealed that D(RGERPPR) exhibited improved tumor-penetrating ability in C6, U87 and A549 tumor-bearing nude mice. Using nude mice bearing A549 xenograft tumors as a model, we found that the rate of tumor growth in the group co-administered with D(RGERPPR) and gemcitabine (Gem) was significantly lower than the gemcitabine-treated group with a tumor suppression rate (TSR%) of 55.4%. Together, our results demonstrate that D(RGERPPR) is a potential tumor-penetrating peptide.

Introduction

Tumor tissues exhibit specific pathological features, including large amounts of random vasculature with high permeability, the absence of lymphatic drainage, a dense extracellular matrix, high interstitial fluid pressure, and an acidic and anoxic environment in the center of the tumor that is caused by poor elimination of metabolites owing to insufficient blood supply1, 2, 3, 4, 5, 6, 7, 8. Because of the dense extracellular matrix and high interstitial fluid pressure, antitumor agents are unable to reach deep into the tumor parenchyma, leading to retention of drugs around the tumor blood vessels. This non-uniform distribution of drugs is one of the crucial reasons for relapse and metastasis, which is serious obstacle in the treatment of tumors.

Neuropilin-1 (NRP-1), a transmembrane glycoprotein, was first discovered in the nervous system of Xenopus tadpoles by Takagi et al. via hybridoma technology in 1987. As a coreceptor of sema3A and VEGF165, NRP-1 played an important role in tumor occurrence, growth, and metastasis10, 11, 12. NRP-1 was highly expressed on many tumor cells including glioma, non-small cell lung cancer, pancreatic cancer, gastric cancer, colon cancer and breast cancer13, 14, 15. A series of peptides possessing a consensus R/KXXR/K motif with l-amino acids were identified that showed high binding affinity to NRP-1. The C-terminal arginine (or rarely lysine) was important for NRP-1 binding activity. Adding another amino acid or blocking the free carboxyl group of this arginine residue by amidation eliminated the binding activity to NRP-1. This strict requirement for a C-terminal residue was termed as the C-end rule (CendR). By binding NRP-1, the CendR peptides could induce enhanced permeability of tumor blood vessels and tumor tissues, which could enhance the tissue-penetrating ability of free drugs and nanodrug delivery systems.

Despite these initial promising results, these linear l-peptides are susceptible to protease degradation in plasma and readily lose their bioactivity. In our previous study, we selected one of CendR peptides (sequence RGERPPR) as a parent l-peptide and synthesized its retro-inverso isomer, D(RPPREGR)16, 17. The retro-inverso peptide was composed of d-amino acids assembled in the reverse order from that in the parent l-sequence. Compared with RGERPPR, D(RPPREGR) not only showed improved biological stability but also increased binding affinity to NRP-1. We synthesized D(RPPREGR) modified highly branched polyethyleneimine as gene vector, which significantly enhanced the therapeutic effect of genes on orthotopic tumor-bearing animal models.

The inverso peptide, which is composed of d-amino acids assembled in the same order as that in the parent l-sequence, shows a mirror-image relationship with the parent l-peptide. In general, the inverso peptide might not have similar bioactivity as its parent l-peptide. However, CendR peptides are a series of specific peptides with a consensus R/KXXR/K motif and exhibited high binding affinity to NRP-1 and improved tumor-penetrating ability. The inverso isomer of CendR peptides also has a consensus D(R/KXXR/K) motif. Since the C-terminal amino acid of the inverso isomer is arginine or lysine, we speculate that the inverso isomer of CendR peptides also have NRP-1 binding activity and tumor-penetrating ability. To verify our hypothesis, in this paper, we chose RGERPPR as the parent l-peptide and synthesized its inverso isomer, D(RGERPPR). We then investigated the activities of D(RGERPPR), including binding affinity to NRP-1, tumor cell internalization and tumor-penetrating ability.

Methods

Materials

Boc-protected amino acid was purchased from Peptide International (Japan). Hydrogen amine resin and 4-methylbenzene was supplied by Bachem AG (Switzerland). Diisopropylethylamine (DIEA) was obtained from Fluka (USA). O-(Benzotriazol-1-yl)-N,N,N′,N′-tetramethyluronium hexafluorophosphate (HBTU) was from American Bioanalytical (USA). iRGD (CRGDKGPDC) was purchased from ChinaPeptides (Shanghai, China). Fluorescein-5-maleimide (Mal-FITC) was purchased from Fanbo Biochemicals (Beijing, China). Phosphate-buffered saline (PBS), plastic cell culture dishes and plates were supplied by Corning Incorporated (Corning, USA). Dulbecco׳s modified Eagle׳s medium (DMEM), fetal bovine serum (FBS), trypsin-EDTA (0.25%) and penicillin–streptomycin were from Gibco (Life Technology, USA). Evans Blue was purchased from Ourchem (Shanghai, China). Gemcitabine hydrochloride and carboplatin were obtained from Meilun Biological Technology Co. Ltd. (Dalian, China). Thiazolyl blue tetrazolium bromide (MTT) and all other chemical reagents were analytical pure grade and obtained from Sigma–Aldrich (St Louis, MO, USA).

Synthesis of peptides and FITC-labeled peptides

RGERPPR, D(RPPREGR), D(RGERPPR), CRGERPPR, D(CRPPREGR) and D(CRGERPPR) were synthesized using Boc-protected solid phase peptide synthesis strategy on 4-methylbenzene and hydrogen amine resin. Resin was swelled with N,N-dimethylformamide for 20 min. The Boc-protected amino acid was activated by dissolving in HBTU and DIEA for 2 min. The activated amino acid solution was then added to the resin and mixed by stirring for 20 min. The amino group was deprotected by approximately twice the resin volume of trifluoroacetic acid and reacted with the carboxyl group of the next Boc-protected amino acid. This cycle was repeated to form the desired peptide chain. The peptide was then cut from resin by hydrogen fluoride. To prepare fluorescein-labeled peptides and decrease the interference of fluorescein, the N-terminals of RGERPPR, D(RPPREGR) and D(RGERPPR) were extended with cysteine. After purification to homogeneity via high-performance liquid chromatography (HPLC), all synthetic peptides were analyzed by ultraperformance liquid chromatography (UPLC) and electrospray ionization-mass spectrometry (ESI-MS). RGERPPR, D(RPPREGR), and D(RGERPPR) were also characterized by circular dichroism spectrum.

To synthesize FITC-CRGERPPR, a 1.2 × molar excess of Mal-FITC was dissolved in 1 mL of PBS with 5 mg CRGERPPR and stirred for 1 h in the dark. The reaction solution was purified via HPLC. The syntheses of other FITC-labeled peptides were the same as described for FITC-CRGERPPR. All FITC-labeled peptides were analyzed by UPLC and ESI-MS.

Surface plasmon resonance studies

To evaluate the binding affinity of RGERPPR, D(RPPREGR) and D(RGERPPR) with NRP-1, SPR interaction analyses were conducted using a Biacore T200 system (GE Healthcare, Uppsala, Sweden). First, the NRP-1 protein was coupled to the CM5 sensor chip via the standard amine coupling procedure. The peptides were then dissolved in HBS-EP (pH 4.0) buffer from low to high concentration (31.25 nmol/L to 1 μmol/L) and injected for detecting resonance changes. Data were analyzed by Biacore T200 evaluation software (GE Healthcare).

Computer-aided simulated design

The crystal structure model of NRP-1 protein was downloaded from www.pdb.org (pdb code: 5j1x). The missing atoms in the protein structure were added by the protein prepare package in the software Schrodinger Suites (2015, Schrödinger, LLC, New York, USA). The structure model of CendR peptides was built with the software Schrodinger and optimized in OPLS2005 force field. The design peptides were docked to NRP-1 protein by the Glide package in Schrodinger with default parameters19, 20, 21, 22. The binding energy was predicted by Glide. The resulting binding modes were analyzed, and the figures presented in paper were made by the software PyMOL (Version 1.3, Schrödinger, LLC, New York, USA).

Cell culture

Human lung adenocarcinoma epithelial cells (A549), human umbilical vein endothelial cells (HUVEC), human glioblastoma-astrocytoma cells (U87), and C6 rat malignant glioma cells were obtained from the Shanghai Institute of Cell Biology. Cells were cultured in DMEM containing 10% FBS, 100 U/mL penicillin and 100 mg/mL streptomycin at 37 °C under 5% CO2 atmosphere.

Cellular uptake and flow cytometry tests

HUVECs were seeded in 12-well plates at a density of 8 × 104 cells/well in 1 mL of DMEM with 10% FBS and incubated overnight (5% CO2) at 37 °C. FITC-CRGERPPR, FITC-D(CRPPREGR) and FITC-D(CRGERPPR) were added to the plates at a final concentration of 5 × 10−6 mol/L and incubated for 4 h. HUVECs were washed with PBS (0.01 mol/L, pH 7.4) twice and fixed with 10% formaldehyde for 10 min. The cells were washed with PBS twice and then overlaid with 10% buffered glycerol saline before observation using an inverted fluorescence microscope (Leica, DMI4000B, Wetzlar, Germany).

For flow cytometry evaluation, following the incubation described above, the HUVECs were trypsinized, centrifuged, resuspended in 300 μL PBS and analyzed by flow cytometry (FACSAria, BD, USA). Cellular uptake and flow cytometry tests of C6 cells, U87 cells and A549 cells were performed as the same as described for HUVECs.

Cytotoxicity assay

A549 cells were seeded into 96-well plates (5 × 103 cells/well) with 200 μL/well DMEM containing 10% FBS, 100 U/mL penicillin and 100 mg/mL streptomycin and incubated at 37 °C, 5% CO2 for 24 h. Various concentrations of gemcitabine or carboplatin in PBS were added into the wells, and plates were incubated at 37 °C, 5% CO2 for 48 h. Control cells were treated with fresh medium. Twenty microliters of 5 mg/mL MTT solution was added to each well and plates were incubated for 4 h. The medium was then removed, and 150 μL of DMSO was added to dissolve the formazan crystals. After 15 min, measurements were performed using a microplate reader (wavelength, 490 nm; PowerWave XS, Bio-TEK, Winooski, VT, USA). The measurements for each sample were conducted in triplicate.

Animal model establishment

Male nude mice (6–8 weeks old, 18–22 g weight) were purchased from Shanghai Sippr-BK Laboratory Animal Co. Ltd. (Shanghai, China). Mice were kept under specific pathogen-free conditions and provided with sufficient food and water. All animal experiments were conducted in accordance with the guidelines evaluated and approved by the Ethics Committee of Fudan University. The tumor model was established by subcutaneous injection of A549 cells (7 × 106 cells), C6 cells (8 × 106 cells) or U87 cells (9 × 106 cells) in 100 μL PBS into the right flank region of nude mice.

Tumor-penetrating ability study

The C6, U87, and A549 tumor-bearing nude mice were intravenously injected with 1 mg of Evans blue in PBS, followed another injection 5 min later with 4 μmol/kg of D(RGERPPR) or iRGD in PBS. The control mice were treated with an equal volume of PBS. The animals were perfused through the heart with PBS 30 min later, and then tumor tissues were collected and cut into pieces. Evans blue was extracted by incubating the tissue pieces in N,N-dimethylformamide at 37 °C in a water-bath shaker kept for up to 24 h. After centrifugation (4000 rpm, 10 min, Centrifuge H1650-W, Xiangyi, Changsha, China), the absorbance of the supernatant was measured via an ultraviolet spectrophotometer at 600 nm.

In vivo treatment study

The A549 tumor-bearing nude mice were randomly divided into four treatment groups (n = 5) when the tumors reached an average volume of 170 mm3. Mice were injected intravenously via a tail vein twice a week for a total of 5 times with the following: 1) equal volume of PBS (the control group); 2) D(RGERPPR) (total dose 20 μmol/kg) in PBS; 3) the free form of gemcitabine (total dose 75 mg/kg) dissolved in PBS; 4) gemcitabine (total dose 75 mg/kg) in PBS combined with D(RGERPPR) (total dose 20 μmol/kg) in PBS. Every two days, tumor size was measured across its two perpendicular diameters using a vernier caliper. The volume (V) was calculated according to formulawhere a and b are the longest and shortest diameters of the tumors, respectively. Tumor suppression rate was calculated as follows:where Vc is the tumor volume of the control group and Vt is the tumor volume of the treatment group.

Statistical analysis

All data are represented as mean ± SEM for all quantitative data. Statistical analysis was performed with the GraphPad Prism software. Two-tailed Student׳s t-test was used to determine significant differences between two groups (*P < 0.05). Two-way analysis of variance (ANOVA) was used to determine significant differences between multiple groups (*P < 0.05, ****P < 0.0001). ns indicates no significant differences (P > 0.05).

Results

Characterization and evaluation of peptides

Using a solid phase peptide synthesis strategy, we successfully obtained adequate quantities of highly pure RGERPPR, D(RPPREGR), D(RGERPPR) and FITC-labeled peptides. The HPLC and ESI-MS data (Supplementary Information Fig. S1) confirmed high purity levels and the molecular weight fit the theoretical value. RGERPPR, D(RPPREGR) and D(RGERPPR) were also characterized by circular dichroism spectrum (Fig. 1). All three peptides had an absorption peak at about 200 nm, indicating that the three peptides are disordered linear peptides. Moreover, the absorption peak of D(RGERPPR) was enantiomorphous with RGERPPR, indicating that the optical rotation of D(RGERPPR) was opposite to RGERPPR, proving that the peptides had different configurations. To evaluate the binding affinity of peptides to the NRP-1, we conducted SPR interaction analyses using a Biacore T200 optical biosensor. The D(RGERPPR) peptide registered a Kd value of 44.3 nmol/L, which was approximately 7.7 and 2 times lower than that of RGERPPR and D(RPPREGR), respectively (Table 1 and Supplementary information Fig. S2), indicating D(RGERPPR) possessed higher binding affinity to NRP-1.

Figure 1

Circular dichroism spectrum of RGERPPR, D(RPPREGR) and D(RGERPPR). The peptides were dissolved in deionized water at a concentration of 0.5 mg/mL for detection.

Table 1

Experimentally determined and predicted free energies of binding of the peptides with NRP-1 (ΔG = 1.3636 logKd).

Peptide	Kd (nmol/L)	Experimental ΔG (kcal/mol)	Predicted ΔG (kcal/mol)
RGERPPR	580	–8.50	–8.03
D(RPPREGR)	88.5	–9.61	–8.58
D(RGERPPR)	44.3	–10.02	–9.31

Computational analyses

Here we designed a peptide comprised of d-amino acids, D(RGERPPR). Notably, although the side chain of residues in D(RGERPPR) and D(RPPREGR) adopted opposite orientations compared with the reported active l-peptide, D(RGERPPR) and D(RPPREGR) showed higher binding affinity to the NRP-1. We performed a molecular docking study with the Schrodinger software package to explore the binding modes of the active peptides and determine the reason for the higher binding affinity of the two peptides comprised of d-amino acids. The three peptides were docked into the binding pocket of NRP-1 with the peptide docking package in Schrodinger. Examination of the binding modes of the peptides revealed that RGERPPR adopted the classic binding mode of the C-end peptide, in which the Arg7 at the C-terminal formed an ion interaction with Asp320 in NRP-1, and the COO– at the C-terminal formed a hydrogen bond network with Tyr353, Ser346, Thr349 and Trp301. Arg1 and Arg4 formed hydrogen bonds with the backbone of NRP-1. Arg7 and Arg1 formed cation–π interactions with Tyr353 and Tyr297 (Fig. 2 A and B).

Figure 2

Interaction of CendR peptides with receptors. (A), (C) and (E) NRP-1 is represented by the electrostatic potential surface. The designed peptides are shown in stick models and colored as green, blue and yellow for RGERPPR, D(RPPREGR) and D(RGERPPR), respectively. (B), (D) and (F) NRP-1 is represented by cartoon models. The interacting residues in the binding pocket and the peptides are shown in stick models.

In the binding mode of D(RPPREGR), the DArg1 at the N-terminal formed an ion interaction with Glu319 in NRP-1, and the NH3+ at the N-terminal interacted with Asp320. DArg4 formed a cation–π interaction with Tyr297 and Tyr353. DArg7 formed an ion interaction with Glu348 (Fig. 2C and D).

In the binding mode of D(RGERPPR), the three DArg residues also formed a strong hydrogen network and ion interaction with the residues of NRP-1. DArg7 interacted with negatively charged Glu324 and Glu319. DArg4 formed hydrogen bonds with the backbone of NRP-1. DArg1 in D(RGERPPR) formed ion interactions with Asp320. DArg1 and DArg4 also formed cation–π interactions with Tyr353 and Tyr297 (Fig. 2E and F).

In the sequence of the peptides, the Arg or DArg residues contributed to the primary binding affinity of the peptides. Comparing the binding modes of the three peptides, the three DArg residues in D(RGERPPR) and D(RPPREGR) fit well with the binding pocket and effectively interact with residues in NRP-1 by hydrogen bond, ion interaction and cation–π interaction. In the sequence RGERPPR, only Arg7 at the C-terminal could form strong interaction with the residues of NRP-1, and Arg1 and Arg4 formed weak interaction with the backbone of NRP-1. We predicted the binding affinity of the three peptides (Table 1), which indicated that D(RGERPPR) and D(RPPREGR) were more potent than RGERPPR, and D(RGERPPR) was slightly stronger than D(RPPREGR) in terms of binding affinity. The computational predicted result was also in good agreement with SPR-determined results. These results demonstrate that D(RGERPPR) and D(RPPREGR) possessed stronger binding affinity to NRP-1 compared with RGERPPR.

In vitro cellular uptake of peptides

The cellular uptake and flow cytometry data for FITC-CRGERPPR, FITC-D(CRPPREGR), and FITC-D(CRGERPPR) were shown in Fig. 3. In the cell uptake experiments, higher fluorescent signals of FITC-D(CRGERPPR) were observed in HUVEC, C6, U87 and A549 cells compared with FITC-D(CRPPREGR) and FITC-CRGERPPR. The mean fluorescent intensity of FITC-D(CRGERPPR) in HUVECs was nearly 13.1 times higher than that of FITC-CRGERPPR, and 5.7 times higher than that of FITC-D(CRPPREGR). In C6 and U87 cells, the mean fluorescent intensity of FITC-D(CRGERPPR) was 5.4 times and 3.3 times higher than that of FITC-CRGERPPR, and 2.9 times and 3.5 times higher than FITC-D(CRPPREGR), respectively. FITC-positive ratio of FITC-D(CRGERPPR) in A549 cells was 1.3 times and 1.6 times higher than that of FITC-D(CRPPREGR) and FITC-CRGERPPR. These data indicated that the cellular uptake efficiencies of D(RGERPPR) by HUVEC, C6, U87 and A549 cells were higher than that of both RGERPPR and D(RPPREGR), which was consistent with the result of SPR interaction analyses above. As HUVECs can simulate tumor neovascularization, we speculated that D(RGERPPR) possessed the capacity to target tumor vessels. Thus, D(RGERPPR) might have broad-spectrum tumor-penetrating ability.

Figure 3

Representative images of cellular uptake and flow cytometry data for FITC-CRGERPPR, FITC-D(CRPPREGR), and FITC-D(CRGERPPR) in various cell lines as indicated. Percentages and numbers below in the flow cytometry data represent FITC-positive ratios and mean fluorescent intensities, respectively.

Tumor-penetrating ability of(RGERPPR)

Based on aforementioned investigations in vitro, we can conclude that D(RGERPPR) possessed higher NRP-1 binding affinity and better cellular uptake efficiencies than other two peptides. Next, we investigated the tumor-penetrating ability of D(RGERPPR) in vivo.

Evans blue is a water-soluble albumin-binding dye, which was used to assess the integrity of blood vessels as well as the tumor-penetrating effect mediated by NRP-1. Sugahara et al. designed a cyclic peptide containing the CendR motif, iRGD (CRGDKGPDC), possessing tumor-penetrating ability. The iRGD peptide was cleaved by protease to expose the cryptic CendR motif: CRGDK. The CendR peptide then mediated binding to NRP-1 resulting enhanced permeability of tumor blood vessels and tumor tissues, which was confirmed via quantification of tumor-specific entry of Evans blue into tumor tissues. In this study, D(RGERPPR) co-administrated with the Evans blue dye, caused tumor-specific accumulation of the dye in all 3 tumor models we tested (Fig. 4). Compared with the control PBS group, D(RGERPPR) increased the Evans blue permeation by 36.4%, 30.5% and 46.4% in C6, U87 and A549 subcutaneous tumor models, respectively. By contrast, iRGD increased the Evans blue permeation by 49.7%, 17.7% and 32.9% in C6, U87 and A549 subcutaneous tumor models, respectively. According to a previous report, iRGD could remarkably increase the permeability of tumor tissues. In our research, we found that the tumor-penetrating ability of D(RGERPPR) was comparable to that of iRGD in C6, U87 and A549 tumor-bearing nude mice. Furthermore, the Evans blue average permeation of PBS group in A549 tumor model was 3.04 times and 12.25 times higher than that of C6 tumor model and U87 tumor model, indicating that tumor tissue of A549 tumor model has better permeability. And the tumor-penetrating ability of D(RGERPPR) was the highest in A549 tumor-bearing mice among three tumor models. Based on the above considerations, the in vivo treatment study was then performed on the A549 subcutaneous tumor model.

Figure 4

Quantification of tumor-specific entry of Evans blue into different subcutaneous tumor tissues in d(RGERPPR)-injected tumor-bearing mice (n = 3).

In vivo anti-tumor efficiency

We examined gemcitabine and carboplatin, as these are mainly used as first-line therapy in clinical treatment for non-small-cell lung cancer. The in vitro cytotoxicity of gemcitabine and carboplatin were measured by MTT assays. The IC50 values were 4.83 × 10–6 mol/L for gemcitabine and 1.38 × 10–4 mol/L for carboplatin after 48 h treatment (Supplementary Information Fig. S3), indicating that gemcitabine exhibited higher growth inhibitory effects in A549 cells than carboplatin.

The in vivo antitumor efficiency was next investigated on nude mice bearing A549 cell-derived xenograft tumors. PBS, D(RGERPPR), gemcitabine, or gemcitabine combined with D(RGERPPR) were injected via the tail vein. To evaluate tumor growth, tumor volumes (V) were divided by their primary tumor volume (V0) (Fig. 5). Compared with the PBS treatment group, the gemcitabine combined with D(RGERPPR) treatment was effective in retarding tumor growth, as shown in Fig. 5. The rate of tumor growth in the co-administration group was significantly lower than that of the gemcitabine alone group (P < 0.05), with a tumor suppression rate (TSR%) of 55.4%. The enhancement of anti-tumor effects could be attributed to the tumor-penetrating ability of D(RGERPPR), as evidenced above.

Figure 5

Tumor volume analysis of A549 tumor-bearing nude mice after i.v. injection of PBS, D(RGERPPR), free gemcitabine and gemcitabine combined with D(RGERPPR) (n = 5). Administration started on Day 0 and was repeated twice a week for a total of 5 times. Tumor volume was measured every other two days until Day 15. V indicates tumor volume at Day n and V0 indicates tumor volume at Day 0.

Discussion

Peptides possess low immunogenicity and high specificity, and are typically easy to synthesize with high quality. Consequently, peptides are widely applied in the treatment of tumors as well as cardiovascular, immune, and hematological diseases, and are also used as recognition molecules for drug delivery systems. However, natural bioactive peptides are susceptible to protease degradation and readily lose bioactivity in vivo. In addition, the chains of natural peptides are able to freely rotate and change shapes. Hence, they can completely or partially bind to various receptors, leading to low selectivity. The modification of the structure of bioactive peptides to enhance their stability, specificity, and selectivity has been a research hotspot. Generally, structural transformation of peptides mainly involves the following approaches: 1) substitution of unnatural amino acid, including d-amino acid, N-methyl-α-amino acid, β-amino acid, etc.; 2) cyclization of the peptide sequence; 3) isosteric, or not, amide bond replacement between two amino acids; 4) modification of N- or C-terminal ends; and 5) retro-inverso peptides. Among these methods, substitution of d-amino acids has remarkable superiority. First, because of the lack of protease that can degrade d-peptides in vivo, substitution of d-amino acid can enhance the biological stability and prolong the plasma half-life. Second, the substitution of d-amino acid also changes the conformation of peptides, which is beneficial to screen the bioactive peptides with high specificity and affinity to target protein.

In 1979, Goodman and Chorev established a retro-inverso peptide method that involved d-amino acids assembled in the reverse order from that in the parent l-sequence. Since the side chain topology of the retro-inverso peptide was similar to that of the parent l-peptide, the retro-inverso peptide might have similar bioactivities as the parent l-peptides. In our previous studies, we selected a CendR peptide (sequence RGERPPR) as the parent l-peptide, and designed its retro-inverso isomer, D(RPPREGR). Our results confirmed that D (RPPREGR) not only exhibited improved cellular uptake efficiency, but also possessed higher tumor-penetrating ability and better biological stability compared to RGERPPR.

Since d-peptide shows opposite chirality to the parent l-peptide with the same sequence, the inverso peptide might not have similar bioactivity. However, CendR peptides were a series of specific peptides with R/KXXR/K sequence, of which inverso isomer owned a D(R/KXXR/K) sequence. Since the C-terminal amino acid of the inverso isomer of CendR peptides was also arginine or lysine, we thus speculated that the inverso isomer also had NRP-1 binding activity and tumor-penetrating ability. In this study, we selected RGERPPR as the parent l-peptide and assembled d-amino acid residues in the order of the parent l-peptide sequence to synthesize the inverso peptide D(RGERPPR). As expected, the D(RGERPPR) possessed improved bio-stability (Supplementary Information Fig. S4). And the circular dichroism spectrum results showed that the absorption peak of D(RGERPPR) was enantiomorphous with RGERPPR, indicating that the optical rotation of D(RGERPPR) was opposite to RGERPPR. Surprisingly, the binding affinity of D(RGERPPR) with NRP-1 was significantly higher than that of RGERPPR. This result could be explained by computer simulation analysis as following: the three DArg residues in D(RGERPPR) fit well with the binding pocket and effectively interact with residues in NRP-1 by hydrogen bond, ion interaction and cation–π interaction. However, in the sequence RGERPPR, only Arg7 at C-terminal could interact well with the residues of NRP-1, and Arg1 and Arg4 formed weak interaction with the backbone of NRP-1. Therefore, D(RGERPPR) possessed stronger binding affinity to NRP-1 than RGERPPR. In the cellular uptake test, the cell uptake amount of D(RGERPPR) was significantly higher than that of FITC-labeled RGERPPR. We also investigated the tumor-penetrating ability by Evans blue test and found that the tumor-penetrating ability of D(RGERPPR) was comparable to iRGD in C6, U87 and A549 tumor-bearing nude mice. Due to its excellent biological activity, D(RGERPPR) could accumulate at tumor tissues and bind to NRP-1 on the surface of tumor cells. By binding NRP-1, D(RGERPPR) enhanced the permeability of tumor blood vessels and tumor tissues, which could facilitate anti-tumor agents to infiltrate into the tumor tissues.

Conclusions

By binding NRP-1, the CendR peptides can induce enhanced permeability of tumor blood vessels and tumor tissues. In this study, we selected a CendR peptide (sequence RGERPPR) as the parent l-peptide and synthesized its inverso isomer D(RGERPPR). The results indicated that D(RGERPPR) possessed high binding affinity to NRP-1 and improved tumor cell internalization, which could enhance the permeability of tumor blood vessels and tumor tissues, as well as facilitate anti-tumor agents to penetrate deep into tumor parenchyma. In conclusion, our findings demonstrate that D(RGERPPR) is a potential tumor-penetrating peptide and should be explored in further studies.