Impact of different spacers on the conjugation between Anderson-Evans polyoxometalates and peptides

The Anderson-Evans polyoxometalates (POM) display a promising anticancer activity. The conjugation with the GRP-receptor antagonist peptide Demobesin (fQWAVGHL-NHEt) was exploited to impart cell targeting capabilities and improve the selectivity of such polyanions. However, the POM interacts with the grafted peptides, inducing chains folding and self-assembly of the resulting hybrids, thus decreasing their recognition ability. Within this context, a tailored spacer, including two domains, i.e., a hydrophilic one (1,13-diamino-4,7,10-trioxatridecan-succinamic acid, Ttds) and a tetra-anionic one (Glu-Glu-Glu-Glu-β Ala, EEEE-β A) was previously utilized to mitigate such interaction. In this work, hybrid POMs containing only Ttds or EEEE-β A were prepared and the contribution of the two spacers was separatelystudiedbyusing 2DNMR, ﬂ uorimetry and circular dichroism (CD). Transmission electron microscopy (TEM) was also used to observe the impact of the different spacers on self-assembly. Owing to the relevant effects observed for EEEE-β A, MD calculations were ﬁ nally performed to elucidate its behavior when incorporated in the hybrid POM. Our results show that, despite the stronger impact of EEEE-β A spacer, only when both spacer are present together it is possible to observe a signi ﬁ cant effect on the retention of peptide's secondary structure and recognition capability.


Introduction
Polyoxometalates (POMs), a class of polynuclear oxo-bridged transition metal complexes (Blazevic and Rompel, 2016), have received extensive attention due to their rich topology and tunable chemical and physical properties.In addition to their application in multidisciplinary fields (Omwoma et al., 2015), such as material science (Yu et al., 2023) and catalysis (Wang and Yang, 2015), their biomedical activity (Čolović et al., 2020), as antitumor (Bijelic et al., 2019), antiviral (Dan et al., 2020), and antibacterial (Bijelic et al., 2018) agents, has been highlighted.The main advantage of POMs is that their properties (especially dimensions and charge) can be tailored to optimize the interaction with biological macromolecules (Arefian et al., 2017;Lentink et al., 2023).Moreover, they can be rationally synthesized on a multi-gram scale at a low cost.
POMs were generally proven to cross cell membranes, although with high toxicity and instability under physiological pH conditions (Wu and Liang, 2017).To overcome these drawbacks, different delivery systems were devised (Bijelic et al., 2019), while the POM architectures were successfully modified by introducing transition metals (such as Co, Ru or Mn) (Wang et al., 2014;Fabbian et al., 2022;Carvalho and Aureliano, 2023) or by grafting suitable organic pendants (Cameron et al., 2022).This last approach, consisting of the hybridization of inorganic POMs with organic or biological moieties, by exploiting weak interactions and/or covalent bonds, increases both stability and bioavailability (Dolbecq et al., 2010;Chang et al., 2022), and provides interesting opportunities for tracking (Modugno et al., 2018) and targeting (Zamolo et al., 2018;Lentink et al., 2023).
In this work, the Mn-Anderson-Evans was functionalized with two 1 analogs containing one of the two portions of the previously described spacer.Peptide 2 contains only the Ttds group, while peptide 3 only the EEEE-βA anionic domain (Figure 1).The corresponding POM constructs (2-POM and 3-POM) were prepared (Figure 1 and Supplementary Figure S1) and compared with spacer-free, and Ttds-EEEE-βA containing POMs (1-POM and 4-POM, respectively), to investigate the contribution of different spacers.POM hybrids were characterized by various spectroscopic techniques, including ESI-MS, FT-IR, 2D-NMR, UV-Vis, CD and fluorimetry, DLS, TEM, and tested in vitro on cancer cells.Owing to its interesting properties, 3-POM was also investigated by MD calculations to explain the observed behavior.

