Introduction

Front. Chem.

Frontiers in Chemistry

Front. Chem.

2296-2646

Frontiers Media S.A.

10.3389/fchem.2019.00170

Chemistry

Original Research

Rational Design of Antiangiogenic Helical Oligopeptides Targeting the Vascular Endothelial Growth Factor Receptors

Zanella

Simone

¹ ^† Bocchinfuso

Gianfranco

² ^† De Zotti

Marta

³ ^† Arosio

Daniela

⁴ Marino

Franca

⁵ Raniolo

Stefano

² Pignataro

Luca

¹ Sacco

Giovanni

¹ Palleschi

Antonio

² Siano

Alvaro S.

⁶ Piarulli

Umberto

⁵ Belvisi

Laura

¹ ⁴ Formaggio

Fernando

³ ^* Gennari

Cesare

¹ ⁴ ^* Stella

Lorenzo

² ^*

¹Department of Chemistry, University of Milan, Milan, Italy ²Department of Chemical Science and Technologies, University of Rome Tor Vergata, Rome, Italy ³Padova Unit, Department of Chemistry, Institute of Biomolecular Chemistry, CNR, University of Padova, Padova, Italy ⁴National Research Council, Institute of Molecular Science and Technologies, Milan, Italy ⁵Center for Research in Medical Pharmacology, University of Insubria, Varese, Italy ⁶Departamento de Química Organica, Facultad de Bioquímica y Ciencias Biologicas, Universidad Nacional del Litoral, Santa Fe, Argentina

Edited by: Alessandro Contini, University of Milan, Italy

Reviewed by: Maria Jesús Pérez De Vega, Instituto de Química Médica (IQM), Spain; Nicolas Inguimbert, Université de Perpignan Via Domitia, France

*Correspondence: Fernando Formaggio fernando.formaggio@unipd.it Cesare Gennari cesare.gennari@unimi.it Lorenzo Stella stella@stc.uniroma2.it

This article was submitted to Chemical Biology, a section of the journal Frontiers in Chemistry

†These authors have contributed equally to this work

29 03 2019

2019

170

11 01 2019 05 03 2019

2019

Zanella, Bocchinfuso, De Zotti, Arosio, Marino, Raniolo, Pignataro, Sacco, Palleschi, Siano, Piarulli, Belvisi, Formaggio, Gennari and Stella

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Tumor angiogenesis, essential for cancer development, is regulated mainly by vascular endothelial growth factors (VEGFs) and their receptors (VEGFRs), which are overexpressed in cancer cells. Therefore, the VEGF/VEGFR interaction represents a promising pharmaceutical target to fight cancer progression. The VEGF surface interacting with VEGFRs comprises a short α-helix. In this work, helical oligopeptides mimicking the VEGF-C helix were rationally designed based on structural analyses and computational studies. The helical conformation was stabilized by optimizing intramolecular interactions and by introducing helix-inducing C^α,α-disubstituted amino acids. The conformational features of the synthetic peptides were characterized by circular dichroism and nuclear magnetic resonance, and their receptor binding properties and antiangiogenic activity were determined. The best hits exhibited antiangiogenic activity in vitro at nanomolar concentrations and were resistant to proteolytic degradation.

helical folded peptides protein-protein interactions C^α,α-disubstituted amino acids VEGF-C angiogenesis

Ministero dell’Istruzione, dell’Università e della Ricerca10.13039/501100003407

Associazione Italiana per la Ricerca sul Cancro10.13039/501100005010

Partnership for Advanced Computing in Europe AISBL10.13039/501100001943

Introduction

Angiogenesis—i.e., the formation of new blood vasculature from the established blood vessel network—can be associated to both physiological and pathological processes (e.g., inflammation, tumor growth, and metastasis). Among the latter, tumor angiogenesis is essential for cancer development, since neovascularization provides a steady supply of oxygen and nutrients, supporting the proliferation of cancer cells (Mizejewski, 1999; Danhier et al., 2012; Johannessen et al., 2013). For this reason, antiangiogenic agents are used to impede or retard cancer progression and metastasis (Ferrara and Adamis, 2016).

Different receptors are involved in angiogenesis regulation, and their expression depends on the conditions of the cell environment (e.g., pH, oxygen or supply of nutrients) (Mizejewski, 1999). Angiogenic processes are also mediated by cross-talk mechanisms that trigger direct association and cluster formation between specific receptors. For example, the cooperation between integrins and vascular endothelial growth factor receptors (VEGFRs) was found to be crucial in pathological processes such as tumor growth and development (Somanath et al., 2009; Desgrosellier and Cheresh, 2010).

VEGFRs are receptor tyrosine kinases (RTKs) that have a central role in tumor angiogenesis and progression. Indeed, the hypoxic conditions characteristic of the tumor environment induce both up-regulation of VEGFRs and gene expression of vascular endothelial growth factors (VEGFs) (Ferrara et al., 2003; Hoeben et al., 2004; Olsson et al., 2006; Carmeliet and Jain, 2011; Shibuya, 2013). In addition to the tumor promoting effects of neovascularization, autocrine VEGF/VEGFR signaling favors growth, proliferation and migration of cancer cells (Su et al., 2007; Simon et al., 2017).

Dimerization and activation of VEGFRs are triggered by binding of VEGFs to the extracellular domain of the receptors (Ferrara and Adamis, 2016). The VEGF family consists of five members [VEGF-A, VEGF-B, VEGF-C, VEGF-D, and PlGF (placental growth factor)], and there are three distinct receptors (VEGFR-1, VEGFR-2, VEGFR-3). The various growth factors have differential selectivity: for instance, VEGF-A binds to VEGFR-1 and VEGFR-2, while VEGF-C binds to VEGFR-2 and VEGFR-3. In addition, the three receptors have different functions, with VEGFR-2 being mainly involved in angiogenesis and VEGFR-3 regulating lymphangiogenesis (Ferrara and Adamis, 2016; Nasir, 2019).

Three main approaches targeting VEGF-A/VEGFR-2 signaling in human cancer have been approved for clinical practice (Nasir, 2019). One strategy involves inhibition of the tyrosine kinase activity of VEGFR-2 by small molecules interacting with the intracellular segment of the receptor, such as sorafenib and sunitinib (Musumeci et al., 2012). Alternatively, VEGF-mediated angiogenesis can be impaired by blocking the VEGF-A/VEGFR-2 interaction. This has been accomplished by binding and neutralizing circulating VEGF-A with monoclonal antibodies (bevacizumab) or recombinant proteins mimicking the receptor (aflibercept) (Ferrara et al., 2004; Ferrara and Adamis, 2016). Alternatively, VEGFR-2 has been targeted with therapeutic monoclonal antibodies, such as ramucirumab. These molecules are now a standard of care for the treatment of several metastatic cancers (Ferrara and Adamis, 2016), even though clinical results have not met in all cancer types the initial high hopes for this therapeutic strategy, underlining the need for more effective antiangiogenic drugs (Vasudev and Reynolds, 2014; Ronca et al., 2017). Anti-VEGF-A/VEGFR-2 therapy has proven more effective in ophthalmology, in the treatment of intraocular neo vascular disorders, such as age-related macular degeneration (Ferrara and Adamis, 2016). However, each of the above mentioned approved drugs has its own limitations: kinase inhibitors have limited selectivity, while antibodies suffer from poor pharmacokinetics, limited tissue penetration and high costs (Chames et al., 2009; Howard et al., 2015). Their use in the therapy of eye diseases requires repeated intravitreal injections, and new drugs that can be administered by simpler and safer routes are highly desirable (Sidman et al., 2015).

