Targeting transdiﬀerentiated hepatic stellate cells and monitoring the hepatic fibrogenic process by means of IGF2R-specific peptides designed in silico †

The lack of accurate and easily applicable methods for the diagnosis of liver fibrosis, a disease characterized by an accumulation of the extracellular matrix released by activated hepatic stellate cells (HSCs), has been a major limitation for the clinical management of liver diseases. The identification of biomarkers specific to liver microstructure alterations, combined with a non-invasive optical imaging modality, could guide clinicians towards a therapeutic strategy. In this study, structural information of the insulin-like growth factor 2 receptor (IGF2R), an overexpressed protein on activated HSCs, was used for in silico screening of novel IGF2R-specific peptide ligands. Molecular dynamics simulations, followed by computational alanine scanning of the IGF2R/IGF2 complex, led to the identification of a putative peptide sequence containing the most relevant amino acids for the receptor–ligand interaction (IGF2 E12-C21). The Residue Scan tool, implemented in the MOE software, was then used to optimize the binding affinity of this sequence by amino acid mutations. The designed peptides and their associated scrambled sequences were fluorescently labelled and their binding affinity to LX-2 cells (model for activated human HSCs) was tested using flow cytometry and confocal microscopy. In vitro binding was verified for all sequences ( K D r 13.2 m M). With respect to the putative binding sequence, most mutations led to an increased affinity. All sequences have shown superior binding compared to their associated scrambled sequences. Using HPLC, all peptides were tested in vitro for their proteolytic resistance and showed a stability of Z 60% intact after 24 h at 37 1 C in 50% v/v FBS. In view of their prospective diagnostic application, a comparison of binding affinity was performed in perpetuated and quiescent-like LX-2 cells. Furthermore, the IGF2R expression for different cell phenotypes was analysed by a quantitative mass spectrometric approach. Our peptides showed increased binding to the perpetuated cell state, indicating their good selectivity for the diagnostically relevant phenotype. In summary, the increased binding affinity of our peptides towards perpetuated LX-2 cells, as well as the satisfactory proteolytic stability, proves that the in silico designed sequences offer a new potential strategy for the targeting of hepatic fibrosis.


Introduction
Fibrosis is a progressive dysregulated pathway of the natural wound healing response, defined as an excessive accumulation and deposition of the extracellular matrix (ECM) released by activated myofibroblasts after acute tissue injury of any origin. 1he replacement of healthy functional tissue with collagenproducing cells forming scar tissue and the resulting architectural deformation lead to fibrosis, cirrhosis, and hepatocellular carcinoma in its final and pathological manifestation. 2 According to recent epidemiological studies, 10 million European citizens with undiagnosed liver disease are bearing significant liver fibrosis and approximately 1.5 million have silent cirrhosis. 3At present, a definite diagnosis of early or severe fibrosis can only be achieved by means of liver biopsy, an invasive and painful methodology that often has a high risk for complications and misjudgement. 4ence, to improve prognosis for the patients, it is of paramount importance to develop a more convenient method for comprehensive and early diagnosis of hepatic fibrosis. 5As a first assessment, the ratios of combinatorial markers like aspartate transaminase/alanine aminotransferase (AST/ALT) and the ASTto-platelet ratio index (APRI), or various direct and indirect markers, such as collagens and collagenases, glycoproteins, or cytokines released by inflamed tissue, are only preliminary indicators monitored prior to proceeding with further diagnoses. 6hese serum markers are mostly unspecific to liver tissue or liver fibrosis and do not provide any information concerning the grade and stage of the disease. 5,6Another approach to determine fibrosis-related transformations in hepatic tissue is to use clinical imaging techniques such as ultrasound and contrast-enhanced ultrasonography, which enable the alteration of the liver surface and of the intrahepatic microcirculation to be evaluated by assessing changes in the hepatic architecture with a satisfactory specificity (95%) for severe liver fibrosis. 6However, due to inadequate sensitivity and high inter-observer variability, these techniques offer limited potential for the evaluation of early fibrosis progression, staging, and grading. 7FibroScan TM s is a widely used, well-validated, and user-friendly elastographic technique that can be performed at the bedside with moderate reproducibility.As a function of increased liver stiffness, which occurs in hepatic fibrosis, the transfer of shear waves will be attenuated. 8The applicability depends on the experience of the operator and on the physical constitution of the patient, bringing along substantial limitations for obese patients suffering from ascites or a low diaphragm. 8Magnetic resonance imaging has the potential to become a reference technique to determine the stage, grade, and progress of liver fibrosis, but it suffers from a lack of appropriate and selective contrast agents, especially at early stage fibrosis. 4ver the last decade, the development of peptide-based contrast agents for optical medical imaging of acute and chronic inflammations and their underlying pathways has become increasingly prominent. 9The peptide contrast agents are produced by linking the peptide sequence directly to a fluorescent probe, to magnetic resonance imaging and photo-acoustically active agent, or onto the surface of nanoscale lipid carrier systems incorporating contrast agents as an additional payload. 10,11Direct binding of these labelled peptides to biomarkers overexpressed on the surface of a disease-specific cell type would result in sitespecific labelling and a reduction in background noise. 12n liver fibrosis, hepatic stellate cells (HSCs) are central disease mediators and thus considered optimal targets. 13HSCs undergo a mesenchymal transdifferentiation during fibrosis progression from a quiescent cell state with vitamin A storing capacity to an activated myofibroblast-like phenotype. 13,14on activation, it has been observed that the insulin growth factor 2 receptor (IGF2R), also known as the cationindependent mannose-6-phosphate receptor, is selectively overexpressed on the HSC surface. 15][18] Owing to the ability of the receptor to internalize rather rapidly extracellular molecules and its ubiquitous expression on activated HSCs, IGF2R represents an optimal candidate for the development of a functional diagnostic platform for liver fibrosis. 17,19By means of a combinatorial bio-panning strategy on immobilized IGF2R and rat HSCs, a dodecamer peptide (VHWDFRQWWQPS, P431) was identified by Chen and co-workers as the best ligand for IGRF2R and was validated for its specificity to HSCs using LX-2 cells, an immortalized human HSC line. 20In the present study, the structural information of IGF2R derived from the PDB database was used for an in silico screening of novel peptide ligands specific to this receptor.Subsequently, the sequences showing the best in silico binding affinity for IGF2R were chosen to assess their selectivity to target activated HSCs over quiescentlike HSCs in comparison to peptide P431.The selected peptides and their scrambled sequences were labelled with fluorescein isothiocyanate (FITC) (Table 1).Their binding to different IGF2R-expressing human cell lines as well as cell uptake was assessed using flow cytometry and confocal microscopy.Furthermore, the expression of IGF2R was analysed for all tested cell lines by means of a quantitative mass spectrometric approach.Unlabelled sequences were tested in vitro for their proteolytic resistance in view of their prospective diagnostic application.