Instruments
Nuclear Magnetic Resonance (NMR) Spectroscopy was performed on a Bruker DMX-400 instrument (Bruker Corporation, Billerica, MA, USA) operating at 399.92 MHz for 1 H. Proton chemical shifts, in parts per million (ppm), are referred to the residual 1 H-DMSO solvent signal in DMSO-d 6 (δ = 2.49 ppm).About 4 mg of each sample were dissolved in 600 μL DMSO-d 6 for analysis.All NMR experiments have been carried out at temperature of 298 K. One-dimensional (1D) 1 H-NMR spectra were acquired using typically 32 scans with 32K data size.For the two-dimensional (2D) 1 H NMR spectra, including COSY (Correlation Spectroscopy), TOCSY (Total Correlation Spectroscopy), and ROESY (Rotating Frame Overhauser Effect Spectroscopy), homonuclear pulse programs of the standard Bruker library were used (cosygpmfqf for COSY experiments, mlevetgp for TOCSY experiments, roesyetgp for ROESY experiments).Typically, 512 experiments of 50 scans each were performed: relaxation delay T2 equal to 1; size 4K; 5,597 Hz spectral width in F2.Spectral processing was performed using the software TopSpin.
Infrared Spectra were collected by preparing KBr pellets, on a Nicolet 5700 FI-IR instrument (Thermo Fisher Scientific, Waltham, MA; USA).
Electrospray Ionization Mass Spectroscopy (ESI-MS) was performed on a LC/MSD Trap SL Agilent instrument (Agilent Technologies, Santa Clara, CA, USA).The samples were dissolved in CH 3 CN and analyzed in negative mode.
Circular dichroism and UV-vis spectra were acquired on a Jasco J-1500 CD spectrometer combined with a Jasco PTC-423S temperature controller (JASCO Corporation, Tokyo, Japan).For Circular Dichroism (CD) analysis, 1 mg/mL stock solutions were prepared in TFE and diluted to 0.1 mg/mL in the investigated TFE/ water mixtures (%TFE= 10,20,40,60,80,100), which were prepared in 2 mL flasks.Far-UV CD spectra were collected (at least 16 scans per sample) in the 190-260 nm range in 0.1 cm pathlength quartz cuvettes (Hellma GmbH & Co., Müllheim, Germany) at 25 °C at 50 nm/min scanning speed, 1 s response time, 1 nm bandwidth, 200 mdeg sensitivity, 0.5 nm datapitch.The secondary structure estimation (SSE) analysis of CD spectra was performed both using the Jasco Spectra manager application and the CDApps software developed at B23 beamline of the Diamond Light Source synchrotron (Didcot, UK) (Hussain et al., 2015), enabling signal deconvolution by means of the CONTINLL algorithm (encompassing α-helix, distorted α -helix, β-strand, distorted β-strand, turns and unordered reference structures).
Emission fluorescence experiments were collected on a Perkin Elmer LS50B spectrofluorimeter, equipped with FL-WinLab software (PerkinElmer, Inc., Waltham, MA, USA).Samples were measured in 10 mm pathlength quartz cuvettes at 25 °C.Emission spectra were recorded in the range 300-550 nm, with excitation at λ=295 nm, 3.5 nm excitation and emission slits and 400 nm/min scan speed.The spectra were the average of 4 scans.
Dynamic Light Scattering (DLS) was monitored by a Zetasizer Nano ZS instrument (Malvern Panalytical Ltd, UK), collecting 3 series of measurements from 0.2 mM solutions in phosphate buffer (10 mM, pH 7) with 3% DMSO.
Standards S.r.l., Sesto San Giovanni, Italy), cultured according to the manufacturer's instructions, were used.On day 1, 10,000 HeLa cervical cancer cells/well in a volume of 200 µL were plated in 96wells plates, as previously reported (Tagliavini et al., 2021).The cells were exposed for 24, and 48 h to different concentrations of POM-TRIS, 2-POM, 3-POM (5, 25, 50, 100, and 150 μM).Stock solutions with 100 mM POM in water, containing 5% DMSO, were prepared by prior solubilization of POMs in DMSO, followed by dilution in water.The results were expressed as the mean ± SD percentage of living cells normalized to the control.