Peptides and peptidomimetics with a well-defined conformation (foldamers) represent a promising alternative to biological therapeutics for the inhibition of protein-protein interactions, and are currently experiencing a revival of interest from the pharmaceutical industry (Henninot et al., 2018). Among those able to interact with the extracellular domain of VEGFRs, it is worth mentioning the peptoid ligands described by Kodadek and colleagues (Udugamasooriya et al., 2008), the helical peptides developed starting from VEGF-A and Vammin hotspots (García-Aranda et al., 2013), and cyclopeptides isolated by phage display technique (Zilberberg et al., 2003), or rationally designed (Gautier et al., 2010). In 2011, D'Andrea and coworkers developed a α-helical decapentapeptide based on the natural sequence of VEGF-A, with potent inhibitory activity against VEGF-stimulated angiogenesis in vivo (Basile et al., 2011; Diana et al., 2013). Because of its synthetic accessibility and antiangiogenic properties, recently some of us selected that peptide to prepare a dual-action compound able to interfere with the integrin α_Vβ₃-VEGFR-1 cross-talk (Zanella et al., 2015). This conjugate was able to bind in vitro both integrin α_Vβ₃ and VEGFR-1, and exerted a strong antiangiogenic effect in VEGF-stimulated morphogenesis assays on human umbilical vein endothelial cells (HUVECs).

Herein, we report the results of our efforts to develop new VEGFR antagonists based on a helical fragment of VEGF-C, with promising activity against VEGF-mediated angiogenesis.

Results Design A Short Helix Is an Important Element of the VEGF/VEGFR Interface

The extracellular domain of VEGFRs consists of seven Ig homology domains. VEGF binding takes place mostly on domain D2. Figure 1 shows the structure of the VEGF-A/VEGFR-1 (D2 domain) interface, whose main features are conserved in all VEGF/VEGFR complexes (Leppänen et al., 2013). We started our analysis from this complex, because all currently approved antiangiogenic drugs are targeted to VEGF-A or to VEGFR-1. VEGF residues are colored from green to red in order of increasing penalty in the standard binding free energy caused by their substitution to Ala, as predicted by in silico alanine scanning. This analysis indicates that a short helix, located at the N-terminal region of the growth factor, represents a significant portion of the interaction interface. Previous experimental studies have shown that this helix is involved in receptor specificity and protein dimerization of VEGFs (Siemeister et al., 1998; Robinson and Stringer, 2001; Leppänen et al., 2010). In the case of VEGF-A, the helix comprises residues 17–25 (a nonapeptide). Table 1 reports the numerical results of the in silico Ala scan for those residues, showing that the interaction with the receptor is mediated principally by amino acid residues at positions 1, 2, 5, 6, and 9 of the helical nonapeptide. This finding is consistent with experimental Ala scan data (Muller et al., 1997; Li et al., 2000). A similar analysis performed on the structure of the VEGF-A/VEGFR-2 complex confirmed the N-terminal helix as an important interacting element (Figure S1).

Figure 1

Structure of a VEGF/VEGFR complex (PDB code 1FLT). Ig-homology domain 2 of the receptor (VEGFR-1) is shown in surface representation (light blue and blue for the two subunits), while the two subunits of VEGF-A are shown as ribbons of two different shades of green. VEGF residues are colored depending on the standard binding free energy penalty associated with their mutation to Ala, as predicted by DrugScorePPI (values from 0 to 1.5 kcal/mol are reported in shades of color going from green to red; values 1.5 kcal/mol and above, up to the maximum of 2.4 kcal/mol, are all reported in red).

Table 1

In silico Ala scan of the VEGF-A helix interacting with VEGFR-1 (PDB code 1FLT)^a.

Position in the helix	Residue	FoldX	DrugScorePPI	KFC2	PCRPi
		ΔΔG binding° (kcal/mol)	ΔΔG binding° (kcal/mol)	Hot spots	Probability (%)
1	Phe 17	1.6	0.4	X	48
2	Met 18	2.0	0.4	X	8
3	Asp 19	0.1	–	–	0
4	Val 20	0.1	–	–	0
5	Tyr 21	2.4	2.4	X	14
6	Gln 22	0.7	0.3	–	2
7	Arg 23	0.2	0.1	–	0
8	Ser 24	0.2	–	–	0
9	Tyr 25	0.5	1.4	X	4

FoldX and DrugScorePPI predict the penalty in standard binding free energy for substitution of each residue to Ala, the KFC2 server simply predicts hotspots and PCRPi reports the probability that a given residue is a hotspot. Substitutions that are predicted to have a strong effect on the binding are colored in red (i.e., those with a v ΔΔGbinding° higher than 1.0 kcal/mol or indicated by KFC2 as hot spots or with a PCRPi probability greater than 10). Substitutions that are predicted to have a mild effect on the binding are colored in light red (i.e., those with ΔΔGbinding° comprised between 0.3 and 1.0 kcal/mol or with a PCRPi probability comprsed between 1 and 10)

The Helix of VEGF-C Presents the Most Promising Inter- and Intra-molecular Interactions

Figure 2 shows the corresponding helical sequences in various growth factors, together with the available crystallographic structures of their complexes with receptors. In most cases, the interaction is largely based on hydrophobic effects. The only significant exception is provided by VEGF-C, which forms an intermolecular salt bridge with Arg164 of VEGFR-2 through the Asp residue at position 2 of the helix. This specific interaction is predicted to improve the selective recognition of the VEGF-C helix by VEGFR-2. Figure S1 shows the results of an in silico Ala scan performed on the VEGFR-2/VEGF-C complex, confirming the centrality of the N-terminal helix in the recognition process.

Figure 2

Structures of the helix-interacting interface in VEGF/VEGFR receptor complexes. The receptor is shown in light blue and the growth factor helix in green. The following crystallographic structures were used: 1FLT for VEGF-A/VEGFR-1, 3V2A for VEGF-A/VEGFR-2, 2XAC for VEGF-B/VEGFR-1, 2X1X for VEGF-C/VEGFR-2, 4BSK for VEGF-C/VEGFR-3, 3VSB for VEGF-E/VEGFR-2, 1RV6 for PlGF/VEGFR-1. Sequences of the helix correspond to residues 17–25 for VEGF-A, -B and -E, residues 122–130 for VEGF-C, and residues 25–33 for PlGF. In the sequences, reported on top in the single letter code, interacting residues are underlined, bolded, and colored, according to the following code: green, purple, red, and blue for hydrophobic, polar, anionic, and cationic residues, respectively. The ion bridge formed in the VEGF-C/VEGFR-2 complex is indicated by a purple line.

The VEGF-C helix presents also another interesting property. Short helical segments are usually largely disordered when separated from the protein that contained them. Without the stabilization provided by the rest of the protein structure, competition by water molecules breaks the intramolecular H-bond network that holds the helical conformation together. On the other hand, binding to the receptor requires a helical conformation. Therefore, a disordered peptide pays a significant entropic cost for association to the receptor. Thus, stabilization of the helical conformation in such short peptides is advisable. Interestingly, the helix of VEGF-C presents three pairs of side chains that potentially form interactions stabilizing a helical conformation, being located at an i - i+3 or i - i+4 distance. In particular, Ile1 and Trp5 can form a hydrophobic cluster, while both Asp2 and Arg6 and Glu4 and Lys7 can form salt bridges (numbers correspond to positions in the helix).

Replica exchange molecular dynamics (REMD) simulations were used to verify the stability of the VEGF-C helix, when free in solution. Figure 3 reports the secondary structure of the peptide, a representative conformation and the histograms of the distances between the side-chains mentioned above during the simulation trajectory. The data confirm that the ion pairs and hydrophobic cluster are significantly populated and stabilize the helical conformation of the peptide in solution.