Methods
Model parameterization and simulation protocol.A molecular model of IGF2R (domains 11-13) bound with IGF2 was taken from the PDB database (PDB code: 2V5P). 17Cysteine disulphide bonds were manually imposed according to the PDB file.The ff14SB force field 21 was chosen for the protein complex; protonation states were assigned to mimic conditions at physiological pH (pH 7.4).The structure was solvated with about 39 000 explicit TIP3P 22 water molecules and 14 Na + ions were added to assure electroneutrality; ion parameters optimized for the TIP3P water model were taken from Joung and Cheatham. 23,24lecular dynamics (MD) simulations were carried out as follows: first, energy minimization was performed applying a harmonic restraint on the solute (harmonic constant equal to 500 kcal mol À1 Å À2 ) in order to remove bad solute/solvent and solvent/solvent contacts due to the random placement of explicit water molecules; energy minimization was subsequently repeated on the overall system, removing the restraint.The temperature was raised to 310 K by means of 20 ps in an NVT ensemble (i.e., at a constant number of particles N, volume V, and temperature T) and system density was equilibrated through 1 ns in NPT ensemble (i.e., at a constant number of particles N, pressure P, and temperature T) at 310 K and 1 atm.Finally, 100 ns of MD simulations were carried out in NPT ensemble at 310 K and 1 atm.Temperature and pressure were kept at desired values using Langevin dynamics (with a collision frequency equal to 1 ps À1 ) and a Berendsen barostat with isotropic position scaling, respectively.
Simulations were performed adopting periodic boundary conditions; a cut-off value equal to 10 Å was chosen for nonbonded interactions.The SHAKE 25 algorithm was employed to restrain all covalent bonds involving hydrogen atoms: this allowed adopting a time step equal to 2 fs to propagate dynamics by integrating Newton's equation of motion using the Leap-Frog algorithm.Molecular coordinates were collected every 20 ps.The attainment of an equilibrated complex structure was verified by evaluating root mean square displacement (RMSD) and interaction energy over time by means of Molecular Mechanics Generalized Born Surface Area (MMGBSA, vide infra).Plots are presented in the ESI † (Fig. S1 and S2).
Table 1 Peptide sequences as well as HPLC and MS analysis results.Tc is a decamer representing the putative natural peptide sequence containing the most relevant amino acids for the interaction of IGF2 with IGF2R.Mutations in the Tc sequence resulted in peptides T1-T5.Amino acids that are in bold and underlined indicate mutated sites.Scrambled sequences are labelled SCR.Representative chromatograms and representative mass spectra are available in the ESI (Fig. S3A-D Molecular Mechanics Generalized Born Surface Area.MMGBSA 26 was employed to compute the interaction energy DE int between IGF2 and IGF2R: where DE complex , DE receptor and DE ligand are the energy change values related to the complex, the receptor (IGF2R) and the ligand (IGF2), respectively.According to MMGBSA formalism, the energy change of each component DE is the sum of two contributions: where E gas is the gas phase energy computed with the force field, which accounts for both bonded (bonds, angles, dihedrals) and non-bonded (van der Waals, electrostatic) interactions, while DG solv is the solvation free energy, which is, in turn, the sum of two terms: where DG polar and DG np are the polar and the non-polar contribution to solvation free energy, respectively.The first term is computed through an implicit solvent method; in this framework, generalized Born method 27 was chosen for the sake of computational efficiency, adopting a salt concentration equal to 0.137 M, in order to account for the physiological environment.The second term is evaluated as a function of the solvent accessible surface area (SASA): where a is a constant, equal to 0.0072 kcal mol À1 Å À2 .Per-residue decomposition was subsequently performed in order to highlight the specific contribution to the interaction energy of IGF2 amino acids involved in the binding, aiming at highlighting the most relevant ones that lead to the formation of the non-covalent complex. 26The role of the amino acids identified through perresidue decomposition was further investigated by means of computational alanine scanning (CAS). 26CAS is defined as the evaluation of the variation in binding free energy (DDG b ) caused by the implementation of mutations to Ala of the selected interfacial residue (amino acid). 28Briefly, a single given amino acid is mutated to alanine and the interaction energy is recomputed using the original molecular trajectory, adopting the MMGBSA formalism discussed above.The resulting change in interaction energy provides a better quantification of the importance of the mutated amino acid in the attainment of the complex.The obtained results were the starting point for the design of peptides (vide infra).Calculations were performed considering the last 20 ns of the MD trajectory, averaging over 1000 frames.Design of peptides.The last frame of the MD trajectory described above was chosen as the starting point to design inhibitor peptides.The MOE software was used to perform all the calculations described hereafter. 29The starting peptide was obtained from IGF2 (chain D of 2V5P PDB file) by deleting all residues, except the helical sequence spanning from E12 to residue C21.Indeed, this sequence contained the highest number of contiguous hot and warm spots, according to CAS results.In the IGF2 structure C21 forms a disulphide bridge with C60 and does not interact with the receptor.For this reason, it was replaced by a serine to increase the hydrophilic character of the resulting peptide.The peptide was then capped at the N-and C-termini with acetyl and N-methyl groups, respectively.The obtained complex, formed by IGF2R (chain A of 2V5P PDB file) and the E12-S21 peptide (hereafter referred to as WT), was subjected to geometry minimization using the Amber10:EHT force field and the Reaction Field solvation model for water, implemented in MOE. 29Minimization was initially performed on the peptide alone, keeping the receptor restrained.A second minimization step included the receptor, but kept the backbone restrained.The last step consisted in an unrestrained minimization of the whole system up to a gradient of 0.1 kcal mol À1 Å À2 .The model was then used for the design of optimized peptides via the Residue Scan (RS) tool implemented in MOE.RS allows the evaluation of the differences in binding energy (dAffinity, DG bind ) and in thermodynamic stability (dStability, DG tot ) following the mutation of a given residue.In the chosen setup, DG bind and DG tot were evaluated on an ensemble of low-energy conformations generated, after the mutation, using the Low Mode molecular dynamics method implemented in MOE. 29ITC labelling.Peptide P431 and peptide SCR431 were labelled with fluorescein isothiocyanate (FITC).First, a 4 mg mL À1 stock solution of FITC in anhydrous DMSO was freshly prepared and kept in the dark.Peptide stock solutions were prepared at a concentration of 5 mg mL À1 in anhydrous DMSO.Subsequently, 100 mL peptide stock solution and 120 mL FITC stock solution were added to 280 mL borate buffer (pH 9) to obtain a final peptide concentration of 1 mg mL À1 and a FITC-to-peptide molar ratio of 5 : 1.The reaction mixture was stirred at room temperature overnight in the dark and purified by preparative high performance liquid chromatography (HPLC, vide infra).The purity was confirmed by analytical HPLC (vide infra) and the compounds were characterized by mass spectrometry (LC-MS, vide infra).
Peptide purification.For the purification of P431_FITC and SCR431_FITC, a HPLC (1260 Infinity II, Agilent Technologies, Santa Clara, CA, USA or Ultimate 3000, Thermo Fisher, Ma ¨nnedorf, Switzerland) equipped with a quaternary pump, an autosampler, and a DAD-UV detector was used.A MN Hypersil column (C18, 4.0 Â 250 mm, 5 mm; Macherey Nagel, Du ¨ren, Germany) was kept at 30 1C.The injection volume was 100 mL and the flow rate was set to 1.0 mL min À1 .The elution of the peptides started at t = 0 min with 75% v/v water (+0.1% v/v TFA) and 25% v/v ACN (+0.1% v/v TFA), followed by a linear gradient of ACN (+0.1% v/v TFA): 25-50% v/v over 25 min, 50-95% v/v over 0.01 min.The isolated product was then lyophilized for further use.
Mass spectrometry.Peptides (100 nM) in a mixture of ACN and water (1 : 5 v/v) + 0.07% formic acid were analyzed using a QTRAP 5500 mass spectrometer (Sciex, Darmstadt, Germany) equipped with an electrospray ionization source and operating in the positive ion mode.For full scans (m/z 600-900), the ion spray voltage was set to 3000 V, the sheath gas pressure to 16 psi, and the declustering potential to 50 V.In contrast to the settings described above, for product ion scans (m/z 50-1200), a sheath gas pressure of 13 psi and a declustering potential of 150 V were used, and the collision energy was ramped from 10 to 60 eV in 2 eV steps.Mass spectra were processed using Analyst 1.6 (Sciex, Darmstadt, Germany).
Purity.The purity of peptides was evaluated by analytical HPLC (1260 Infinity II, Agilent Technology, Santa Clara, CA, USA).A MN Hypersil column (C18, 4.0 Â 250 mm, 5 mm; Macherey-Nagel, Du ¨ren, Germany) was used and kept at 60 1C.The flow rate was set to 1.0 mL min À1 and the absorbance was recorded at l = 220 nm.The injection volume was set to 20 mL.At the beginning, the mobile phase contained 25% v/v of mobile phase A (ACN + 0.1% v/v TFA) and 75% v/v of mobile phase B (ultrapure water + 0.1% v/v TFA).After 16 min mobile phase A was increased to 90% using a linear gradient, followed by a slight increase to 95% at 18.01 min.
Serum stability.Peptides (25 mM) were incubated with 50% v/v FBS in 500 mL PBS (pH 7.4) for 24 h at 37 1C.Aliquots (100 mL) were sampled after 0.5, 1, 4, 8, and 24 h and inactivated with 200 mL ice cold methanol containing 25 mg mL À1 ferulic acid as the internal standard, followed by 30 s of vortexing and 20 min at 4 1C.Inactivated enzymes were removed by centrifugation for 15 min at 12 000 Â g in a MiniSpin centrifuge (Eppendorf, Hamburg, Germany).The supernatant was collected.The analysis was performed using an Agilent 1260 Infinity II HPLC.The injection volume of the purified peptides was 20 mL and the temperature of the autosampler was adjusted to room temperature (RT).For the analysis, a previously described method was applied (peptide purification).Peptide T4 was analysed differently due to an overlap with an FBS-derived peak in the chromatogram.Briefly, the analysis was started with 40% v/v of mobile phase A and 60% v/v of mobile phase B, which were kept constant for 25 min, followed by a rapid increase to 95% v/v of mobile phase A. After 28 min, the mobile phase contained 25% v/v of mobile phase A. These parameters were kept constant until the end of the run at 35 min.In order to determine the amount of intact peptide, the following equations were used to calculate the percentages at certain time points: where C t represents the ratio of the area under the curve (AUC 220 nm ) of the sample divided by the AUC 220 nm of the internal standard observed at different time points and C 0 represents the ratio of the AUC 220 nm of the sample divided by the AUC 220 nm of the internal standard at time zero.Cell culture.LX-2 cells were kept in DMEM (containing 4.5 g L À1 glucose) supplemented with 1% v/v L-glutamine (LGlu), 2% v/v FBS, and 1% v/v penicillin/streptomycin solution (P/S).HepG2 cells were maintained in RPMI supplemented with 1% v/v LGlu, 10% v/v FBS and 1% v/v P/S.HEK293 cells were grown in DMEM supplemented with 1% v/v LGlu, 10% v/v FBS and 1% v/v P/S.All cells were maintained at 37 1C in a 5% CO 2 -humidified atmosphere.Subcultivation was performed at a cell confluence of about 90%.The initiation of the perpetuated or quiescent-like state of the LX-2 cells was performed following the protocol published by Valentino et al. 30 Briefly, to obtain quiescent LX-2 cells, a solution of retinol (Rol) and palmitic acid (PA) in DMSO was intensively mixed with FBS-free DMEM containing LGlu and P/S to reach a final concentration of 10 mM Rol and 300 mM PA.For the perpetuation of LX-2 cells, TGF-b 1 dissolved in PBS was diluted in FBS-free DMEM containing LGlu and P/S to reach a final concentration of 10 ng mL À1 .Cells were then incubated with these solutions for 24 h at 37 1C in a 5% CO 2 -humidified atmosphere.After the incubation, the cells were prepared for experiments as described for naive LX-2 cells.For all experiments, cells from three different and independent subcultures with the same cell passage numbers were used.LX-2 cells were used at passage numbers P11 and P12.HepG2 cells were used at passage number P7 and HEK293 were used at passage number P6.Cells were regularly checked for the absence of mycoplasma.
Cellular metabolism assay.MTT was dissolved in PBS to obtain a final concentration of 5 mg mL À1 .LX-2 cells were seeded into 96-well plates with 100 mL per well culturing medium (vide supra) at a density of 1.2 Â 10 4 cells per well.After 18 h at 37 1C under 5% CO 2 , cells were rinsed twice with PBS and then incubated with 100 mL per well unlabelled peptide dissolved in phenol red-free cell culture medium (5, 10 mM) for 8 h at 37 1C under 5% CO 2 .MTT dissolved in PBS (5 mg mL À1 ) was then added (10 mL per well) and the plates were left for an additional 4 h at 37 1C under 5% CO 2 .Afterwards 100 mL of anhydrous DMSO was added to each well to solubilize the water-insoluble formazan crystals.Absorbance was then measured using a plate reader (Spark 10M, Tecan, Ma ¨nnedorf, Switzerland) at l = 540 nm at room temperature.The following equation was used to calculate the cell metabolic activity in percent: Cell metabolic activity (%) = (OD sample /OD control ) Â 100 (7)   where OD sample represents the optical density of the cells treated with the peptides at different concentrations and OD control is the optical density of the cells exposed to DMEM.
Cell binding and uptake of labelled peptides in different cell lines and determination of the equilibrium dissociation constant (K D ) for LX-2 cells.After subcultivation by means of non-enzymatic cell dissociation solution, 2 Â 10 5 cells resuspended in culture medium were pipetted into a reaction tube and washed twice by centrifugation with 500 mL DPBS without calcium and magnesium, supplemented with 0.02% w/v EDTA and 2% v/v FBS (FACS buffer).Cells were then centrifuged at 300 Â g for 5 min at RT.The supernatant was discarded and 500 mL of the peptides in FACS buffer (10 mM) were added to each tube and the cells were thoroughly resuspended upon gentle vortexing.For the evaluation of K D , LX-2 cells were treated with different concentrations of FITC-labelled peptides (1, 5, 10, 25, and 50 mM), followed by 1 h incubation in the dark under light shaking at 37 1C.Cells were then centrifuged (300 Â g, 4 1C, 5 min).All supernatants were discarded, and the cell pellets were resuspended and washed three times by centrifugation with 500 mL ice-cold FACS buffer.Subsequently, 350 mL ice-cold FACS buffer was used for the final resuspension.Cells were kept in the dark at 4 1C prior to flow cytometry analysis (vide infra).To determine the percentage of FITC-positive cells, cells treated with DMSO in FACS buffer were used as the mock control.
Flow cytometry.A volume of 350 mL cell dispersion in FACS buffer (vide supra) was filtered through a 70 mm mesh to exclude agglomerated cells.After the filtration, 2 mL of 50 mg mL À1 propidium iodide (PI) solution in Milli-Q water was added to discriminate between living and dead cells.Pipetting gently, the dispersion was transferred into a 5 mL conical tube.The data were acquired using a FACSCanto TM Plus flow cytometer (BD Biosciences, Franklin Lakes, USA) and a FACS LSR II SORP (BD Biosciences, Franklin Lakes, USA).A minimum of 20 000 events in total per sample were recorded.A 488 nm laser was used for fluorescence excitation.The FITC channel (530/30 nm band-pass) was used to acquire FITC emission on the peptides and the PerCP channel (695/40 nm band-pass) was used to detect PI emission.The acquisition software was BD FACS Diva 8.0.1 (BD Biosciences, Franklin Lakes, USA).The data were analysed using FlowJo V10.0 software (BD Biosciences, Franklin Lakes, USA).The fluorescence intensity of the peptides bound to the cell surface and/or internalized by the cells was assessed by means of a gating strategy, as displayed in Fig. 1, in order to ensure that unwanted events and artefacts were minimized.Cellspecific signals were used for doublet discrimination.FSC-W/FSC-H was used to detect double positive signals.
It is well documented that cells with compromised membrane integrity, such as apoptotic cells, are prone to uncontrolled and unspecific uptake of proteins and peptides, with a resultant higher fluorescence emission due to the increased uptake. 31Consequently, an exclusion of dead cells from the final population was performed using propidium iodide.After excluding all unwanted events that could be responsible for artefacts, the percentage of FITC-positive labelled cells was assessed by gating all cells exceeding a certain value.To determine the threshold, unstained cells were used, and the threshold was set at 0.25% positive events in the control.In every set of experiments, the gating strategy and the compensation were separately adjusted optimizing the signal for each cell type to exclude cell type-specific differences in auto-fluorescence and complexity that could lead to shifts in the detected signal.
Confocal microscopy.LX-2 cells were seeded into a m-slide 8-well ibiTreat (ibidi GmbH, Gra ¨felfing, Germany) chambered coverslip in 300 mL culturing medium (vide supra) at a density of 2.5 Â 10 3 cells per chamber.After 18 h at 37 1C under 5% CO 2 , cells were washed twice with PBS and then incubated with 200 mL per well FITC-labelled peptides and their scrambled sequences (10 mM) dissolved in PBS with 2% v/v FBS for 1 h at 37 1C under 5% CO 2 .After incubation, the supernatant was removed and cells were gently washed 3 times with ice-cold PBS.Cells were then fixed with 200 mL per chamber ROTI s Histofix 4% for 10 min at RT in the dark.After fixation and two washing steps with PBS, cell nuclei were stained using 200 mL DAPI dissolved in ultrapure water (300 nM) for 3 min.After additional two washing steps, cells were transferred to a Zeiss LSM 800 confocal microscope (Carl Zeiss AG, Jena, Germany).For image acquisition, a 401Â/1.3oil objective with DAPI and GFP Filter was used.The acquisition software was Zeiss ZEN 2.6 (blue edition).Images were analysed using ImageJ (Laboratory for Optical and Computational Instrumentation, University of Wisconsin, Bethesda, USA) and the mean of the total corrected cell fluorescence (CTCF) was determined.
Detection of IGF2R by a quantitative mass spectrometric approach.LX-2 cells were kept and activated in DMEM (containing 4.5 g L À1 glucose) supplemented with 1% v/v LGlu, 2% v/v FBS, and 1% v/v penicillin/streptomycin solution (P/S).HepG2 cells were maintained in RPMI supplemented with 1% v/v LGlu, 10% v/v FBS and 1% v/v P/S.HEK293 cells were grown in DMEM supplemented with 1% v/v LGlu, 10% v/v FBS and 1% v/v P/S.All cells were maintained at 37 1C in a 5% CO 2 -humidified atmosphere.Subcultivation was performed at a cell confluence of about 90%.For the experiments, cells from three different and independent subcultures with the same cell passage number were used.Per independent replicate, 1 Â 10 7 cells were analysed.For the final analysis via mass spectrometry, plasma membrane proteins were isolated according to the protocol by Lin et al. 32 The obtained cell pellets were re-suspended (1 : 30 w/v) in hypotonic buffer (50 mM mannitol and 5 mM HEPES, pH 7.4, supplemented with protease inhibitor) and left in an ultrasound bath for 1 min (in ice).Subsequently, 1 M calcium chloride solution was added to a final concentration of 10 mM and the suspension was vortexed for 10 min at RT.After centrifugation for 15 min at 3000 Â g (4 1C), the supernatant was transferred into a new tube and ultra-centrifuged at 48 000 Â g (4 1C, 30 min).The pellet was suspended in 8 M urea/100 mM Tris (pH = 8), and a BCA assay was performed to quantify the absolute amount of proteins.Subsequently, 10 mg of protein was reduced, alkylated, and digested by Lys C for 2 h and then by trypsin overnight.The digests were analyzed by liquid chromatography (LC)-MS/MS (PROXEON coupled to a Q Exactive HF mass spectrometer, Thermo Fisher Scientific, Ma ¨nnedorf, Switzerland) with one injection of 5 mL digests.Peptides were trapped on a mPrecolumn PepMap100 C18 (5 mm, 100 Å, 300 mm Â 5 mm, Thermo Fisher Scientific, Ma ¨nnedorf, Switzerland) and separated by backflush on a C18 column (5 mm, 100 Å, 75 mm Â 15 cm) by applying a 60 minute gradient of 5-40% v/v ACN (+0.1% v/v formic acid) in ultrapure water (+0.1% v/v formic acid) at a flow rate of 350 nL min À1 .A parallel reaction monitoring approach was used with the inclusion list of MPRI human peptides, the resolution being set at 30 000, with an automatic gain control target of 2E05, a maximum ion injection time of 130 ms, HCD collision energy to 27, a fill ratio of 1%, charge exclusion of unassigned and 1+ ions, and peptide match preferred, respectively.The results were analyzed using Skyline v20.1.0.155 software (MacCoss Lab Software, University of Washington, Seattle, USA) at the MS2 fragment level.
Statistical analysis.Results are expressed as mean AE S.D. (n = 3).If not stated otherwise, differences between groups were determined with an ordinary one-way ANOVA combined with a multiple comparisons test after performing a normality test.Results are considered statistically significant if p r 0.05 (*), p r 0.01 (**), p r 0.001 (***) and p r 0.0001 (****).For results based on calculations with experimental triplicates, the propagated uncertainty of error is displayed.

Results and discussion
Molecular dynamics simulations of the IGF2/IGF2R complex A 100 ns MD simulation was performed in order to obtain a reasonably equilibrated structure of the complex; this was successfully verified by computing the RMSD and complex energetics as a function of simulation time (Fig. S1 and S2, ESI †).Interaction energy was obtained through the MMGBSA approach considering the last 20 ns of the molecular trajectory; this led to a DE int value equal to À48.54 AE 7.58 kcal mol À1 .Specific contributions are listed in Table 2 where DE elec and DE vdW are the contributions to the binding due to electrostatic and van der Waals interactions, respectively, and DDG solv is the change in solvation free energy after the attainment of the noncovalent complex.
Trajectory post-processing revealed that both electrostatic (which include salt bridges) and van der Waals (that also account for hydrophobic effects) interactions contribute to the formation of the complex.The positive value of DDG solv indicates that the complex attainment slightly penalizes the overall solubility; this is due to the charged moieties that are no more exposed to the solvent when the complex is formed.Per-residue decomposition allowed highlighting the relevant contribution of the fragment E21-S21 of IGF2, which is almost entirely included in the binding pocket; detailed outcomes are reported in Table S1 (ESI †).
Alanine scanning was subsequently performed with the selected amino acids from the previous analysis, in order to underline the impact of their mutation.The results are summarized in Table 3.
As discussed in detail below (vide infra), alanine scanning highlighted two hot spots for the binding: E12 (that leads to a salt bridge with R1571) and F19, whose contribution is related to hydrophobic effects.Although the analysis also identified E6 as relevant for the binding, the E12-S21 segment is better enclosed in the binding pocket and thus constitutes an ideal starting point for the in silico design of new peptides.
In silico design of peptides Fig. 2 presents the WT peptide (sequence E12-S21) bound to IGF2R.According to CAS analysis, two principal hot-spots mediate the interaction: E12, which makes a salt-bridge with R1571, and F19, whose benzyl group makes p-p interactions with F1567 and hydrophobic contacts with the pocket formed by Y1542, F1567, Y1606, L1626, L1629 and L1636.Additionally, hydrogen bonds are observed between the D15 sidechain carboxy group and the Y1542 hydroxy group, and between the Q18 sidechain carbonyl and the S1543 hydroxy group.Starting from the WT peptide a set of 68-point mutations, aimed to improve binding affinity (DG bind ) and complex thermodynamic stability (DG tot ), were evaluated according to Table 4.
The results of the RS analysis, summarized in Table 5, show that the T16W mutation positively affects both the binding affinity (DG bind ) and the thermodynamic stability (DG tot ) by À7.9 and À0.1 kcal mol À1 , respectively.A T16F mutation is also beneficial for both DG bind (À3.9 kcal mol À1 ) and DG tot (À0.2 kcal mol À1 ), while the T16Y slightly deteriorates the predicted total stability (DG tot = 0.1 kcal mol À1 ), but not the binding affinity (DG bind = À5.9 kcal mol À1 ).
Mutations at S21 (S21W and S21R) also provide improved affinity, although at the expense of overall stability (DG tot = 0.4 kcal mol À1 ).In particular, the highest decrease in binding affinity (coupled with the highest loss in thermodynamic stability, within the selected peptides) is observed for the S21W mutation (DG bind = À10.8kcal mol À1 ).All the peptides described herein were selected for experimental evaluation.
Production and characterization of FITC-labelled, in silico identified IGF2R peptide ligands The previously described dodecamer (P431) identified by Chen et al. has been selected as a positive reference for the binding affinity to the IGF2R in the screening of the newly in silico identified IGF2R ligands. 20Randomly generated scrambled sequences of all peptides were implemented into the binding assay for a comprehensive evaluation and to identify unspecific binding to the cell membrane. 33All peptides were labelled with FITC (Table 1), and in order to minimize the intramolecular hydrophobic interactions between the fluorophore and the amino acids of the peptides essential for IGF2R binding, a 3-glycine spacer was included in the sequence at the C-terminus. 34A reactive primary amine was introduced in the peptides by including a terminal lysine at the C-terminus, which was then used to attach FITC.The N-terminus was blocked by acetylation to prevent cyclization reactions, while N-methyl amine was used as a protection group for the C-terminus.Amino acid sequences and masses of all peptide ligands are reported in Table 1.A purity of Z94% was confirmed for all ligands.

Serum stability
When aiming at identifying novel targeting ligands to functionalize nanoparticles or biopolymers for diagnostic or theragnostic purposes, testing the serum stability of the targeting sequence is one of the pivotal parameters to investigate.To simulate the biological environment, non-FITC-labelled peptides were incubated in 50% v/v FBS in PBS and the amount of intact peptide was analyzed over time using analytical HPLC.A minimum proteolytic stability of 80% intact after 8 h and 60% intact after 24 h was set as the satisfactory threshold in this experiment. 35,36n contrast to the previously described P431 for which only traces (approx.3%) could be detected after 8 h of incubation in the serum-containing medium, the six in silico designed peptides (Tc, T1-T5) showed almost no proteolytic degradation after 8 h (intact peptide 485%; Fig. 3).After 24 h, P431 was completely degraded, while Z60% remained intact for most of the other candidates (except T3: 58.2 AE 2.7%), indicating their potential as targeting ligands for in vivo investigations.

Cellular metabolism assay
The in silico derived sequences were investigated for their role as targeting ligands for HSCs.Optimally, a targeting ligand for Fig. 2 Representation of the WT peptide in complex with the IGF2 receptor.Hot-spot and warm-spot residues, identified by CAS analysis, are coloured in red and orange, respectively.optical imaging should show a low intrinsic activity (competitive agonism) and should not induce apoptosis or alterations in metabolic activity. 37Therefore, the mitochondrial dehydrogenase activity was analyzed as a measure of metabolic activity after incubation with the peptide (5 and 10 mM) over a period of 8 h to reflect an expected timeframe of the receptor-ligand interaction in vivo.Untreated cells were used as the 100% activity control.All other results were calculated as the percentage of metabolic activity in relation to untreated cells.In this experiment, values Z70% were considered as indicating unchanged or slightly changed metabolic activity (Fig. 4).No remarkable effects on metabolic activity could be detected for the concentrations of 5 and 10 mM for all sequences.The designed peptides are suitable targeting ligands that do not induce substantial cytotoxicity.

Cellular uptake by the hepatic stellate cell line LX-2
In order to evaluate the potential of our peptide sequences as targeting ligands to be anchored to the surface of diagnostic liposomes, nanoparticles, or gadolinium complexes, it was of crucial importance to evaluate their binding efficiency to different cell lines.Accordingly, naı ¨ve LX-2 cells were used as a model of activated HSCs that highly express IGF2R (vide infra).In healthy liver tissue, HSCs are surrounded by hepatocytes forming liver sinusoids.Hence, HepG2 cells, a model for polarized human hepatocytes, were used to demonstrate superior peptide binding on LX-2 cells compared to other liver tissue.HEK293, immortalized human embryonic kidney cells, were used as a negative control.Cells were incubated for 1 h with the fluorescent peptides (10 mM) or their scrambled sequences (10 mM) and further analysed via flow cytometry.Percentages of FITC positive labelled cells (%FITC+) were obtained for LX-2 cells (Fig. 5A), HepG2 and HEK293 (data shown in the ESI, † Fig. S4B and C, respectively).Furthermore, the binding was performed by three different observers in two different labs (Fig. S5, ESI †).In order to display the superiority of the binding peptides to their scrambled sequences on LX-2 cells (Fig. 5B) and their intercellular selectivity compared to HepG2 (Fig. 5C) or HEK293 (Fig. 5D), the flow cytometry data were plotted as the ratio R SCR between the percentage of sequencespecific %FITC+ LX-2 cells and the percentage of sequence-specific %FITC+ of the scrambled peptides, HepG2 or HEK293.
Furthermore, to compare the strength of the bimolecular interaction with LX-2-cells expressing IGF2R, the equilibrium dissociation constant (K D ) was determined for all sequences.B max [mM] is defined as the maximum specific binding extrapolated to higher concentrations (4C max ) than those tested in the experiment.To determine K D , the %FITC+ living LX-2 cells were tested after incubation with different concentrations, i.e. 1, 5, 10, 25, and 50 mM, of the peptide.The results were calculated by means of nonlinear regression of the specific binding with Hill slope calculation as displayed in Fig. 6.The obtained K D , B max , coefficient of determination (R 2 ), and standard error of regression (SxÁy) values are shown in Table 6.
In contrast to cells treated with vehicle control (DMSO), all cells incubated with FITC-modified peptides generated a detectable fluorescence signal higher than the auto-fluorescence of untreated cells.The sequence Tc_FITC, a decamer representing the putative natural peptide sequence containing the most relevant amino acids for the interaction of IGF2 with IGF2R, was used as the reference peptide.Favourably, by replacing one amino acid in the sequence Tc, the binding affinity towards the target should be increased without impairing the natural binding to IGF2R.The selectivity toward LX-2 cells compared to HepG2 or HEK293 is expressed as R HepG2 or R HEK293 and should be 4 1 (Fig. 5C and D, respectively).As shown in Fig. 5 and 6, the reference ligand Tc_FITC showed the lowest binding affinity with K D = 13.03 mM on LX-2 cells (46.2% positively labelled cells (%FITC+)) while demonstrating the highest selectivity toward LX-2 cells relative to HepG2 and HEK293 (R HepG2 = 3.79, R HEK293 = 3.97).The highest peptide uptake was observed after treatment with the previously described sequence P431_FITC (%FITC+ = 82.8%).Interestingly, also its scrambled sequence SCR431_FITC, a negative control not included in the original publication but generated in this study to validate the binding properties of this allegedly IGF2R-specific peptide, bound in a similar fashion to LX-2 cells (%FITC+ = 85.2%,R SCR = 0.97, Fig. 5B), suggesting ultimately the lack of specificity of P431 for these cells.The binding affinity of P431_FITC taken from the original publication (K D = 6.19 mM) could be confirmed (K D = 6.26 mM) in this set of experiments, proving the validity of our testing method. 20However, for both sequences P431_FITC and its scrambled version SCR431_FITC, no specific selectivity for HSCs could be detected as the R HepG2 (0.94 and 0.92, respectively) and R HEK293 (0.89 and 0.91, respectively) values prove (Fig. 5C and D).Contrariwise, the in silico generated target-specific peptides were characterized by a 1.2-to 1.9-fold increase in binding compared to their corresponding scrambled sequences (R SCR , Fig. 5B) which was demonstrated as significantly effective in vivo. 38Thus, the increased binding of our peptides compared to their equivalent scrambled sequence indicates their potential as IGF2R-selective targeting ligands.
The sequence T1 was generated by mutating the hydrophilic amino acid threonine (log P = À0.26) at position 5 replacing its hydroxyl side group with an aromatic indole-functionalized amino acid such as tryptophan (log P = 1.46) in the middle of the sequence.A similar mutation was performed for peptide T3.At position 10 in sequence T3, the hydrophilic serine (log P = À0.45) was replaced with a tryptophan.The resulting increase in binding and decrease in selectivity were nearly identical to those of T1.Replacing the tyrosine at position 5 with a threonine (log P = 0.55) in sequence T2 (K D = 13.15 mM) did not substantially improve the binding affinity.However, T2 (R HepG2 = 2.73, R HEK293 = 2.41) is characterized by a noteworthy intercellular selectivity compared to the other mutated sequences.
The hydroxyphenyl side group of tyrosine is more lipophilic than the hydroxyl group in threonine but may also form hydrogen bonds which ultimately leads to only a small change in molecular charge. 39It can be hypothesized that the replacement of the non-aromatic hydrogen bond-forming amino acid with a lipophilic hydroxyl group-free amino acid at position 5 or 10 of the sequence led to an increased IGF2R affinity but to decreased selectivity, underlying the primary role of a hydrogen bond at position 5 to achieve pronounced selective binding.Replacing threonine with phenylalanine (log P = 1.16) at position 5 in sequence T4 resulted in a considerable 68% increase of the uptake.Consistent with the results obtained for peptides T1 and T3, a mutation towards a more lipophilic amino acid without being able to form hydrogen bonds resulted in higher binding and lower selectivity of the peptide.Besides the loss of hydrogen bonds, a decreased steric interaction associated with minor changes in the 3D structure compared to sequence T1 could be a pivotal parameter for the observed differences.An exception to the previously described hypothesis is represented by peptide T5, in which the amino acid serine at the extremity of the sequence was replaced by arginine (log P = À0.79).In this case, however, the mutation introduces a positive charge of the side chain, resulting in an overall neutral charge of the sequence at physiological pH, likely responsible for unspecific interactions with the cell membrane, as confirmed by the lower selectivity of peptide T5 (R HepG2 = 1.53,R HEK293 = 1.24). 40This could also be an explanation for the high uptake of FITC_P431 and its scrambled version in comparison to the other sequences.Both sequences (P431 and SCR431) consist of the amino acid arginine and histidine, ultimately resulting in a positive net charge at physiological pH leading to R SCR = 1 due to charge-related interactions between the cell membrane and the peptide. 37 negative net charge and a mutation towards a more lipophilic amino acid not able to form hydrogen bonds at position 5 or 10 of the sequence resulted overall in the best binding with respect to activated HSCs.However, these mutations negatively impacted cell selectivity.To generate a peptide with improved binding and selectivity, additional mutations need to be performed on the tested sequences.
To further confirm the binding affinity observed with FACS and to analyse the cellular location of the FITC-labelled peptides after incubation with cells, confocal microscopy imaging was carried out.LX-2 cells were incubated with 10 mM of FITClabelled peptides or their corresponding scrambled sequences for 1 h.To determine the absolute number of observed cells, the nucleus was stained with DAPI.In order to quantify the fluorescence signal of FITC, the corrected total cell fluorescence (CTCF) was calculated.As displayed in Fig. 7A and B, the majority of the FITC+ signal after 1 h is mainly located in the cellular lumen.Representative confocal microscopy images of LX-2 cells after the incubation with the scrambled sequences are available in the ESI † (Fig. S6A and B).This result confirms the ability of the receptor to rapidly internalize the specific ligand.Furthermore, the analysis of the CTCF (Fig. S7, ESI †) confirmed the results obtained by flow cytometry.The obtained pattern of the fluorescence intensity for all peptides (Fig. S7, ESI †) is comparable to the flow cytometry data displaying the %FITC+ events (Fig. 5A).The ability of being rapidly internalized makes our peptides promising candidates to decorate theragnostic

Cellular uptake in LX-2 cells in different states of activation
The progression and resolution of hepatic fibrosis is characterized as a complex dynamic process of interactions between parenchymal and non-parenchymal liver cells.HSCs, represented as a quiescentlike phenotype in a non-fibrotic liver, undergo a (reversible) transition towards an activated proliferative myofibroblast-like phenotype after being exposed to chronic apoptosis of hepatocytes associated with immune cell infiltration and the resulting release of proinflammatory and pro-fibrogenic cytokines. 41Furthermore, activated HSCs are characterized by attenuated proliferation and aberrant production of proteins, such as collagen type I and type IV, desmin, vimentin, and a-smooth muscle actin that form the major components of the ECM. 42articularly relevant for the present investigation are the different expression levels of IGF2R as a function of the activation state, with a higher amount of IGF2R expressed on transdifferentiated HSCs in comparison with the ones in a quiescent state. 43LX-2 cells retain crucial features of activated HSCs, but it is possible to perpetuate their activation through supplementary treatment with TGF-b1. 30On the contrary, by exposing LX-2 cells to a mixture of retinol and palmitic acid (Rol + PA), their quiescent-like phenotype can be promoted. 44n line with the work published by Chen et al., 20 we also shared the final goal to identify a ligand specific for activated HSCs.Aiming at following the progression and/or regression of hepatic fibrosis, it was important to prove that the uptake of the identified peptide sequences would be higher in naive and perpetuated LX-2 cells in comparison to quiescent-like LX-2 cells.The results of the ratio R fibrogenic (%FITC+ perpetuated LX-2/ %FITC+ quiescent LX-2) are displayed in Fig. 8.
A R fibrogenic value 4 1 is defined as increased selectivity towards the activated state.R fibrogenic r 1 characterizes no (=1) or reciprocal selectivity (o1) which is not desired for a future approach regarding the diagnosis of hepatic fibrosis.The treatment-specific %FITC+ for the quiescent-like and perpetuated LX-2-phenotypes are displayed in the ESI † in Fig. S8.Upon treatment of naı ¨ve LX-2 with TGF-b1, an increase of peptide uptake was observed for all the tested in silico derived peptides with respect to untreated LX-2 cells (Fig. S8 and Fig. 5A).As expected, after the treatment of naı ¨ve LX-2 with Rol + PA, a reduced binding could be detected for all in silico sequences, except T5 and P431.For P431, we obtained a R fibrogenic value of 0.99, indicating the same affinity for all states of activation and thus a lack of selectivity.Tc (R fibrogenic = 2.21), T1 (R fibrogenic = 1.58), and T2 (R fibrogenic = 1.53) showed a significantly decreased uptake after cells were treated with Rol + PA and an increased uptake after the treatment with TGF-b1 (p r 0.0001), making them the most promising sequences for a phenotype-specific uptake into HSCs.Sequence T4 (R fibrogenic = 1.25) was only characterized by a moderate difference towards the quiescent-like cells after perpetuation.T3 (R fibrogenic = 1.18) showed a significant but minor selective uptake to the perpetuated state.
For T5 (R fibrogenic = 0.95), no significant difference in uptake in naı ¨ve, perpetuated, or quiescent-like LX-2 was detected.This supports the hypothesis that T5 and P431 mainly or fully have unspecific interactions with the cellular membrane due to the presence of multiple arginines in the sequence, which elevate the positive charge of the molecule at physiological pH.Although some of our mutated peptides show selectivity towards the activated LX-2 cells, when compared to the putative binding ligand Tc, the implemented amino acid mutations seem to negatively impact (slightly for T1 and T2, and moderately for T3 and T4) selectivity.

Detection of IGF2R by a quantitative mass spectrometric approach
In order to provide experimental evidence of the IGF2R expression levels in LX-2 cells at different states of activation and to explain the results observed with peptide P431 proposed by Chen et al., 20 the expression of the receptor was analysed by a quantitative mass spectrometric approach.To the best of our knowledge, no information regarding the exact protein expression levels in HSCs and their comparison between the different activation-based phenotypes is available in the literature.In Fig. 9, the relative receptor expression of different phenotypes normalized to the expression of naive LX-2 cells is displayed.As indicated by the RNA-expression analysis performed by Beljaars et al., an activation-specific expression pattern of the IGF2Rpeptide could be observed compared to untreated LX-2 cells. 45,46redictably, compared to naive LX-2 cells, the perpetuated phenotype is characterized by an elevated receptor expression of approximately 16.8% (Fig. 9).
The quiescent-like phenotype is defined by a decrease of approximately 19.9%.Thus, an activation-specific targeting of our peptides would demonstrate specificity not only for HSCs with respect to other cell types but also for their state of activation, and therefore determine areas with higher disease progression (activated HSCs).In LX-2 cells, IGF2R is not exclusively expressed on the cell surface, but compared to our negative controls, HepG2 and HEK293, a larger proportion of the protein should be found. 45,46Therefore, we investigated the IGF2R-expression levels of our negative controls HepG2 and HEK293.These cells were chosen based on their relatively low IGF2R RNA-expression levels derived from the Human Protein Atlas. 47However, in an apparent contradiction to the assumption, the relative protein expression of IGF2R in HepG2-cells was 7.5-fold elevated compared to naı ¨ve LX-2 cells.A possible explanation can be found in the specific cell line and in the main localization of the receptor expression.Hepatocytes are characterized by a generally increased metabolization rate positively affecting the protein translation and their release in the endoplasmic reticulum. 48Furthermore, HepG2-cells are used as a model for polarized hepatocytes but represent a well-differentiated hepatocellular carcinoma. 49It has been described that protein and gene expression is increased in an uncontrolled and unique manner with respect to cultured primary human hepatocytes. 50herefore, it can be hypothesized that the majority of IGF2R is not Fig. 8 Ratio R fibrogenic (%FITC+ perpetuated LX-2 /%FITC+ quiescent LX-2 ) after incubation of perpetuated and quiescent-like LX-2 cells with FITC-labelled peptides at a concentration of 10 mM for 1 h at 37 1C.LX-2 cells were transdifferentiated into the perpetuated state with 10 ng mL À1 TGF-b1 or into the quiescent-like state with Rol (10 mM) and PA (300 mM) for 24 h at 37 1C prior to incubation with peptides.R fibrogenic 4 1 indicates selectivity for the perpetuated phenotype over quiescent-like LX-2 cells.In Table S3 (ESI †), the results of the multiple comparison tests (Holm-Sidak) between all peptides are displayed.expressed on the surface of the cell but in the Golgi apparatus and vesicles located in the cellular lumen. 51Thus, only a small proportion of IGF2R is located on the surface of HepG2 cells.The same explanation can be applied for HEK293 showing a 10% increase in the relative IGF2R-peptide expression with respect to naive LX-2 cells. 52

Conclusions
In this study, we used an efficient in silico approach for the development of targeting ligands selective for IGF2R.Structural information of the binding site of the receptor was used for the screening of specific peptide ligands.Molecular dynamics simulations, followed by computational alanine scanning of the IGF2R/IGF2 complex, led to the identification of a putative peptide sequence (Tc) containing the most relevant amino acids for the receptor-ligand interaction (IGF2 E12-C21).The Residue Scan tool, implemented in MOE software, was then used to optimize the binding affinity by replacing amino acids at position 5 or 10 of the putative binding sequence in order to implement mutations that increase affinity without affecting proteolytic stability.Selective binding of the peptides to IGF2R-expressing cells was successfully verified for all sequences (K D r 13.2 mM).With respect to Tc, the majority of mutations led to an increased receptor affinity (T1, T3, T4, T5), except for T2 where the affinity remained unchanged.In apparent contrast to the previously published P431, all sequences have shown a 2-to 4-fold increase in binding for LX-2 cells compared to the negative controls, HepG2 or HEK293 cells.They also showed superior binding compared to their associated scrambled sequences.All unlabelled sequences were tested in vitro for their proteolytic resistance by HPLC-DAD and were characterized by a stability of Z60% intact after 24 h incubation in 50% FBS (except T3, with 58%).In view of their prospective diagnostic application, a comparison of binding affinity was performed in perpetuated and quiescent-like LX-2 cells.Increased binding was detected in the perpetuated cell state, indicating selectivity for the diagnostically relevant phenotype.In summary, the increased binding affinity of the identified peptides towards perpetuated LX-2 cells compared to other LX-2 states or other cell types, and the satisfactory proteolytic stability show that the in silico designed sequences offer a new potential strategy for the targeting of hepatic fibrosis.

Fig. 3
Fig. 3 Proteolytic stability of the different peptides (25 mM) over 8 and 24 h in the presence of 50% v/v FBS in PBS at 37 1C.

Fig. 4
Fig.4Cell metabolic activity in % after an 8 h incubation of LX-2 cells with purified peptides at a concentration of 5 mM (A) or 10 mM (B) in complete growth medium at 37 1C.

Fig. 6
Fig. 6 Determination of the equilibrium dissociation constants (K D ) of the FITC-labelled peptide sequences after incubation with LX-2 cells.The population of %FITC+ for the different peptide concentrations was acquired after 1 h incubation at 37 1C.

Table 2
Interaction energy and specific contributions for the IGF2/IGF2R complex from MMGBSA analysis

Table 3
Changes in interaction energy predicted by alanine scanning

Table 4
Point mutations evaluated by the Residue Scan (RS) analysis, starting from the E12-S21 WT peptide

Table 5
Peptides selected according to the RS analysis.The differences in binding affinity (DG bind ) and complex stability (DG tot ) following the mutation are reported in kcal mol À1 .Only the results with DG bind o À3.5 and DG tot o 0.5 kcal mol À1 are reported Mutation DG bind [kcal mol À1 ] DG tot [kcal mol À1 ] Final sequence

Table 6
Equilibrium dissociation constant (K D ), calculated receptor saturation concentration (B max ) after incubation with LX-2 cells, correlation coefficient (R 2 ) and standard error of regression (SyÁx) of the correlation function carriers designed for early detection of alterations in HSCs and treatment of hepatic fibrosis.