Molecular Dynamics (MD) investigations
The coordinates of the POM moiety, that was kept frozen during the simulations, were taken from the published crystal structure (Rosnes et al., 2012).For the peptide moieties, a fully extended starting structure was used.The peptide chains and the POM peptide complexes were build using standard molecular manipulation tools.For the peptides, the GROMOS 54a7 force field was used (Schmid et al., 2011).The parameters for non-standard residues were obtained by analogy with similar already parameterized chemical groups.As the POM structure is kept fixed during the MD simulations, only non-bonded parameters were derived, i.e., partial atomic charges and Lennard-Jones parameters.The partial atomic charges were evaluated from quantum chemical calculations on the isolated POM moiety bound on both sides to an amino group (see Figure 1).RESP charges (Bayly et al., 1993) were calculated at the density functional theory (DFT) level (Parr and Yang, 1995) with the B3LYP functional (Becke, 1993).The atomic basis sets were as follows: (i) for the Mn and Mo atoms, the LANL2DZ effective core potential for the inner electrons and a double Gaussian basis set of (5S,5P,5D)/[3S,3P,2D] quality for the valence electrons were used (Hay and Wadt, 1985); (ii) for the hydrogen, carbon, oxygen and nitrogen atoms, a standard 6-31+G(d) Gaussian basis set was used (Krishnan et al., 2008).Quantum chemical calculations were performed with the Gaussian09 package (M. Jea Frisch et al., 2009).The Lennard-Jones parameters for the POM were taken from the universal force field (UFF).
The POM-peptide hybrids were solvated in a dodecahedral box, large enough to contain the solute (with the peptides in the extended conformation) and at least 1.0 nm of solvent on all sides.Water was modeled by the simple point charge (SPC) model (Berendsen et al., 1987).The LINCS algorithm (Hess et al., 1997) was used to constrain bond lengths and a time step of 2 fs for numerical integration of the equations of motion was used.The particle mesh Ewald method (Darden et al., 1993) was used for the calculation of the long-range interactions and a cut-off of 1.1 nm was used.After a solute optimization and a subsequent solvent relaxation, each system was gradually heated from 50 to 300 K using short MD simulations.Then, a short (100 ps) equilibration simulation was performed for the two systems in the NPT ensemble with the velocity rescaling temperature coupling (Bussi et al., 2007) to keep the temperature constant at 300 K and the Berendsen pressure coupling (Berendsen et al., 1984) to keep the pressure constant at 1 bar.Three productive independent 100 nslong MD simulations for each system were then carried out in the NVT ensemble at 300 K.For both compound 1-POM and compound 3-POM, the same starting structure was used for each of the three simulations, varying the initial velocities randomly generated conforming to a Maxwell velocity distribution at 300 K.The MD simulations were performed with the GROMACS (Van Der Spoel et al., 2005) software package.The secondary structure of the peptides was analyzed with the DSSP tool (Kabsch and Sander, 1983) of the GROMACS package.

Peptide and POM-peptide synthesis
Peptides 1-4 were prepared by solid phase synthesis adopting the Fmoc/HBTU strategy, according to the previously described procedure (Rubini et al., 2010;Ventura et al., 2018;Tagliavini et al., 2021).The N-ethylamide group was introduced ab initio using an ethyl-indole AM Resin.The peptides were purified by HPLC and characterized by ESI-MS, FT-IR and 2D-NMR (see Supplementary Figures S2-S11, showing the characterization of peptides 2 and 3).
The syntheses of the 1-4-POM hybrids were performed in air, at room temperature for 24 h, under continuous stirring (Supplementary Figure S1) (Ventura et al., 2018;Tagliavini et al., 2021).To this aim, POM-NHS and 2 equivs. of peptides (1-4) were dissolved in DMF, then DIPEA was added.The obtained solutions were placed in a closed vessel under diethyl ether vapors, for 24 h, to allow the precipitation of the products, which were finally washed with diethyl ether and dried under air, to achieve pale orange solids (yields >85%).
FT-IR and ESI-MS are quick and useful tools to assess the successful conjugation.FTIR show the replacement of the C=O signals of NHS group, around 1740 cm −1 , by a band at 1,650-1,685 cm −1 due to the amide C=O bond and by a peak ascribed to N-H bending at 1,535-1,543 cm −1 .The integrity of the inorganic POM scaffold was confirmed by the retention of the bands at 942, 920-921, 900-902 and 663-670 cm −1 (Supplementary Figures S12, S20).ESI-MS (-) spectra show the bis-functionalized POMs, detected as di-or tri-anionic species (Supplementary Figures S13, S21).