Figure 3

Stability of intramolecular interactions and secondary structure in the REMD simulation of peptide 1 (300 K replica). Upper: distributions of the minimum distances between the charged groups of Asp123 and Arg127 (red) and of Glu125 and Lys128 (orange) and between the side chains of Ile122 and Trp126 (green). Residue numbers refer to the VEGF-C sequence. The first ns of the simulation was excluded from the analysis. Lower: secondary structure according to DSSP. Blue: alpha helix, grey: 3₁₀ helix, yellow: turn, green: bend, white: coil.

For all these reasons, we decided to focus on the VEGF-C helix sequence (peptide 1) (Figure 4) to develop VEGF/VEGFR interaction inhibitors.

Figure 4

Sequences of the VEGF-C natural portion 1 and of the derived peptides 2–6. Sequences are colored in green, purple, red, and blue for hydrophobic, polar, anionic, and cationic residues. The introduced modifications are underlined.

Substitution of Two Residues in the VEGF-C Sequence Is Predicted to Increase Binding Affinity

In order to optimize peptide affinity for VEGFRs, in silico mutagenesis was performed, concentrating on the VEGF-C/VEGFR-2 complex (PDB code 2X1X). All residues of the helix previously identified as hot spots were mutated to all possible coded amino acids, evaluating the effect of the substitution on the standard binding free energy. This analysis identified substitution of the first residue from Ile to Trp, and of the last one from Gln to Trp or Val as the only mutations that would lead to a significant predicted improvement in standard binding free energy (>1 kcal/mol). Either hydrophobic substitution at position 1 would maintain the helix-stabilizing cluster identified above, but Trp has a higher intrinsic helix propensity than Ile (Chakrabartty et al., 1994; Pace and Scholtz, 1998). Based on these findings, sequence 2 was designed (Figure 4). Figure 5 shows how the N- and C-terminal Trp residues increase the interactions between helix and receptor.

Figure 5

Structures of peptides 1 (left) and 2 (right), interacting with VEGFR-2. The structure of peptide 1 was taken from the VEGF-C structure (PDB code 2X1X), while that of peptide 2 was modeled with FoldX.

A more comprehensive computational analysis on the proposed substitutions, carried out by molecular mechanics—Poisson Boltzmann surface area calculations (MM-PBSA) (Hou et al., 2010), confirmed the prediction, showing that the presence of Trp residues at each peptide terminus should indeed lead to an increased binding affinity (Table 2).

Table 2

Results of the MM-PBSA calculations of the interaction of peptides 1 and 2 with VEGFR-2^a.

Peptide	ΔE_pot (kcal/mol)	ΔG°_PB(kcal/mol)	ΔSAS (nm²)	ΔG°_binding (kcal/mol)
1	−32 ± 1	26 ± 1	−3.59 ± 0.02	−8 ± 2
2	−39 ± 1	27 ± 1	−4.42 ± 0.02	−14 ± 2

The various terms are described in the Methods section.

C<sup>α,α</sup>-Disubstituted Amino Acids Can Be Inserted to Increase Peptide Helicity and Stability

C^α,α-disubstituted α amino acids strongly stabilize helical structures: due to steric interactions between the gem alkyl and methyl groups linked to the α-carbon, the accessible conformational space of such residues is extremely limited and is located in the region of the Ramachandran plot corresponding to helical structures. For this reason, such C^α,α-disubstituted amino acids constitute suitable building blocks to synthesize oligopeptide with a stable helical folding (Toniolo et al., 2001). In addition, insertion of non-proteinogenic amino acids strongly reduces peptide susceptibility to proteolysis (De Zotti et al., 2009). The simplest and most studied C^α,α-disubstituted amino acid residue is α-amino isobutyric acid or Aib. Aib is a natural amino acid featuring two methyl groups on its α-carbon. It is non-ribosomally included in peptide sequences by fungal synthases. This achiral residue is a well-known helix inducer. One peptide sequence modified to include Aib is the well-known commercial drug semaglutide (Al Musaimi et al., 2018). Positions not involved in the interaction (i.e., 3, 4, 7, and 8 in the helix, see above) were considered for possible modifications. Position 2 was analyzed as well, as the C^α,α-disubstituted amino acids analogue of Asp, α-methyl-aspartic acid (α-Me)Asp, is commercially available. Although literature reports on the conformational properties of this specific residue are lacking, the preference for helical conformations is a common feature of methyl-containing disubstituted residues, and its homologue (α-Me)Asn was reported to promote type III β turn (i.e., a portion of a 3₁₀ helix) in short peptides (Hopkins et al., 2000). Analysis of the structure of the VEGF-C/VEGFR-2 complex showed that addition of a methyl group on the alpha carbon at these positions, or addition of a residue at the N-terminus, would not cause any intermolecular clashes.

Based on these considerations, four analogues were designed, based on sequence 2 (Figure 4). In analogues 5 and 6, Glu4 was substituted by (α-Me)Asp, and Lys7 by the cationic amino acid Api (4-aminopiperidine-4-carboxylic acid). Such cyclic residue can promote a helical conformation to some extent when incorporated into peptides (Cho et al., 2010; Dalzini et al., 2016). In addition, in peptide 5 Asp2 was substituted by (α-Me)Asp, too, while in 6 an Api residue was added at the N-terminus. In analogue 4, Aib was inserted at positions 4 and 8, and also added at the N-terminus. For comparison, analogue 3 has been synthesized, with the same sequence of 4, but without the Trp substitutions.

Synthesis of VEGF-C Derived Peptides 1-6

Peptide sequences 1–6 (Figure 4) were conveniently prepared by solid-phase peptide synthesis (SPPS) on Rink Amide 4-methylbenzhydrylamine (MHBA) resin using the 9-fluorenylmethoxycarbonyl/t-butyl ether (Fmoc/tBu) strategy. Each step of the SPPS was performed with a semi-automatic synthesizer, assisting coupling reactions with microwaves. First, the resin was swelled in dimethylformamide (DMF) and treated with a 25% solution of piperidine in DMF to remove the Fmoc-protecting group, releasing the reactive amino moiety on the beads (step a). The Fmoc-amino acid to be attached to the solid support was activated with N,N′-diisopropylcarbodiimide (DIC) and 1-hydroxy-7-azabenzotriazole (HOAt) coupling reagents in the presence of N,N-diisopropylethylamine (iPr₂NEt) in DMF: after stirring for 25 min at 0°C, the mixture was added to the resin and a cycle of coupling, capping and deprotection was performed (step b). This procedure was repeated until the sequence was completed, then the N-terminal residue was acetylated and the beads were treated with TFA in the presence of thioanisole, 1,2-ethanedithiol (EDT) and anisole as scavengers. Subsequent purification of the crude peptide with reversed phase high-performance liquid chromatography (RP-HPLC) and freeze-drying from glacial acetic acid gave the pure compound as a fluffy solid.

Due to the high steric congestion, α,α-disubstituted residues are poorly reactive even after activation. For this reason, the coupling reaction was performed twice whenever a quaternary amino acid had to be attached to the resin: according to this synthetic protocol, after the first coupling step, a second aliquot of the activated quaternary amino acid was added to the resin and another condensation reaction was carried out. This procedure allowed to minimize the number of unreacted amino moieties on the beads.

Further details are provided in the Supplementary Information (Figures S2–S14, Tables S1–S6, Scheme S1).

Structural Investigation