NMR investigation
The resonance assignment of amino acid residues in peptides 2 and 3 and their corresponding POM derivatives were obtained by collecting COSY, TOCSY and ROESY spectra for all samples (Supplementary Tables S1-S4).Due to the limited solubility of POM hybrids in water, the analysis was performed in DMSO-d 6 .The latter is also useful for the study of peptides that interact with membrane receptors, having a dielectric constant (ε 47.5) very close to that of the membrane interface (ε 40).
Information on the secondary structure of peptides was obtained from the secondary chemical shift (Δδ) values of the H α , calculated by subtracting representative random coil values from the observed H α shift.In DMSO-d 6 , a minimum of four values of Δδ less than −0.16 ppm or greater than 0.22 ppm are indicative of the presence of α-helix or β-sheet secondary structures, respectively (Tremblay et al., 2010).According to these threshold values, it is possible to assess how the C-terminal regions of all peptides have a propensity to adopt a helical structure (see Val, Gly, His, Leu residues in Figure 2A), with the most negative Δδ values observed for peptide 2. Once conjugated on the POM, the peptides generally show even more negative Δδ values.The signals of peptide 1 show the biggest changes (see Δδ values for Gln and Trp residues before and after grafting, Figure 2), with a tendency to form the helical structure also at the N-terminus.2-POM, instead, displays its most negative Δδ value at the His residue, towards the C-terminus.
To evaluate the spacers' contribution on the interactions between POM and different peptides, the NH and Hα secondary chemical shift values, determined by subtracting the NH and H α resonance of each peptide to the resonance of corresponding POM construct were also compared (Figure 3).The presence of POM cluster had an evident impact in the resonance of the closer N-terminal peptide region in 1-POM and 2-POM, indicating that the introduction of the negatively charged Glu 4 -βAla sequence (in 3-POM and 4-POM) seems useful to attenuate the interaction between the peptide and the inorganic cluster.

Circular dichroism spectroscopy
The far-UV (185-260 nm) CD spectrum of peptide is largely determined by the electronic transitions of the amide chromophores and provides useful information on its secondary structure.To  overcome the high absorbance of DMSO in the far-UV region, that hinders its use as solvent in CD spectroscopy, 2,2,2-trifluoroethanol (TFE) was used as the solvent.TFE, in addition to being an excellent solvent for both peptide and the POM, is widely used to mimic the cell membrane environment (Ma et al., 2015).The influence of TFE on the secondary structure of peptides 2 and 3 was evaluated by recording CD spectra at increasing TFE percentages, so to monitor the behavior of the compounds in an environment with decreasing polarity (CD spectra in Figure 4; Supplementary Figures S28, S29, S32, S33), without the need of screening other solvents, where they display low solubility.
At a low percentage of TFE (10% v/v), the CD spectra of peptides 2 and 3 show a negative band at 197-200 nm (n → π* electronic transition) and a positive band at λ<190 nm (π → π*), ascribed to an unordered conformation.In 80% TFE, the CD spectra of peptides 2 and 3 are characterized by a positive band at 192 nm and a negative band around 206 nm, generated by the presence of more ordered structures, involving an increased contribution of the αhelical structure.
Parallel experiments on POM-peptide constructs showed that, at low percentages of TFE (10% v/v, Figure 4), the introduction of the POM moiety influences the position of the dichroic bands, whose negative minima display a shift to λ = 202-204 nm.As far as the bands intensity is concerned, the double value of ellipticity is due to the presence of two peptide chains in the hybrids.Moreover, 3-POM shows a stronger signal than 2-POM, (with Θ = 2.7×10 5 deg cm 2 dmol -1 at λ = 204 vs. 1.6×10 5 deg cm 2 dmol −1 at λ = 202 nm, respectively, see Supplementary Figures S15, S23), likely due also to the presence of additional aminoacids in 3-POM.In such condition, β-strand and random coil are the main conformations (37%-38% and 32%-35%, respectively).At higher TFE percentage (≥80% v/v), the bands further increase their intensity, especially the positive one at 191 nm (Θ = 5.4 × 10 5 and 7.5 × 10 5 deg cm 2 dmol -1 for 2-POM and 3-POM, respectively, see Supplementary Figures S29, S33), highlighting a strong tendency of the peptide to adopt to α-helix (up to 59%-68%) when exposed to the less hydrophilic medium and conjugated with the POM.The behavior observed in non-aqueous system is in agreement with the NMR analysis in DMSO-d 6 , and it is the result of a double effect: on one hand, TFE acts as secondary structure stabilizer, primarily inducing α-helical conformation (Vincenzi et al., 2019), on the other hand, the conjugation with POM may also assist the evolution to α-helix, as previously reported by (Yvon et al., 2014) for a different POM-peptide conjugate dissolved in CH 3 CN.Considering the estimated ratio between α-helix and βstrands in different TFE/H 2 O mixtures, hybrid POMs appear much more sensitive than free peptides to solvent change, being 2-POM the one showing the highest amount of α-helix (Supplementary Figures S30, S31, S34-S36).Nevertheless, the differences in the intensity of the bands level off when the compounds are situated in TFE ≥80% v/v.
For a sake of comparison, 1-POM and 4-POM showed the highest and the lowest impact on DB-1 folding, favoring, respectively, α-helix and β-strand conformations (Ventura et al., 2018;Tagliavini et al., 2021).In summary, the two spacers, Ttds and EEEEβA, have a relatively similar effect in terms of secondary structure evolution in H 2 O/TFE, and their behavior is intermediate with respect to the spacer-free and doublespacer POMs.
Interestingly, when POM and peptide are not conjugated (see Supplementary Figure S37, reporting the example for peptide 3 in the presence of 0.5 equiv. of POM-succ) no significant difference was observed in the CD spectra of peptide with or without POM at both the examined TFE percentages (10 and 40% v/v).This suggests that non-covalent interactions between POM and peptides are of intramolecular type and that the increase of dichroic signal observed for the hybrids can be due to an induction of chirality to the POM core in a less polar solvent, where intramolecular hydrogen bonds are promoted.

Fluorescence spectroscopy
Additional information on the interaction between peptides and POM was obtained analyzing the fluorescence emission of the Trp residue, being its emission dependent on the chemical surroundings.Far-UV CD spectra of 2 and 2-POM (A), and of and 3 and 3-POM (B), at 10 and 80% (v/v) TFE percentages (indicated).The spectra were recorded by a Jasco J-1500 spectropolarimeter using a 0.1 cm pathlength Suprasil quartz cuvette, 1 nm bandwidth, 50 nm/min scan speed, 1 nm data pitch, 1 s d.i.t.
Comparing the emission wavelengths and the intensity of the Trp fluorescence in equimolar solutions of peptides alone or graphed to POM we did not observe a relevant shift of the maximum wavelength emission of Trp residue, suggesting that the POM did not alter the chemical surrounding of the Trp side-chain.Nevertheless, the emission intensity is always less intense in graphted peptides than in the corresponding free peptide, highlighting the occurrence of an energy transfer from the peptide to the POM (Zhang et al., 2007).Fluorescence quenching experiments were then performed adding cesium ions, to provide information on the exposure of the Trp side-chain.In Figure 5 the observed fluorescence (F), at λ= 295 nm, is reported, as F/F0, with respect to the initial one (F 0 ).Instead of observing a further quenching of the fluorescence, however, Cs + induced an increase in fluorescence emission, suggesting that the added ions disrupt the interaction between POM moiety and peptide, consequently making less efficient the energy transfer between the two domains (Tagliavini et al., 2021).The increase in fluorescent emission is more evident in 1-POM construct and decreases with increasing spacer length.2-POM shows intermediate behavior, being also sensitive to the addition of Cs + in terms of fluorescence recovery, while 3-POM and 2-POM appeared almost insensitive to Cs + .The behavior observed upon addition of KI has a trend in agreement with previous data, with fluorescence increasing for 1-POM>2-POM and fluorescence quenching being more effective for 4-POM than for 3-POM, that can only be quenched at higher KI concentrations (Supplementary Figure S38).In the same conditions, peptides 1-4 showed the expected quenching of the signal, with efficiency in the order 1>2>3=4, in agreement with the repulsion between negative charges of I − and Glu residues in 3 and 4.
These fluorescence experiments highlight a greater impact of the spacers in terms of peptide-POM interactions rather than POMinduced secondary structure evolution, demonstrating the need to attenuate them to increase peptide availability.