Circular dichroism (CD) experiments in water demonstrated that all analogues populate helical conformations in solution to some extent, as predicted by REMD simulations for the natural sequence 1. As shown in Figure 6, the positive peak at approximately 190 nm, and the negative bands at about 205 and 220 nm, typical of helical structures, are present in the spectra of all analogues (Kelly and Price, 1997). Quantitative determinations of the degree of helicity is complicated by the presence of multiple aromatic residues, whose side-chains can contribute significantly to the CD in the far UV region. However, several features indicate that analogues 5 and 6 are more helical than the other peptides: compared to the spectra of analogues 1–4, the positive peak is higher, the negative band at short wavelengths is red-shifted and the ratio of the two negative bands at approximately 220 and 205 nm is closer to one (Kelly and Price, 1997). Actually, for 5 this ratio is even higher than one, possibly indicating that some peptide aggregation is taking place (Dai et al., 2004). A stabilizing effect of the introduced C^α,α-disubstituted amino acids was less obvious in the spectra of analogues 3 and 4.

Figure 6

CD spectra of the peptides in distilled H₂O (10⁻⁴ m), in the far UV (left) and near UV (right).

In near UV CD spectra, the induced dichroism band of Trp was more intense for analogues 3–6 than for peptides 1 and 2, which comprise coded amino acids only. This finding indicates that the C^α,α-disubstituted amino acids rigidified the helical conformation. To confirm this conclusion also in the case where far UV CD spectra were less informative, analogues 1 and 3 were directly compared in 2D nuclear magnetic resonance (NMR) measurements in water (Figure 7). For both peptides, all sequential NH-NH cross peaks were present in the ROESY (Rotating-frame nuclear Overhauser Effect correlation SpectroscopY) spectrum, indicating the onset of a helical structure. However, the spectra did not show any long range connectivity for 1, while the presence of two long-range cross-peaks (Asp3HA-Aib5HN and Trp6HA-Aib9HN) for peptide 3, including an αH_i → NH_i+2 interaction, is a clear indication that this analogue adopts a well-developed 3₁₀-helical conformation in water.

Figure 7

2D proton NMR ROESY spectra in H₂O:D₂O 9:1 (600 MHz, 298K) for compounds 1 (left) and 3 (right) (concentrations 1.1 and 1.4 mm, respectively). Top panels report NH-NH sequential correlations, while the bottom panels correspond to the fingerprint region showing long-range αH-NH correlations. In the case of compound 3, long range cross-peaks, including the αH_i → NH_i+2 diagnostic of 3₁₀-helical structures, are highlighted in blue.

Overall, these data indicate that the introduction of C^α,α-disubstituted amino acids stabilized the helical conformations, particularly in analogs 5 and 6, where chiral (α-Me)Asp residues were inserted in the sequence.

Biological Studies Inhibition of VEGF-A/VEGFR-1 Complex Formation <italic>in vitro</italic>

A chemiluminescent assay to determine the inhibition of complex formation between VEGF-A (isoform VEGF₁₆₅) and the extracellular domain of VEGFR-1 (Goncalves et al., 2007; García-Aranda et al., 2013) had been previously employed in our group (Zanella et al., 2015). VEGF-C (from which analogue 1 was derived) is not selective for VEGFR-1 (Su et al., 2007), but unfortunately the same assay protocol turned out to be ineffective with VEGFR-2, at least in our hands. The assay based on VEGFR-1 was therefore applied here as a very stringent test of the affinity of our peptides for VEGFRs, and to assess the different peptides comparatively. All peptides were able to inhibit complex formation (Table 3), although in the high micromolar range. The introduction of quaternary amino acids resulted in a beneficial effect in terms of affinity toward VEGFR-1, with analogs 3, 5 and 6 being the most effective. Further studies were focused on analogues 5 and 6, considering also the higher stability of their secondary structure, which affects resistance to proteolytic degradation.

Table 3

Inhibition of biotinylated VEGF₁₆₅ binding to isolated VEGFR-1.

Peptide	% Inhibition^a
1	27 ± 1
2	29 ± 7
3	86 ± 1
4	40 ± 10
5	82 ± 1
6	83 ± 3

% of inhibition of biotinylated VEGF₁₆₅ binding to VEGFR-1 at 500 μM.

Resistance to Proteolytic Degradation

Peptides 5 and 6 were very stable against proteases, presumably thanks to the presence of non-natural amino acids in their sequences and to their stable secondary structure. HPLC analysis (Figures S15–S24) demonstrated that peptide 5 was fully stable to trypsin, chymotrypsin and pronase for 90 minutes and it persisted even after many days, although some degradation was slowly occurring. Peptide 6 was fully stable even after 6 days. By contrast, the natural sequence 1 was degraded within 15 min by all three enzymes.

Inhibition of HUVEC Morphogenesis

The ability of peptides 5 and 6 to affect neovessel formation in vitro was investigated on HUVECs according to the previously reported procedure (Fanelli et al., 2014). Experiments were performed in the presence of peptide 5 or 6 under resting (absence of stimuli) or stimulated (VEGF₁₆₅) conditions. As expected, under resting conditions, HUVECs did not show any significant network formation, while the presence of VEGF₁₆₅ induced a strong increase in loop and branches formation (Figure 8). Pre-incubation with either 5 or 6 significantly reduced VEGF₁₆₅-induced loop and branches formation and this effect was concentration-dependent, reaching statistical significance at 10 nm for loop and 1 μm for branches formation (Figure 8). Interestingly, these values are much lower than the concentrations needed to inhibit VEGF-A/VEGFR-1 association, possibly indicating a significant selectivity toward VEGFR-2. It is worth noting that inhibition of loop formation was more marked for analogue 5 than for 6, at all concentrations tested.

Figure 8

Antiangiogenic activity of analogues 5 and 6. The images are representative phase contrast photomicrographs of HUVECs plated on Matrigel under resting conditions (lower left) or after stimulation with VEGF₁₆₅, 10 ng/mL (upper left); the effect of VEGF on loop formation, in comparison to medium alone, is clearly evident in the upper panel. 5 (upper panels) and 6 (lower panels) were added at different concentration on VEGF-stimulated cells; both 5 and 6 are able to reduce, in a concentration-dependent manner, the ability of VEGF to induce loops formation. The graphs describe the effect of addition of 5 (purple) and 6 (green) on VEGF-induced neovessel formation measured as total number of loops (left graph) or total branch length (right graph). Data are expressed as mean ± SEM of 5 separate experiments. Asterisks indicate significant differences with respect to the positive control (basal medium + VEGF), with p < 0.05, according to a two-tailed paired t-test. Negative controls (basal medium only) are also shown.

Discussion

Several peptides able to bind to VEGFRs have been reported in the literature. Some of them were identified by library screening or phage display approaches, while others were developed by rational design, reproducing different epitopes of the VEGFR interacting surface of VEGFs. However, most of these studies were focused on mimicking VEGF-A, which is traditionally considered the master regulator of angiogenesis through VEGFR-2 binding (Ferrara and Adamis, 2016). In particular, all helical foldamers mimicking the N-terminal helix were developed based exclusively on the interacting elements of VEGF-A, i.e., hydrophobic residues only. Examples include the MA peptide (Diana et al., 2013) and peptides developed by Pérez De Vega and coworkers (García-Aranda et al., 2013; Balsera et al., 2017). Due to their hydrophobic driving force for association, these peptides are prone to selectivity issues, since binding sites for amphipathic helical peptides are extremely common at protein-protein interfaces (Bonache et al., 2014). It is also worth mentioning that some of the peptides modeled on the VEGF-A helix exhibited pro-angiogenic, rather than anti-angiogenic, activity (De Rosa et al., 2018).

In the present study, from an analysis of the available crystallographic structures of VEGFs/VEGFRs complexes, we noted that the same helix in VEGF-C exhibits several interesting features, including an intermolecular salt bridge with VEGFR-2, and multiple intramolecular helix-stabilizing interactions. For these reasons, we focused our study on the VEGF-C sequence.

VEGF-C binds to both VEGFR-2 and VEGFR-3 (Su et al., 2007; Leppänen et al., 2010; Chen et al., 2013; Wang and Tsai, 2015). Traditionally, regulation of lymphangiogenesis through association with VEGFR-3 has been considered its prevalent activity. Indeed, VEGF-C expression is closely related to lymphangiogenesis and lymphatic metastasis in a variety of human tumors (Chen et al., 2013). However, more recently, VEGF-C has been demonstrated to play an important role also in the regulation of physiological and pathological angiogenesis (Chen et al., 2013). In addition, it has fundamental functions in the autocrine signaling of cancer cells. VEGF-C is expressed in a number of human malignancies, and high levels of expression correlate with poor prognosis (Chen et al., 2013). Autocrine VEGF-C signaling regulates cell invasion, proliferation, and resistance to chemotherapy (Su et al., 2007). VEGF-C can also modulate the immune system so that tumor cells more easily escape immune surveillance (Wang and Tsai, 2015). A very recent, groundbreaking article (Michaelsen et al., 2018; Niclou, 2018), identified VEGF-C, rather than VEGF-A, as the main responsible for autocrine VEGFR-2 activation and cell proliferation in glioblastoma. In this study, VEGF-C silencing was superior to bevacizumab therapy in improving tumor control. Interestingly, bevacizumab treatment increased VEGF-C expression (Michaelsen et al., 2018), and up-regulation of VEGF-C has been observed in tumor cells that have acquired resistance to anti-VEGF-A therapy (Wang and Tsai, 2015), suggesting that VEGF-C may compensate for VEGF-A depletion. For all these reasons, molecules inhibiting VEGF-C signaling, rather than, or in addition to VEGF-A interactions, might find important therapeutic applications. In principle, our peptides, based on the VEGF-C helix, can compete with both growth factors for binding to the receptors.

The VEGF-C nonapeptide helix sequence has some intrinsic helix propensity, but CD and NMR studies demonstrated that insertion of C^α,α-disubstituted amino acids led to a significant stabilization of the secondary structure. In particular, analogues 5 and 6, containing the chiral (α-Me)Asp amino acids, had a more stable helical conformation than peptides 3 and 4, where achiral Aib residues were introduced. Peptides comprising these modifications were highly resistant to proteolytic degradation, and stabilization of the helical conformation probably contributed to favor receptor binding, by reducing the entropic cost of the association process. Indeed, these molecules were able to bind even receptor VEGFR-1, where specific ion-pair interactions are not possible, although in the high micromolar range. More importantly, they inhibited VEGF-induced morphogenesis in the low nm range, possibly through their interaction with the VEGFR-2 receptor.

One point worth mentioning is that the N-terminal VEGF helix, mimicked by our peptides, is involved also in the interface of the growth factor dimer (see Figure 1). Therefore, inhibition of VEGF dimerization could contribute to the antiangiogenic activity of the compounds developed here.

Further studies are warranted to characterize the interaction of the designed peptides with different VEGFRs and the molecular mechanisms of their antiangiogenic activity. However, our data clearly indicate that the development of antiangiogenic folded synthetic peptides inspired by the VEGF-C N-terminal helix might open the way to a novel class of anticancer agents.

Methods <italic>In silico</italic> Mutagenesis

in silico alanine scanning was performed with FoldX 4.0 (Schymkowitz et al., 2005) and the DrugScorePPI (Krüger and Gohlke, 2010), KFC2 (Zhu and Mitchell, 2011), and PCRPi (Segura Mora et al., 2010) servers. In silico mutagenesis was performed with FoldX, by first repairing the PDB file (PDB code 2X1X), then introducing the mutations and finally calculating the complex standard binding free energy, and subtracting the value of the WT complex. Structural images were obtained using the Chimera software (Pettersen et al., 2004).

Molecular Dynamics Simulations

Molecular dynamics (MD) simulations were carried out by using the GROMACS 4.6.3 software package (Hess et al., 2008). The ff53a6 parameters (Oostenbrink et al., 2004) were adopted for the complexes and the SPC model (Berendsen et al., 1981) was used for the water molecules. Short-range electrostatic interactions were cut-off at 1.2 nm and long range electrostatic interactions were calculated using the particle mesh Ewald (PME) algorithm (Essmann et al., 1995). Simulations were run with a 2 fs time step. The Berendsen algorithm was used to keep temperature and pressure constant (Berendsen et al., 1984). Bond lengths were constrained with the LINCS algorithm (Hess et al., 1997). The structure with PDB code 2X1X from the Protein Data Bank was used as a starting point for the simulations. After solvating the proteins with roughly 6500 water molecules and an opportune number of Na⁺ and Cl⁻ ion to ensure electroneutrality in a box of about 180 nm³, the potential energy of the systems was minimized. A multi-step procedure was adopted for equilibration, first constraining both the protein and the peptide, then releasing the constraints on the peptide conformation, and finally removing all constraints. In each step, the temperature was gradually raised from 50 K to 300 K in 8 ns, and the conformation with the most favorable peptide-protein interaction energy was selected as the starting structure for the successive step. Finally, 50 ns production runs were performed.

MM-PBSA Calculations

The protein-peptide binding energies were obtained from MD simulations by using the MM-PBSA protocol. The binding free energies were calculated according to the following equation

ΔGbinding°=Gcomplex°-Greceptor°-Gligand°

In these free energy values, the potential energy terms were obtained by using the GROMACS g_energy tool. The electrostatic solvation terms were calculated with the APBS software (Baker et al., 2001). The dielectric constant was set equal to 78.54 for water and 2 for the solute. The input files for APBS were generated with the PDB2PQR server (Dolinsky et al., 2004).

The non-polar solvation energies were computed as

Gsurf°=γSAS

in which γ was set equal to 0.543 kcal mol⁻¹ nm⁻² (Fogolari et al., 2003), solvent accessible surface (SAS) values were obtained by using the g_sas tool of GROMACS.

REMD Simulations

REMD simulations were performed using the GROMACS software package. An extended conformation was used as the starting structure. After minimization and equilibration, 16 replicas were simulated with temperatures ranging from 260 to 650 K, chosen to ensure an exchange probability between replicas equal to 55% (Patriksson and van der Spoel, 2008); exchanges were attempted every 2 ps and the simulation time of each replica was 50 ns. The Amber FF99-SB parameters were used to describe the peptide (Hornak et al., 2006) and the OBC (Onufriev, Bashford, and Case) GBSA implicit solvent model (Onufriev et al., 2004) was used to mime the solvent. Bond lengths were constrained with the SHAKE algorithm (Ryckaert et al., 1977) and simulations were run with a 2-fs time step. Temperatures were kept constant by using the velocity-rescale algorithm (Bussi et al., 2007). Secondary structure was assigned by means of Dictionary of Protein Secondary Structures (DSSP) (Kabsch and Sander, 1983). The g_mindist tool in GROMACS was used to calculate the minimum distances between groups of atoms.

SPPS

Semi-automatic SPPS was accomplished through the Biotage® Initiator™ synthesizer, assisted by microwave (MW) irradiation; Fmoc/tBu strategy and Rink Amide MHBA Resin (100–200 mesh; loading: 0.5 mmol/g) were used. Peptides were obtained with yields up to 30% and purities between 97% and >99%. Refer to the Supporting Information for a detailed description of materials, procedures, and methods.