Transmission Electron Microscopy
The formation of aggregates may also have a strong impact on the accessibility of peptide chains.TEM was employed to achieve information about size and shape of the nanoparticles (Figure 6).In water with 3% DMSO, 1-POM forms particles with diameter up to 50 nm diameter (Tagliavini et al., 2021), 2-POM assembles into fibers of 200-300 nm in length, whereas 3-POM, in agreement with DLS data (showing particles with hydrodynamic diameter of 30-50 nm, see Supplementary Figure S39), forms amorphous particles with diameter smaller than 50 nm.Compound 4-POM shows entangled fibers in the 300-500 nm diameter range.To understand the role of POM, with respect to that of the peptides, in driving the morphology of the aggregates, we have also monitored the behavior of POM-TRIS, POM-succ and of a mixture of both.POM-TRIS (casted form acetonitrile solution) displays a tendency to align, forming fibers between 500 nm and 3 μm in length, as shown in Supplementary Figure S40, suggesting that the 1D morphology is likely dictated by POMs.On the contrary, POMsucc forms polyhedral crystals-like structures (Supplementary Figure S41).Moreover, fibers longer than 500 nm, differently from the ones formed by POM-TRIS, can be obtained by mixing equimolar amounts of POM-TRIS and POM-succ in acetonitrile, as represented in Supplementary Figure S42.Complementary interactions (amino groups/POM surface and amino groups/ carboxylic groups) are thus also useful to establish directional interactions.Highly negatively charged POMs, as POM-succ and 3-POM, instead, only form nanoparticles.
Noteworthy, in DMF as solvent, 2-POM showed a gelator effect, likely mediated by the polar Ttds chains (Supplementary Figure S43), as previously observed also for 4-POM (Tagliavini et al., 2021).These results underly how hybrid POMs may represent versatile building blocks, capable of exploiting both the polyanionic surface and the organic pendants to control the interactions.