Jasco J-715 (Tokyo, Japan) spectropolarimeter—equipped with a Haake thermostat (Thermo Fisher Scientific, Waltham, MA)—was used to collect circular dichroism spectra. Bidistilled water was used as solvent. Fused quartz cells of either 1 or 10 mm pathlength (Hellma, Mühlheim, Germany) were used. Spectra were baseline subtracted as expressed in terms of [θ]_T, total molar ellipticity (deg × cm² × dmol⁻¹).

NMR

2D-NMR experiments for conformational studies were carried out on a Bruker Avance DMX-600 instrument, using TOPSPIN 1.3 software package. For peptide 3, COSY (COrrelation SpectroscopY), CLEAN-TOCSY (TOtal Correlation SpectroscopY), and ROESY spectra (all with watergate for water suppression) were acquired. COSY was phase sensitive. TOCSY spectrum (spin-lock pulse, 70 ms) was acquired by collecting 400 recordings of 76 scans each. For the ROESY spectrum (200 ms mixing time), 512 experiments, each one consisting of 80 scans, were acquired. The full assignment was achieved exploiting the procedure proposed by K. Wüthrich (Wüthrich, 1986).

Inhibition of VEGF/VEGFR Complex Formation <italic>in vitro</italic>

The chemiluminescent screening assay for the detection of VEGFR-1 ligands was accomplished according to the procedure described in the Supporting Information. Unlabeled VEGF₁₆₅, employed as reference compound, showed an IC₅₀ value of 146 pm (Zanella et al., 2015), comparable with the previously reported value (Goncalves et al., 2007).

Inhibition of HUVEC Morphogenesis