MD simulations
From the data collected, it is evident that the spacer based on the tetra glutamic acid has a bigger impact than Ttds on improving the accessibility of the peptide chains.Indeed, POM-3 showed better behavior in terms of: • retention of POM-free peptide's secondary structure in DMSO-d 6 , • higher accessibility of the peptide chain with respect to quenching by ions, • lower tendency to aggregate.
MD simulations were thus performed to elucidate the structural features of POM-3, comparing the outcomes with those of 1-POM, in order to focus only on the effects of the addition of the anionic spacer.The structure and dynamics were studied by means of three independent 100 ns-long MD simulations for each compound, starting from an extended conformation for the peptide moiety.Such a starting conformation allows the secondary structure of both peptides to evolve in an unbiased way, highlighting possible differences.The analysis of the secondary structure sampled in Fluorescence emission of POM-peptide hybrids (indicated) at increasing amounts of cesium ions (CsCl).Fluorescence is referred to as the ratio between the fluorescence at a given quencher concentration (F) and the fluorescence in its absence (F 0 ).
the MD simulations is reported in Figure 7, where the time evolution of the fraction of β-structure, helical structure and unstructured peptide is shown for both compounds.It must be noted that a relevant amount of unstructured peptide is present in the MD simulations (≈50% on average), while in the experimental CD data a lower percentage of unstructured peptides was observed  Time evolution of the relative abundance of β-structure (red), helical structure (blue) and unstructured peptide (gray) along the three independent MD simulations of compound 1-POM (upper panels) and compound 3-POM (lower panels).In the right panels, the same fractions, averaged along the three simulations, are reported for both compounds.To calculate such average values, the first 10 ns of each MD simulation have been neglected.
(32-35%).This can be ascribed to the relatively short simulation time (100 ns for each MD) that does not allow the peptides to assume a stable folded conformation.However, relevant structural differences between the two compounds can be already observed: in both compounds the peptides assume a β-structure, but in compound 1-POM a non-negligible fraction of helical structure is also present.The formation of helical structure is observed in two of the three simulations of 1-POM: in the third simulation, a particularly high unstructured fraction is present, and no helical structure is observed.3-POM shows a bigger amount of βstructures.Similarly to the previous case, in the first simulation of compound 3-POM a more relevant percentage of helical structure is present with respect to the other two simulations of the same compound.The statistical sampling of different secondary structures is not surprising in MD simulations starting from an extended conformation.However, the difference in the average secondary structure content between the two compounds appears to be meaningful.
During the MD simulations, the peptide bent on the POM.In some cases, the POM-peptide interaction was stabilized by the formation of hydrogen bonds between the peptide and the oxygen atoms of the POM.An example of such bending is reported in Figure 8.The bending of the peptide on the POM moiety, monitored along the six MD simulations by plotting the distance between the peptide terminal residue and the POM center of mass, does not seem to be related to the structural differences.In fact, a similar distribution for the terminal residue to POM distance is obtained for both compounds (Supplementary Figure S44).This fact is important to explain the presence of interactions despite the presence of electrostatic repulsions and deserves attention, in terms or length of the peptide and functionality at the terminal end, for further structural optimization.
The interactions between the two peptide chains, belonging to the same POM, were finally investigated for the two compounds, revealing that, in compound 1-POM, more interactions between the two peptide chains are present compared to compound 3-POM.This can be observed from Figure 9, where the time fraction of interchain contacts along the simulations of the two compounds is reported.Such a fraction is defined as follows.Being R i and R j the i th and j th residues of the two peptide chains, respectively, a contact between R i and R j is present at each MD frame in which the minimum distance d between R i and R j is below 0.5 nm. Figure 10 clearly shows that the two peptide chains are more frequently interacting during the MD simulations of compound 1-POM.The increased interchain interactions between the two peptides disfavor the formation of intrachain interactions.As a matter of fact, it can be observed in Supplementary Figure S45 that the number of intrachain hydrogen bonds (HB) in 3-POM is higher with respect to that in 1-POM.Nonetheless, the same Figure also shows that the number of interchain HBs is similar for 1-POM and 3-POM, suggesting that hydrogen bonding is not the main interchain interaction.Inspection of the trajectories reveals that interchain contacts in 1-POM are stabilized by aromatic interactions among His, Phe and Trp residues, that are often arranged in a sort of aromatic cleft (Supplementary Figure S46).
The amount of such transient interactions is likely related to the different conformations observed; in particular, they can make the chains less available for longer range interactions required for β-structures as well as for targeting purposes.Time evolution of the distance between the terminal residue of the two peptides and the center of mass of the POM moiety in one of the MD simulations of compound 3-POM.The black and red lines refer to the two pendants.Some representative snapshots of the POM-peptides conformation are also shown.

Cytotoxicity tests
The toxicity assessments were performed on HeLa cervical cancer cells, which show a moderate overexpression of the bombesin receptor.While the peptides 2 and 3, as well as 2-POM exhibit only a minor anticancer activity, (IC 50 >200 μM), 3-POM led to 58% residual viability at 150 μM after 48 h (IC 50 180 μM at 48 h, Figure 10).With respect to POM-TRIS (IC 50 =134 μM, with 46% residual viability at 150 μM after 48 h) and 4-POM (IC 50 =95 μM, with 33% residual viability at 150 μM after 48 h), the 3-POM is, thus, less active.However, with respect to 1-POM, which was also less active than POM-TRIS (Ventura et al., 2018), it shows a better activity.Moreover, the stronger effect observed at 25 μM of 3-POM, with respect to POM-TRIS, speaks in favor of a higher selectivity, at least, in the low concentration range (Table 1).These outcomes are consistent with the data herein collected, confirming, on one hand, the significant impact of tetra glutamic acid domain, on the other hand, the need of combining the two spacers to achieve a stronger biological activity.

Conclusion
This study delves into the asymmetric environment of novel Mn-Anderson POM-hybrids with potential anticancer activity.Specifically, distinct spacers, EEEEβA and Ttds, were introduced between the POM core and the pendant peptides, with the aim to maintain the peptide chains easily accessible even after grafting, so to facilitate the targeting of cancer cells.Peptide folding is, indeed, a common phenomenon in the presence of salts or other chaotropic agents, including POMs, and could lead to a decrease of biological activity (Soria-Carrera et al., 2023).
The incorporation of spacers could mitigate the influence of the POM cluster.According to 2D NMR, among the two spacers, EEEEβA and Ttds, the first proved to be more useful in minimizing molecular interactions.The negatively charged spacer, indeed, was confirmed to have a positive effect, as 3-POM and 4-POM showed lower CH α shifts compared to the free peptide.In contrast, compounds 1-POM and 2-POM displayed more differences in the secondary structure compared to bare Demobesin-1, favouring the formation of α-helix in DMSO-d 6 .However, the impact of the two spacers on secondary structure seemed not relevant in TFA/H 2 O solutions: CD spectroscopy did not highlight Fraction of interchain contacts in compounds 1-POM (left) and 3-POM (right).The color code represents the fraction of MD frames at which an Ri Rj contact is present.Quenching of fluorescence, instead, underscored the positive effect conferred by the charged spacer, with compounds 1-POM and 2-POM exhibiting no dynamic quenching upon the addition of both KI and CsCl, while 3-POM and 4-POM showed better responsivity of the side chain Trp.MD calculation showed that 3-POM secondary structure is mostly composed by β-sheet and random coil conformations, with minor interchain interactions.Despite the high accessibility of the chains, however, they still tend to fold on the POM structure, exploiting the terminal residues for the interactions.For this reason, the biological activity of 3-POM towards HeLa cells appears still low, and justifies the need for a longer spacer, provided by the addition of Ttds, as for 4-POM.This observation may be crucial to overcome the general problem of undesired POM-triggered peptide folding, and will be useful to drive the design of a next generation of POM-based drugs, with possibility to control targeting and delivery, while enabling mechanistic studies for the still unclear mechanism of action of anticancer POMs.Further studied will evaluate the nature of POM scaffold.

FIGURE 2 Secondary
FIGURE 2 Secondary Hα chemical shift (Δδ) values of DB-1 derivatives (A) or of the corresponding POM-peptide hybrids (B).The dotted line indicates the threshold value for the presence of a helical secondary structure.

FIGURE 3
FIGURE 3 Influence of spacer on NH (A) and Hα (B) resonances of DB-1 in POM-peptide constructs.Secondary chemical shift values were determined by subtracting the NH or Hα chemical shift values of DB-1 analogues (peptide 1-4) to the resonance of corresponding POM constructs (POM 1-4).Values greater than |0.15| ppm (dotted lines) indicate a significant effect on resonances.

FIGURE 10
FIGURE 10Cytotoxicity performances of 3-POM on HeLa cells at increasing concentration.Results are presented after 24 and 48 h.

TABLE 1
IC 50 values ad residual viability, observed for HeLa cells treated with POMs, after 48 h incubation.Frontiers in Chemical Biology frontiersin.orga meaningful difference between 2-POM and 3-POM, which both displayed intermediate behavior with respect to 1-POM and 4-POM.