To assess angiogenic activity, HUVECs (2.5 × 10⁴ cells) were seeded in a 24-well plate coated with 100 μL/well of Matrigel previously polymerized for 1 h at 37°C. Cells were then incubated for 5 h at 37°C in a moist atmosphere of 5% CO₂. The experiments were performed without or in the presence of 5 and 6 under either resting (absence of stimuli, cell cultured in EndoGRO medium alone, without FBS and all the growth factors) or stimulated conditions (addition of VEGF₁₆₅, 10 ng/mL). Network formation was evaluated by phase-contrast microscopy using a fluorescence microscope (AxioVert 40CFL, Carl Zeiss S.p.A. Milan, Italy). Five photos of each well were recorded with 10X magnification. Network formation was finally quantified in terms of mean number of loops per field as topological parameters and the total length of the branches. For the purpose of the analysis, loops were defined as any complete ring formed by HUVECs, while open ramifications were considered as branches. The analysis of the images was performed using the free software ImageJ (https://imagej.nih.gov/ij/).

Resistance to Proteolytic Degradation

The proteolytic stability of peptides 1, 5, and 6 was assessed against three enzymes: pronase, trypsin and chymotrypsin (Sigma-Aldrich). Each peptide was dissolved in the appropriate buffer [i.e., (2-amino-2-hydroxymethyl-propane-1,3-diol (Tris)·HCl 20 mm, containing 20 mm CaCl₂ pH 7.6 for pronase; Tris·HCl 50 mm, pH 7.8 for trypsin and chymotrypsin). Then, it was incubated with and without the enzyme solution (1.25 μg of enzyme in 150 mL buffer) at 37°C for 12 h. The mixture was analyzed by RP-HPLC (column: Phenomenex Jupiter 5 μ C18 300 Å) every ten minutes for the first hour, then every half hour. Gradients: 5–25%B in 10 min for peptide 1; 15–35%B in 10 min for peptide 5; 15–25%B in 10 min for peptide 6. Eluants: A, H₂O:CH₃CN 9/1 + 0.05% trifluoroacetic acid (TFA); B: CH₃CN:H₂O 9/1 + 0.05% TFA.

Data Availability

All datasets generated for this study are included in the manuscript and/or the Supplementary Files.

Author Contributions

SZ, LP, and GS synthesized the peptides. GB, together with SR, AP, and LS, performed the in silico studies. MD carried out the conformational studies. DA designed and performed the binding experiments and FM determined the biological activity of the synthesized compounds. AS performed proteolytic stability assays. LB and UP supervised the binding and activity studies, FF the structural characterization, CG the peptide synthesis, LS conceived the original idea and supervised peptide design. SZ, LB, and LS wrote the manuscript. All authors discussed the results and commented on the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The reviewer NI declared a past co-authorship and collaboration with one of the authors (LS).

We thank Dr. Monica Pinoli for technical assistance in the processing of the optical microscopy images of morphogenesis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2019.00170/full#supplementary-material

References Al Musaimi

Al Shaer

de la Torre

B. G.

Albericio

(2018). 2017 FDA peptide harvest. Pharmaceuticals 11:42. 10.3390/ph11020042

29735913

Baker

N. A.

Sept

Joseph

Holst

M. J.

McCammon

J. A.

(2001). Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl. Acad. Sci. U.S.A. 98, 10037–10041. 10.1073/pnas.181342398

11517324

Balsera

Bonache MÁ

Reille-Seroussi

Gagey-Eilstein

Vidal

González-Muñiz

. (2017). Disrupting VEGF–VEGFR1 interaction: de novo designed linear helical peptides to mimic the VEGF13-25 fragment. Molecules 22:1846. 10.3390/molecules22111846

29143774

Basile

Del Gatto

Diana

Di Stasi

Falco

Festa

. (2011). Characterization of a designed vascular endothelial growth factor receptor antagonist helical peptide with antiangiogenic activity in vivo. J. Med. Chem. 54, 1391–1400. 10.1021/jm101435r

21280635

Berendsen

H. J.

Postma

J. P.

van Gunsteren

W. F.

Hermans

(1981). Interaction models for water in relation to protein hydration, in Intermolecular Forces ed Pullman

(Dordrecht: Springer), 331–342. Berendsen

H. J.

Postma

J. V.

van Gunsteren

W. F.

DiNola

Haak

J. R.

(1984). Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81, 3684–3690. 10.1063/1.448118 Bonache

M. A.

Balsera

López-Méndez

Millet

Brancaccio

Gómez-Monterrey

. (2014). De novo designed library of linear helical peptides: an exploratory tool in the discovery of protein–protein interaction modulators. ACS Comb. Sci. 16, 250–258. 10.1021/co500005x

24725184

Bussi

Donadio

Parrinello

(2007). Canonical sampling through velocity rescaling. J. Chem. Phys. 126:014101. 10.1063/1.2408420 Carmeliet

Jain

R. K.

(2011). Principles and mechanisms of vessel normalization for cancer and other angiogenic diseases. Nat. Rev. Drug Discov. 10, 417–427. 10.1038/nrd3455 Chakrabartty

Kortemme

Baldwin

R. L.

(1994). Helix propensities of the amino acids measured in alanine-based peptides without helix-stabilizing side-chain interactions. Protein Sci. 3, 843–852. 10.1002/pro.5560030514

8061613

Chames

Van Regenmortel

Weiss

Baty

(2009). Therapeutic antibodies: successes, limitations and hopes for the future. Br. J. Pharmacol. 157, 220–233. 10.1111/j.1476-5381.2009.00190.x

19459844

Chen

J. C.

Chang

Y. W.

Hong

C. C.

Y. H.

J. L.

(2013). The role of the VEGF-C/VEGFRs axis in tumor progression and therapy. Int. J. Mol. Sci. 14, 88–107. 10.3390/ijms14010088

23344023

Cho

J. I.

Tanaka

Sato

Kinbara

Aida

(2010). Oligo (4-aminopiperidine-4-carboxylic acid): an unusual basic oligopeptide with an acid-induced helical conformation. J. Am. Chem. Soc. 132, 13176–13178. 10.1021/ja106118w

20812686

Dai

Castellino

F. J.

Prorok

(2004). A single amino acid replacement results in the Ca2+-induced self-assembly of a helical conantokin-based peptide. Biochemistry 43, 13225–13232. 10.1021/bi048796s

15476416

Dalzini

Bergamini

Biondi

De Zotti

Panighel

Fato

. (2016). The rational search for selective anticancer derivatives of the peptide Trichogin GA IV: a multi-technique biophysical approach. Sci. Rep. 6:24000. 10.1038/srep24000 Danhier

Le Breton

Préat

(2012). RGD-based strategies to target alpha(v) beta(3) integrin in cancer therapy and diagnosis. Mol. Pharm. 9, 2961–2973. 10.1021/mp3002733 De Rosa

Di Stasi

D'Andrea

L. D.

(2018). Pro-angiogenic peptides in biomedicine. Arch. Biochem. Biophys. 660, 72–86. 10.1016/j.abb.2018.10.010

30336134

De Zotti

Biondi

Formaggio

Toniolo

Stella

Park

. (2009). Trichogin GA IV: an antibacterial and protease-resistant peptide. J. Pept. Sci. 15, 615–619. 10.1002/psc.1135

19399781

Desgrosellier

J. S.

Cheresh

D. A.

(2010). Integrins in cancer: biological implications and therapeutic opportunities. Nat. Rev. Cancer 10, 9–22. 10.1038/nrc2748

20029421

Diana

Di Stasi

De Rosa

Isernia

D'Andrea

Fattorusso

(2013). Structural investigation of the VEGF receptor interaction with a helical antagonist peptide. J. Pept. Sci. 19, 214–219. 10.1002/psc.2480

23420665

Dolinsky

T. J.

Nielsen

J. E.

McCammon

J. A.

Baker

N. A.

(2004). PDB2PQR: an automated pipeline for the setup of Poisson–Boltzmann electrostatics calculations. Nucleic Acids Res. 32, W665–W667. 10.1093/nar/gkh381

15215472

Essmann

Perera

Berkowitz

M. L.

Darden

Lee

Pedersen

L. G.

(1995). A smooth particle mesh Ewald method. J. Chem. Phys. 103, 8577–8593. 10.1063/1.470117 Fanelli

Schembri

Piarulli

Pinoli

Rasini

Paolillo

. (2014). Effects of a novel cyclic RGD peptidomimetic on cell proliferation, migration and angiogenic activity in human endothelial cells. Vascular Cell 6:11. 10.1186/2045-824X-6-11

25053992

Ferrara

Adamis

A. P.

(2016). Ten years of anti-vascular endothelial growth factor therapy. Nat. Rev. Drug Discov. 15, 385–403. 10.1038/nrd.2015.17

26775688

Ferrara

Gerber

H. P.

LeCouter

(2003). The biology of VEGF and its receptors. Nat. Med. 9, 669–676. 10.1038/nm0603-669

12778165

Ferrara

Hillan

K. J.

Gerber

H. P.

Novotny

(2004). Discovery and development of bevacizumab, an anti-VEGF antibody for treating cancer. Nat. Rev. Drug Discov. 3, 391–400. 10.1038/nrd1381

15136787

Fogolari

Brigo

Molinari

(2003). Protocol for MM/PBSA molecular dynamics simulations of proteins. Biophys. J. 85, 159–166. 10.1016/S0006-3495(03)74462-2

12829472

García-Aranda

M. I.

González-López

Santiveri

C. M.

Gagey-Eilstein

Reille-Seroussi

Martín-Martínez

. (2013). Helical peptides from VEGF and Vammin hotspots for modulating the VEGF–VEGFR interaction. Org. Biomol. Chem. 11, 1896–1905. 10.1039/c3ob27312a

23381088

Gautier

Goncalves

Diana

Di Stasi

Teillet

Lenoir

. (2010). Biochemical and structural analysis of the binding determinants of a vascular endothelial growth factor receptor peptidic antagonist. J. Med. Chem. 53, 4428–4440. 10.1021/jm1002167

20462213

Goncalves

Gautier

Garbay

Vidal

Inguimbert

(2007). Development of a chemiluminescent screening assay for detection of vascular endothelial growth factor receptor 1 ligands. Anal. Biochem. 366, 108–110. 10.1016/j.ab.2007.03.027

17482136

Henninot

Collins

J. C.

Nuss

J. M.

(2018). The current state of peptide drug discovery: back to the future?. J. Med. Chem. 61, 1382–1414. 10.1021/acs.jmedchem.7b00318

28737935

Hess

Bekker

Berendsen

H. J.

Fraaije

J. G.

(1997). LINCS: a linear constraint solver for molecular simulations. J. Comp. Chem. 18, 1463–1472. 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-H Hess

Kutzner

Van Der Spoel

Lindahl

(2008). GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4, 435–447. 10.1021/ct700301q Hoeben

Landuyt

Highley

M. S.

Wildiers

van Oosterom

A. T.

De Bruijin

E. A.

(2004). Vascular endothelial growth factor and angiogenesis. Pharmacol. Rev. 56, 549–580. 10.1124/pr.56.4.3

15602010

Hopkins

S. A.

Konopelski

J. P.

Olmstead

M. M.

Banks

H. D.

(2000). Conformational analysis of peptides containing enantiomerically pure α-methylasparagine: correspondence between computed and solid state structures. Tetrahedron 56, 9733–9737. 10.1016/S0040-4020(00)00880-2 Hornak

Abel

Okur

Strockbine

Roitberg

Simmerling

(2006). Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins Struct. Funct. Gen. 65, 712–725. 10.1002/prot.21123

16981200

Hou

Wang

(2010). Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J. Chem. Inf. Model. 51, 69–82. 10.1021/ci100275a

21117705

Howard

D. H.

Bach

P. B.

Berndt

E. R.

Conti

R. M.

(2015). Pricing in the market for anticancer drugs. J. Econ. Perspect. 29, 139–162. 10.1257/jep.29.1.139

28441702

Johannessen

T. C.

Wagner

Straume

Bjerkvig

Eikesdal

H. P.

(2013). Tumor vasculature: the Achilles' heel of cancer? Expert Opin. Ther. Targets 17, 7–20. 10.1517/14728222.2013.730522

23121690

Kabsch

Sander

(1983). Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 22, 2577–2637. 10.1002/bip.360221211

6667333

Kelly

S. M.

Price

N. C.

(1997). The application of circular dichroism to studies of protein folding and unfolding. Biochim. Biophys. Acta Protein Struct. Mol. Enzymol. 1338, 161–185. 10.1016/S0167-4838(96)00190-2

9128135

Krüger

D. M.

Gohlke

(2010). DrugScore^PPI webserver: fast and accurate in silico alanine scanning for scoring protein–protein interactions. Nucleic Acids Res. 38, W480–W486. 10.1093/nar/gkq471

20511591

Leppänen

V. M.

Prota

A. E.

Jeltsch

Anisimov

Kalkkinen

Strandin

. (2010). Structural determinants of growth factor binding and specificity by VEGF receptor 2. Proc. Natl. Acad. Sci. U.S.A. 107, 2425–2430. 10.1073/pnas.0914318107

20145116

Leppänen

V. M.

Tvorogov

Kisko

Prota

A. E.

Jeltsch

Anisimov

. (2013). Structural and mechanistic insights into VEGF receptor 3 ligand binding and activation. Proc. Natl. Acad. Sci. U.S.A. 110, 12960–12965. 10.1073/pnas.1301415110

23878260

Fuh

Meng

Xin

Gerritsen

M. E.

Cunningham

. (2000). Receptor-selective variants of human vascular endothelial growth factor. J. Biol. Chem. 275, 29823–29828. 10.1074/jbc.M002015200 Michaelsen

S. R.

Staberg

Pedersen

Jensen

K. E.

Majewski

Broholm

. (2018). VEGF-C sustains VEGFR2 activation under bevacizumab therapy and promotes glioblastoma maintenance. Neuro Oncol. 20, 1462–1474. 10.1093/neuonc/noy103

29939339

Mizejewski

G. J.

(1999). Role of integrins in cancer: survey of expression patterns. Proc. Soc. Exp. Biol. Med. 222, 124–138. 10.1046/j.1525-1373.1999.d01-122.x

10564536

Segura Mora

Assi

S. A.

Fernandez-Fuentes

(2010). Presaging critical residues in protein interfaces-web server (PCRPi-W): a web server to chart hot spots in protein interfaces. PLoS ONE 5:e12352. 10.1371/journal.pone.0012352

20808810

Muller

Christinger

Wells

Cunningham

de Vos

(1997). Vascular endothelial growth factor: crystal structure and functional mapping of the kinase domain receptor binding site. Proc. Natl. Acad. Sci. U.S.A. 94, 7192–7197. 10.1073/pnas.94.14.7192 Musumeci

Radi

Brullo

Schenone

(2012). Vascular endothelial growth factor (VEGF) receptors: drugs and new inhibitors. J. Med. Chem. 55, 10797–10822. 10.1021/jm301085w

23098265

Nasir

(2019). Angiogenic signaling pathways and anti-angiogenic therapies in human cancer, in Predictive Biomarkers in Oncology eds Badve

Kumar

(Cham: Springer), 243–262 Niclou

S. P.

(2018). Revival of the VEGF ligand family? Neuro Oncol. 20, 1421–1422. 10.1093/neuonc/noy127

30137575

Olsson

A. K.

Dimberg

Kreuger

Claesson-Welsh

(2006). VEGF receptor signalling — in control of vascular function. Nat. Rev. Mol. Cell Biol. 7, 359–371. 10.1038/nrm1911 Onufriev

Bashford

Case

D. A.

(2004). Exploring protein native states and large-scale conformational changes with a modified generalized born model. Proteins Struct. Funct. Genet. 55, 383–394. 10.1002/prot.20033

15048829

Oostenbrink

Villa

Mark

A. E.

Van Gunsteren

W. F.

(2004). A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6. J. Comput. Chem. 25, 1656–1676. 10.1002/jcc.20090

15264259

Pace

C. N.

Scholtz

J. M.

(1998). A helix propensity scale based on experimental studies of peptides and proteins. Biophys. J. 75, 422–427. 10.1016/S0006-3495(98)77529-0 Patriksson

van der Spoel

(2008). A temperature predictor for parallel tempering simulations. Chem. Chem. Phys. 10, 2073–2077. 10.1039/b716554d Pettersen

E. F.

Goddard

T. D.

Huang

C. C.

Couch

G. S.

Greenblatt

D. M.

Meng

E. C.

. (2004). UCSF chimera—a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612. 10.1002/jcc.20084 Robinson

C. J.

Stringer

S. E.

(2001). The splice variants of vascular endothelial growth factor (VEGF) and their receptors. J. Cell Sci. 114, 853–865. Ronca

Benkheil

Mitola

Struyf

Liekens

(2017). Tumor angiogenesis revisited: regulators and clinical implications. Med. Res. Rev. 37, 1231–1274. 10.1002/med.21452

28643862

Ryckaert

J. P.

Ciccotti

Berendsen

H. J.

(1977). Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J. Comp. Phys. 23, 327–341. 10.1016/0021-9991(77)90098-5 Schymkowitz

Borg

Stricher

Nys

Rousseau

Serrano

(2005). The FoldX web server: an online force field. Nucleic Acids Res. 33, W382–W388. 10.1093/nar/gki387

15980494

Shibuya

(2013). Vascular endothelial growth factor and its receptor system: physiological functions in angiogenesis and pathological roles in various diseases. J. Biochem. 153, 13–19. 10.1093/jb/mvs136 Sidman

R. L.

Lawrence

Musso

G. F.

Giordano

R. J.

. (2015). The peptidomimetic Vasotide targets two retinal VEGF receptors and reduces pathological angiogenesis in murine and nonhuman primate models of retinal disease. Sci. Trans. Med. 7:309ra165. 10.1126/scitranslmed.aac4882

26468327

Siemeister

Marmé

Martiny-Baron

(1998). The alpha-helical domain near the amino terminus is essential for dimerization of vascular endothelial growth factor. J. Biol. Chem. 273, 11115–11120. 10.1074/jbc.273.18.11115

9556597

Simon

Gagliano

Giamas

(2017). Direct effects of anti-angiogenic therapies on tumor cells: VEGF signaling. Trends Mol. Med. 23, 282–292. 10.1016/j.molmed.2017.01.002

28162910

Somanath

P. R.

Malinin

N. L.

Byzova

T. V.

(2009). Cooperation between integrin alphavbeta3 and VEGFR2 in angiogenesis. Angiogenesis 12, 177–185. 10.1007/s10456-009-9141-9

19267251

J. L.

Yen

C. J.

Chen

P. S.

Chuang

S. E.

Hong

C. C.

Kuo

I. H.

. (2007). The role of the VEGF-C/VEGFR-3 axis in cancer progression. Br. J. Cancer 96, 541–545. 10.1038/sj.bjc.6603487

17164762

Toniolo

Crisma

Formaggio

Peggion

(2001). Control of peptide conformation by the Thorpe-Ingold effect (Cα-tetrasubstitution). Pept. Sci. 60, 396–419. 10.1002/1097-0282(2001)60:6<396::AID-BIP10184>3.0.CO;2-7 Udugamasooriya

D. G.

Dineen

S. P.

Brekken

R. A.

Kodadek

(2008). A peptoid “Antibody Surrogate” that antagonizes VEGF receptor 2 activity. J. Am. Chem. Soc. 130, 5744–5752. 10.1021/ja711193x.

18386897

Vasudev

N. S.

Reynolds

A. R.

(2014). Anti-angiogenic therapy for cancer: current progress, unresolved questions and future directions. Angiogenesis 17, 471–494. 10.1007/s10456-014-9420-y

24482243

Wang

C. A.

Tsai

S. J.

(2015). The non-canonical role of vascular endothelial growth factor-C axis in cancer progression. Exp. Biol. Med. 240, 718–724. 10.1177/1535370215583802 Wüthrich

(1986). NMR of Proteins and Nucleic Acids. New York, NY: Wiley. Zanella

Mingozzi

Dal Corso

Fanelli

Arosio

Cosentino

. (2015). Synthesis, characterization, and biological evaluation of a dual-action ligand targeting α_vβ₃ integrin and VEGF receptors. ChemistryOpen 4, 633–641. 10.1002/open.201500062

26491644

Zhu

Mitchell

J. C.

(2011). KFC2: a knowledge-based hot spot prediction method based on interface solvation, atomic density, and plasticity features. Proteins Struct. Funct. Bioinf. 79, 2671–2683. 10.1002/prot.23094

21735484

Zilberberg

Shinkaruk

Lequin

Rousseau

Hagedorn

Costa

. (2003). Structure and inhibitory effects on angiogenesis and tumor development of a new vascular endothelial growth inhibitor. J. Biol. Chem. 278, 35564–35573. 10.1074/jbc.M304435200

Funding. This project was supported by the Italian Ministry of Education, Universities and Research, grant PRIN 20157WW5EH. LS, GB, and AP acknowledge support from AIRC, IG2016 19171 and from PRACE (grant 2017174118), which awarded computational resources at CINECA, Italy.